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FOREWORD 


The Society of Engineering Science, Huntsville Chapter, organized the Symposium, 
the proceedings of which are presented on the following pages. The purpose of 
the two-day meeting was to provide a forum for the presentation and discussion 
of various ideas in connection with the continuation and utilization of the United 
States Apollo space program. In fact, the Symposium had been conceived and 
planned prior to any official definition of aims or goals by the National 
Aeronautics and Space Administration. This is the reason for the great variety 
of problems discussed herein; the connecting thread between the various papers 
is the possibility of future application to what at first was known as Apollo ex- 
tensions, but now has the name of Apollo Applications. 

In the course of preparing for the Symposium, the indefatigable chairman of this 
event, Dr. Rolland G. Sturm, has suffered a great personal loss. This is why 
the following resolution was adopted by the Executive Committee of the Society 
of Engineering Science , Huntsville Chapter: "Be it resolved, that the Proceed- 
ings of this Symposium shall be dedicated to the memory of Mrs. Rolland G. 
Sturm". 
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WELCOME 


by 

Ervin Y. Rodin 

I welcome you to our Symposium on Apollo Applications; a meeting which marks 
at least two ’’firsts". The first step is the fact that what is now called Apollo 
Applications, but which probably will change its name many times in the future, 
is now going to be discussed in the form of a symposium for the first time. The 
second first is for the Huntsville Chapter of the Society of Engineering Science: 
this is our first symposium. I would like to say a few words about both of these 
things. I shall start with the one on which less money was, is being, and shall 
be spent than on Apollo Applications: Our Society. 

In the last few decades advances made both in engineering and research have 
reached such staggering proportions that it became entirely impossible for any 
one person to have an all encompassing knowledge of the kind that 19th century 
scientists used to have. Some of this is of course unavoidable; however, we can 
blame only ourselves for a significant portion of this inability for the assimila- 
tion of more knowledge. In particular, I am sure that all of us here have en- 
countered many times the basic communication and research gap between 
engineers and scientists. There are frustrating situations when the engineer 
cannot interpret the results of the research worker or when the research worker 
comes up with conditions and requirements unattainable by the engineer, or when 
a problem defies idealization and thus transcends and cuts across the boundaries 
of many disciplines. 

In order to meet such problems and to try to find answers for them, a group of 
eminent scientists and engineers from many different fields of endeavor got 
together and founded the Society of Engineering Science. This happened at 
Purdue University about four years ago. 



However, a few of us here in Huntsville felt that what the National Society had to 
offer was not sufficient; we wanted a more concentrated program, right on the 
local scene and level. So, just about a year ago we formed the first Professional 
Chapter of the Society. 

I believe that this first year of ours was a successful one. We had three technical 
meetings, addressed respectively by Dr. Edward Teller of the University of 
California, Dr. Rudolf Hermann of the Alabama Research Institute and Mr. Jim 
Farrior of Lockheed, and, finally, by Dr. Henry Eyring of the University of 
Utah. These men, eminent representatives of physics, aerodynamics, manage- 
ment and chemistry helped us in a large measure towards the attainment of 
better interdisciplinary cooperation. 

This symposium is essentially the culmination of our first year's efforts. We 
are certainly gratified by the response to it, as exhibited in our program; 
government, universities, and industry are all represented here. We take this 
as yet another indication that scientists and engineers welcome the opportunity 
to meet each other and exchange their ideas in a way so that the physicist will 
interest the space physician and the mining engineer will want to integrate his 
efforts with those of the mathematician. 

A final few words about Apollo Applications. What is the philosophy behind our 
symposium? It was, I think, summarized in advanced by what Dr. Teller had 
said during his visit here: 

The growth of technology is exponential. This in part is 
a correct statement based on experience. With some 
optimism one could even call it a mathematical statement 
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More knowledge produces more problems and more problems, 
more demands. It is reasonable to expect and it seems to be 
true that progress at any time is proportional to general 
development up to that time. But such simple description 
of the growth of the topics in which we are all interested 
is really and obviously an improper simplification. The 
growth is very much a structured growth. The field is 
becoming more complex and more specialized. It is also 
becoming more interrelated. 

Thus, an exchange of ideas, both structured and specific, becomes imperative. 

This Symposium, we hope, will be a forum for just such an exchange. We hope 
that all who come to present ideas, will have ample opportunity to discuss them. 
And we hope that those who are here to assess the state of these arts will also 
be able to do so. 

Thus, I now declare this Symposium to be in session. 
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SOME PROBLEMS IN DESIGNING A LUNAR PHOTOGRAPHIC SURVEY MISSION 


By 

W. L. Austin and J. L. Harden ' ' 

/General Electric Company 

INTRODUCTION 


Objective and Scope jy 

The objective of this paper is to discuss some of the problems associated with the de- 
sign of a lunar photographic survey mission. The specific problems to be discussed 
are the lighting conditions on the lunar surface and the translunar trajectory. 


Definitions and Limitations 


Of prime importance to this photographic mission is the illumination of the lunar sur- 
face. The photographs that can be produced by any given camera will vary in detail 
depending on the lunar surface shadows. Of least detail will be those lunar surface 
points that are near the lunar terminator because the angle of incidence is too small 
and those that are near the subsolar point because the angle of incidence is too large. 
The terminator is defined as a great circle along the lunar surface and normal to a line 
drawn from the center of the sun through the center of the moon. The subsolar point 
is the juncture of this line with the lunar surface and is within 1* 34 1 of the lunar equa- 
tor. Between these two extreme points of poor photographic detail is an area of opti- 
mum photographic detail. For this discussion, the area that yields optimum photo- 
graphic detail will be referred to as the preferred lighting ring. (See Figure 1.) 


The purpose in analyzing the preferred lighting ring is twofold: (1) The lighting condi- 
tions on the lunar surface that are required to obtain photographs of acceptable quality 
partly establish the lunar boundary conditions for orbital insertion, and (2) to show the 
relationship of the orbiting vehicle to the preferred lighting ring during a 28-day mis- 
sion in a polar orbit. 

The translunar trajectories to be discussed should not be accepted as general results 
since they are for a particular day in a particular synodic month. However, we may 
accept as a general result that once the lunar lighting conditions and lunar arrival long- 
itude are specified, there is at best one launch window in any given lunar qynodie month. 
(A synodic month is 28.68 days long and is measured from lull moon to full moon.) 
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Terminator 



. General Lunar Lighting Conditions 







MISSION DESCRIPTION 


Mission Objectives 

The photographic mission envisioned will have two functions, mapping and surveying. 
The purpose of the mapping camera is to photograph as much of the lunar surface as 
possible, irrespective of sun angles, as long as there is illumination, A number of 
preselected sites, survey sites, will be photographed at preferred sun angles by a 
high-resolution camera. 


Lunar Orbit Parameters 

An 80-nautical mile lunar circular orbit was chosen for analysis purposes, A polar 
orbit was chosen in order to maximize latitude coverage. There exists a definite trade- 
off between photographic quality, lunar orbit stay-time, and propulsion requirements. 
This is the reason why transearth trajectories were not considered since the time of 
transearth injection from lunar orbit insertion greatly affects the characteristics of the 
transearth trajectory. However, the lighting analysis is based on a 28-day lunar orbit 
duration in order to illustrate the maximum possible surface coverage. 

Arrival Longitude 

In order to maximize the inform? cion gained by photographing the far side in the early 
part of the mission, an approximate arrival longitude of 81 degrees west was chosen. 
This arrival permits several hours for checkout and orbit refinement before initiating 
photography at 90 degrees west. The Russians have photographed the far side; how- 
ever, qualitatively and quantitatively a great deal of information would be gained by 
re -photograph ing the area. Looking down from the lunar North Pole, Figure 2 depicts 
the equatorial longitudes that were photographed by the Russians. Moreover, the 81- 
degree west arrival longitude is near the minimum change of velocity (AV) required for 
translunar injection and deboost into lunar orbit (NASA Lunar Flight Handbook, Vol- 
ume 2, Part 2). This arrival longitude is a goal. How well it is met will be illustrated 
subsequently. Now that the mission objectives are generally defined, the next step is to 
derive the general equations for the preferred lighting ring and to solve them for an 
example mission. 
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Figure 2. Photographed Areas of the Lunar Surface 


LUNAR LIGHTING CONDITIONS 


Preferred Lighting Ring 

The preferred lighting ring was previously defined as the preferred lunar lighting con- 
ditions that are required by a survey camera to produce meaningful photographs. An- 
other way of defining it is the surface area of the moon within which the sun angle at 
every point in the area is within acceptable lighting limits. Thus, when a point in the 
preferred lighting ring is photographed, the resulting picture is of acceptable quality. 
The sun angle, S^, is defined as the elevation of the sun at a point on the lunar 
surface with respect to the local horizontal. For analysis purposes (see Figure 3) 
consider the ring to be bounded by two variables, a 1 and a 2 , which are lengths of 
arc across the lunar surface. o? 1 is the minimum and cn 2 the maximum angular dis- 
tances from the subsolar i oint allowed for photographic purposes. The point where 
= 90° is another definition of the subsolar point. This implies that the sun angles 
of and cv 2 are 90° -o^ and 90°-o; 2 , respectively. Thus, a site is within photographic 
limits when the site sun angle, S^, is: (90°-a 2 ) — ~ (90 # -arj). 


N 



7 = Latitude of site. 

X = Angular distance from 

the subsolar point to 
the site. 

L c = Longitude under con- 

c sideration (measured 

pos. west 0 - 360"). 

| y - LJ = Angular distance along 
the equator of the site 
longitude from the 
subsolar point. 


Figure 3. Geometry Required to Derive the General Equations 
for the Preferred Lighting Ring 
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Given any point on the lunar surface, the time of the synodic month and declination of 
the point with respect to the ecliptic determine the angle the sun makes with the local 
horizontal at the point on the surface. If this point is one of any of the predetermined 
survey sites, then the sun angle is critical. The general equations to determine the 
position of the preferred lighting ring with respect to the lunar surface can easily be 
found. Thus, for specific values of and a latitude coverage may be determined 
simply by substituting the a values for x in Equation 1 and solving for y. 

In deriving the general equations, a spherical moon is assumed with the moon’s equator 
in the plane of the ecliptic. Thus, the declination of the site with respect to the ecliptic 
is equal to the site latitude. Also, the distance from the center of the sun to any point 
on the illuminated side is treated as a constant. The initial time reference point is the 
beginning of the synodic month or full moon. 

The easterly sidereal rotation of the moon is 0. 549 degrees/hour. The terminator ro- 
tates easterly at the rate of 0.041 degrees/hour. If the effects of lunar librations are 
o muted, then the subsolar point, y, is located approximately (0.549 - 0. 041)X or 
0. 508X degrees west of the lunar prime meridian, where X is the time in hours from 
the last full moon. This approximation may be refined by using ephemeris data. 

In order to calculate the sun angle, S^, for any given latitude and longitude, Napier's 
rules apply (see Figure 3). 

Ey Napier's rules for right spherical triangles: 

Cos X = Cos y Cos | y - L q | 

X = Arc Cos [ Cos y Cos | 0. 508X - L» c | ] 

where, L is the longitude under consideration. 

Thus, for any given latitude, longitude, and time of synodic month, the of that lati- 
tude and longitude is given by: 

S A = 90° - X = 90* - Arc Cos [Cosy Cos |0.508X - L c 1 1 (3) 

In determining the boundaries of the preferred lighting ring with respect to the subsolar 
point, it is necessary to consider the two boundary points on X , a , and a g , where 


( 1 ) 

( 2 ) 
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0° < a x < a 2 < 90°. For the o? 2 value, 0° < | 0. 508X - L q | ^ a 2 and for the value, 
0° < | 0. 508X - L c j < cv^. The first inequality defines the outer boundary, and the 
second, the inner boundary of the preferred lighting ring. 


In the example mission which I will present, the preferred lighting ring is bounded by 
or = 30° and cx z = 80°. Thus, the photographs of a site with a sun angle less than or 
equal to 90°-o; i = 60° and greater than or equal to 90 °-q; 2 = 10° would be of acceptable 
quality. The rest of this section will be concerned with the development of the param- 
eters of the preferred lighting ring used for the example mission. Figure 4 is a gen- 
eral description of one quadrant of this preferred lighting ring which relates latitude 
coverage with the number of degrees, I y - L I , from the subsolar point. Note that 
the relation is symmetric about both axes, thus only one quadrant need be drawn. Also, 
the maximum latitude coverage is ±78.5 degrees when photography is initiated 30 de- 
grees west of the subsolar point. This point occurs at I y - L I = 30° which is readily 
seen in Figure 4. The latitudes in the figure were calculated in the following manner. 

In Equation 1 let X = 80°, hence: 


Cos 80° = Cos y Cos | y - L |=» 


(4) 


Y 


= Arc Cos 


0. 17364818 

Cos y - L ’ 
c 


where 0° < y 



(5) 


This results in a projection on the lunar surface of the 10-degree line which shows 
the outer limit of the preferred lighting ring. 


In Equation 1, let X = 30°; therefore 


Cos 30 


= Cos y Cos | y - L c | => 


(6) 


= Arc Cos 


0. 86602540 
Cos I y - L 


c ' -» 


, where 0° < | y - L | s 30° (7) 


This is a projection on the lunar surface of the 60-degree line which is the inner 
limit of the preferred lighting ring. Also, the 20-, 30-, 40-, and 50-degree lines 
were calculated and plotted in Figure 4. 


Polar Orbit and Latitude Coverage 

The vertical dotted lines in Figure 4 represent the relationship of the vehicle's orbit 
to the sun angle lines, and show the latitude coverage with respect to the sun angle lines 
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90° 80° 70° 60° 50° 40’ 30° 20° 10° 0°Subsolar 


iJegrees From Subsolar Point, \y - L Q | Point 


Figure 4. Relationship Between Orbiting Vehicle and an Orbital Insertion 
of 30° Wast of the Subsolar Point 

at the 0-, 10-, 20-, and 28-day time points in the orbit. The reason for the change is 
due to the rotation of the terminator. Table 1 is a summary of the latitude coverage at 
T = 0, 10, 20, and 28 days. Figure 5 shows the relationship between the Earth, Sun, 
and the Moon’s terminator for one synodic month. 

Table 1 

Latitude Coverage 


S A Line 
A 

H 

II 

o 

T = 10 

T = 20 

T = 28 

10° 

+78.5° 

±77.0° 

±73. 5* 

±71.0* 

e 

O 

CM 

±67. 0® 

±63. 5® 

±57.5° 

±51.0* 

30° 

±54. 5* 

±49.0* 

±39.0° 

±21.5° 

e 

O 

±42.0* 

±33.0° 

± 2.5* 

0.0* 

50° 

±28.0* 

± 5.0° 

0.0* 

o 

• 

o 

• 

Oi 

o 

o 

0.0* 

0.0° 

0.0° 

0.0° 
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Figure 5. Relationship Between the Earth, Moon, and Sun for Any Synodic Month 

Direction of Latitude Decay 

Referring to Figure 4, there is no apparent reason why the times stated for the dotted 
lines could not be reversed; i.e., the vehicle could arrive east or west of the subsolar 
point. If the vehicle arrived east of the subsolar point, then the equatorial S^ on the 
illuminated side would be 32. 5 degrees at arrival and 60 degrees at the end of the mis- 
sion. This increase in sun angle with respect to time is due to the easterly rotation of 
the terminator. The maximum latitude coverage at the beginning of the mission would 
be ±71 degrees and +78. 5 degrees at the end of the mission with a minimum sun angle 
of 10 degrees. If the vehicle arrived 30 degrees west of the subsolar point, then the 
equatorial would be 60 degrees at arrival on the illuminated side with +78. 5-degree 
latitude coverage. At the end of the mission, the equatorial would decay to 32. 5 
degrees with +71 degrees in latitude coverage at a minimum sun angle of 10 degrees. 
Thus, an arrival at 30 degrees west of the subsolar point results in a 15-degree greater 
latitude coverage north and south, as opposed to arrival at 57.5 degrees east of the 
subsolar point. Since it was assumed that the far side would be photographed first, the 
vehicle will be assumed to arrive 30 degrees west of the subsolar point in order to max- 
imize latitude coverage in the early part of the mission. This would be particularly de- 
sirable for a 14-day mission. Figure 6 shows the relationship of the vehicle to the 
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moon at arrival. Now that all of the lunar orbit parameters, including lighting, are 
established, the next step is to discuss the survey site coverage and to present an ex- 
ample mission. 




EXAMPLE MISSION 


Survey Site Coverage 

Figure 7 is an equigrid lunar map with the sun angle lines projected on it. These lines 
show the decrease in latitude coverage for the various sun angles over a 28-day mis- 
sion. The outer limits of the lunar area which fall within the photographic limits are 
defined by the 10-degree sun "ingle lines. This is due to the lighting constraints re- 
quired to obtain photographs of acceptable quality and due to the slope of the minimum 
permissible sun angle line, which in this case is 10 degrees. However, all of the lunar 
surface above and below this region is out of the acceptable photographic limits cOosen 
for the example mission. 

If a minimum latitude decay is desired, then a choice must be made between photo- 
graph quality and longitudinal change capability. For if the latitude coverage decay 
over a 28-day mission should be below some preselected value, longitudinal changes 
would be necessary at some point or points in the mission to prevent it. The nature 
of the change would be easterly in the direction of the subsolar point. Such a change 
would place the orbit closer to the subsolar point and then increase latitude coverage. 

If the changes are large and the propellant is not available for these changes for a 28-day 
mission, thenperhaps a shorter mission duration should be considered. Consequently, a 
trade-off exists between mission length, photographic quality, and permissible latitude 
coverage decay. Ideally, the best mode is to avoid making any longitudinal changes at 
all. However, it is interesting to note how much incremental change in velocity (AV) 
would be required per degree of longitudinal change for an 80-n. mi. circular orbit. 

This is illustrated in Figure 8. These values may not appear to be very large, but 
when they are added to all of the other AV requirements, their significance is bet- 
ter enhanced. 


Trajectory Discussion fo*' Juni 23, 1970 

The month of June, 1970, was arbitrarily selected to illustrate the problems involved 
in determining an acceptable translunar trajectory which meets prescribed boundary 
conditions at the moon. For this example mission, the boundary conditions at the 
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Figure 7. Lunar Surface Lighting for 28-Day Mission, 81° West Orbital Insertion 





Figure 8. Velocity Requirements for Longitudinal Changes, AL 

moon are a polar lunar orbit with an 80-n.mi. pericynthion, an 81-degree west seleno- 
graphic insertion longitude, and a 60-degrce equatorial sun elevation angle on the 
illuminated side of the lunar surface at insertion. 

For a 60-degree equatorial sun elevation angle at 81 degrees west longitude, the sele- 
nographic longitude of the sun must be 51 degrees west at lunar orbit insertion. From 
solar and lunar ephemeris data, it was found that the subsolar point is 51 degrees 
west at 8. 82 hours Universal Time on June 23, 1970. 

A set of patched conic non-free- return trajectories with various flight times were 
computed for arrival at the moon on June 23, 1970. Translunar injection is from 
a 100-n. mi. earth parking orbit at the earliest opportunity, i. e. , first parking orbit 
and launch azimuths are restricted to AMR range safety values, 72 to 108 degrees. 

The lunar arrival window resulting from these computations is shown in Figure 9 
In this figure, the longitude of insertion for a 60-degree sun elevation angle is plotted 
as a function of time of arrival and is called the "motion of ax". Superimposed is a 
grid' of launch azimuths and translunar flight times from the patched conic trajectories. 
Thus, the earliest lunar arrival occurs at 9.2 hours Universal Time on June 23, 1970, 
with a launch azimuth of 72 degrees, a translunar flight time of 90. 2 hours, and a 
selenographic longitude of insertion of 81. 2 degrees west. 
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Figure 9. Lunar Arrival Window 

The launch window on June 19, 1970, is shown in Figure 10. In this figure, time 
of launch is plotted as a function of launch azimuth for a 60-degree equatorial sun ele- 
vation angle at lunar orbit insertion on June 23, 1970. Thus, the launch window be- 
gins at 14. 65 hours Universal Time on June 19, 1970, with a launch azimuth of 72 
degrees. 

From Figures 9 and 10, it is apparent that a translunar trajectory may be obtained with 
an 81.2-degree west insertion longitude and a sun elevation angle of 60 degrees. This 
is within 0. 2 degree of the originally specified arrival of 81-degree west longitude. 
Translunar injection occurs on the first parking orbit with a AV of approximately 
10,280 fps, and the lunar orbit deboost AV is approximately 2,730 fps; however, the 
lunar arrival time of 9.2 hours Universal Time on June 23, 1970, is at the beginning of 
the lunar arrival window. If it is necessary to delay translunar injection to the second 
parking orbit, the grid of launch azimuths and flight times in Figure 9 would be shifted 
about 1. 5 hours later, and a more westerly insertion longitude would be necessary to 
meet the 60-degree sun elevation angle. It should also be noted that Figures 9 and 10 
are only valid for the month of June, 1970. A different synodic month would have a dif- 
ferent lunar arrival window and a different launch window. In fact, for another month 
it may not be possible to arrive near the 81-degree west arrival longitude at the cor- 
rect sun elevation angle. 
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Figure 10. Earth Launch Window 



SUMMARY 


The purpose of this paper has been to discuss some of the problems in designing a 
lunar photographic survey mission, specifically the translunar trajectory and pre- 
ferred lighting ring. The preferred lighting ring was discussed in general terms in 
order to show its characteristics over a 28-day lunar orbit. A 28-day orbit duration 
is desirable; since regardless of the polar-orbit altitude, the entire surface would 
rotate under the spacecraft in 27.32 days. 

The trajectories that were discussed were for a particular day in a particular synodic 
month and are not to be taken as general results. It is a general result, however, 
that once an arrival longitude and sun angle are specified, there is at best one launch 
window in any given synodic month. In order to open up this launch window, it would 
be necessary to vary the boundary conditions at lunar orbit insertion. 

In conclusion, the authors would like to emphasize the fact that this paper deals only 
with lunar lighting conditions and translunar trajectories. The transearth trajectories, 
lunar orbit duration, and lunar orbit aborts are complete studies in themselves and 
should be studied in great detail. 
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LUNAR TERRAIN ANALYSIS 


By 

R. L. Sanders and W. M. Greene 
Hayes International Corporation 


ABSTRACT 

N6 7 - 

The purpose of the Lunar Terrain Analysis was to analyze the disfr ibufion 
of lunar topographic and geologic features of scientific interest in order to assess 
surface mobility requirements for early lunar exploration missions. 



The study provided an analysis of lunar geology and geologic age classi- 
fication. It classified, reviewed and analyzed Aeronautical Chart and Information 
Center (ACIC-LAC) charts and United States Geological Survey (USGS) maps. 

A presentation of singular lunar sites of special interest was made, and a study 
of the Apollo Landing Belt was effected. 


The work is described in detail in Technical Report H-MOL-20, prepared 
by Hayes International Corporation, in response to Contract NAS 8-5307 for the 
George C. Marshall Space Flight Center, National Aeronautics and Space 
Administration, Huntsville, Alabama. 
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INTRODUCTION 


Work performed by the Department of Lunar and Planetary Sciences of 
Hayes International Corporation's Missile and Space Support Division at 
Huntsville, Alabama, has been designed primarily to solve problems dealing with 
meaningful exploration programs on the lunar surface. The scope of research 
has not been simply to collect and amass data but to solve specific geological 
problems. To this end, activity has been conducted in geological and geophysical 
studies, Ranger and telescopic photographic studies (with general goals of speci- 
fying astronaut activities on the lunar surface) and in compilation of data from 
existing literature. 

The geology of the moon is an important key to the history of the earth. 
The broad-based lunar exploration program now in progress is designed to pro- 
vide information concerning the origin of the moon, the earth and the solar 
system. 


Part of the program included a review and analysis of Aeronautical Chart 
and Information Center charts. These maps have been and are continually being 
prepared by the Aeronautical Chart and Information Center of the United States 
Air Force at Lowell Observatory in Flagstaff, Arizona. Morphological features 
of the lunar surface are presented at a scale of 1:1 000 000. The ACIC map 
coverage to date is shown in black on Figure 1. Areas currently being mapped, 
in various stages of completion, are shown in gray. Maps completed and scaled 
at 1:500 000 are shown lined diagonally. A typical ACIC-LAC map shown in 
Figure 2 is LAC 74 Grimaldi on the western limb of the lunar disc just below 
the equator. 
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FIGURE 




1. LUNAR CHART INDEX (ACIC LAC) 



FIGURE 2. LUNAR CHART: GRIMALDI LAC 74 



The results of the chart analysis were: 


1. Within coverage of the charts, 18 areas of geologic interest were 
analyzed on the basis of morphology. As geologic coding and interpretation are 
not shown on these maps, as they are on USGS maps, suggested vehicle traverses 
were not drawn. 

2. Features of varying nature were listed by type and by occurrence 
(duplication) . 

3. These features were listed on charts to permit rapid comparison of 
quality and quantity content of features for each area. 

Another part of the study includes a review and analysis of United States 
Geological Survey maps. 

The USGS maps have been and are continually being prepared for NASA by 
the USGS at Flagstaff, Arizona, and show geological features and geological in- 
terpretations on maps scaled at 1:1 000 000. The USGS map coverage to date is 
shown in gray on Figure 3. Areas currently being mapped, in various stages of 
completion, are shown in black. The results of the USGS map study were: 

1. Selection and listing of major geological features in chart form 

2. Suggested scientific traverses for the astronaut manned lunar roving 
vehicle were based on geologic features of interest in 18 areas of feature concen- 
tration. These areas are shown in Figure 4. 

3. It was determined by feature density analysis that circles of 80 
kilometers (50 miles) radius were optimum. Larger circles provided little or 
no significant increase in feature variety, and smaller areas usually decreased 
the variety considerably. 
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FIGURE 3. LUNAR CHART INDEX (USGS) 



FIGURE 4. AREAS OF CONCENTRATED GEOLOGICAL FEATURES 
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A typical USGS map is seen in Figure 5. The circles have radii of 80 
kilometers (50 miles) and after a thorough analysis of the map, have been 



FIGURE 5. ARISTARCHUS 


drawn to include as many geologic features as possible. Of the total 18 areas 
of concentration, these first three are quite typical of all. Area number one 
shows two traverses, each approximately 480 kilometers (300 miles) long, 
drawn along the route of a maximum quantity of geological features of interest. 


Singular or non-reoccurring lunar sites of special interest were also 
studied. Figure 6 shows 125 of these areas that have significant individual 
interests. A priority value was not established because many known factors 
remain to be solved before this analysis can be completely accomplished. The 
sites include areas which have been reported and described by various astrono- 
mers as showing, among other things, activity resembling fog, mist, smoke, 
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FIGURE 6. SINGULAR LUNAR SITES 

steam , light and coloration changes. Others were selected for unusual mor- 
phology and still others because they were extremes, such as Grimaldi which is 
the darkest area of the moon, significantly darker than any other crater. The 
sites were all selected as potential areas of scientific interest which could con- 
tribute to the geologic knowledge of the moon. 

The next five slides are telescopic photos taken from the Lunar Atlas, 
which was compiled by the United States Air Force under the supervision of 
Dr. G. P. Kuiper and were selected to present typical examples of the 125 
significant lunar features. 
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FRIGORIS, TOP NORTH CENTRAL REGION OF MOON 

Figure 7 shows the northern limb of the moon. Mare Frigoris extends 
across the center of the picture, and Plato is shown as the prominent crater in 
the foreground. 



FIGURE 7. FRIGORIS 

Philolaus (top center with two central peaks) , a crater 68 kilometers 
(46 miles) in diameter with walls 3650 meters (12 0C0 feet) high, was reported 
a scene of recent activity when a red glow within the crater was observed. This 
was reported by R. M. Baum and referenced in two books, "Our Moon" and 
"Moon Maps". 



Plato, a Post-Maria crater with a diameter of 100 kilometers (62 miles) 
and walls 1000 to 2000 meters high (3000 to 6000 feet) , has a smooth floor darker 
than the surrounding mare material, and contains several craterlets and bright 
spots. Several observers have at various times reported sighting fog or a similar 
substance responsible for obscuring portions of Plato. 

Several valleys in perpetual shadow and peaks in perpetual sunlight are 
found near both the North and South Poles of the moon. These areas seem 
valuable as areas for solar power stations and sites for making heat flow esti- 
mates from surface temperature measurements. If the lunar magnetic pole is 
located on or near the geographic North Pole, as it is on Earth, the region could 
be considered for magnetic studies. 

The Alpine Valley, 130 kilometers (80 miles) long, is shown in the lower 
right corner of the figure, and is located between Mare Frigoris and Mare 
Imbriim. It is unusually straight with a flat floor on which are a few thin cracks 
and small craterlets. 

PLATO, NORTH CENTRAL PORTION OF MOON 

Figure 8 shows an area south of Mare Frigoris. The crater Plato and 
Alpine Valley are seen at the top of the figure. Mare Imbrium occupies most of 
this photo and is about 1110 kilometers (700 miles) in diameter. It is thought to 
be impact in origin and has been accepted as an important geological time marker. 
The inner ring of mountain peaks seems to represent the original rim of the 
impact area while the outer ring of mountain chains may have been formed by a 
subsequent tectonic adjustment. It appears as though the flooding, which occurred 
after the Imbrium impact, invades the craters Plato and Archimedes (not shown 
here) and marks them as craters formed after the impact and before the flooding 
of this region. 
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FIGURE 8. PLATO 


Mt. Pico and Piton are bright, prominent peaks about 2300 meters high 
(7550 feet) thought to be formed of ejecta material at the time of the original 
Imbrium explosion. Pico, with three peaks, is south of Plato, while Piton, with 
a summit craterlet and two peaks, is west of Cassini. A cloud briefly masking 
Piton has been reported in "Geophysics as Applied to Lunar Exploration" by 
J. Green. 


HYGINUS, CENTRAL PART OF MOON 

The area in Figure 9 is located near the center of the lunar sphere near 
zero latitude and zero longitude. The Ariadaeus rille extends across the center 
of the picture for 240 kilometers (150 miles) and ends with the small, 14- 
kilometer (9-mile) diameter crater Ariadaeus. The rille, appearing to be 
broken by subsequent land slides, is probably a strike slip fault or an en echelon 
fracture. 


30 


"REPRODUCIBILIT Y OF THE O RIGIN AL PAGE IS POOR.” 



FIGURE 9. HYGINUS 

Hyginus, a crater 6 kilometers (4 miles) in diameter, has a small 
craterlet on the north rim. The Hyginus rille is roughly parallel to the 
Ariadaeus rille, leading to speculation that the two could probably be related. 
West of the crater Hyginus are large chain craters associated with the rille. 

Just below Hyginus in Sinus Medii is found the 23-kilometer (14-mile) 
diameter crater Triesnacker and its associated rille system which has an angular 
outline suggesting seismic origin. These rilles have a maximum width of two 
miles, and are believed to be the deepest on the moon. 

Julius Caesar, north of the Ariadaeus rille, is an unusual Pre-Maria 
Age crater having an irregular formation with a gap or burial of the southeast 
wall. The floor of Julius Caesar is rather interesting in that it appears to 
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slope as well as contain different age materials. The darkest portion of the 
floor on the north oide appears similar to other dark areas around Julius 
Caesar. 


SCHICKARD, EXTREME SOUTHWEST QUADRANT 

Figure 10, a photo of the southwestern limb of the moon area, includes 
the crater Schickard, a Pre-Maria ring 200 kilometers (124 miles) in diameter 
with walls about 1276 meters (4200 feet) high. The floor includes rilles, 
craterlets and areas of discoloration. Cloud-like appearances have occasionally 
been reported. 
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Wargentin, which is undiscernible in this picture but is southwest of 
Schickard on the extremity of the limb, is an important crater 90 kilometers 
( 50 miles) in diameter and flooded to the lowest point in the rim. The floor, 
containing a few rilles and craterlets, is almost 300 meters (1000 feet) above 
the outside level. The evidence of overflow has been the object of various 
theories regarding its source and the origin of lunar features in general. 

Hainzel, almost due east of Schickard, another singular site in this 
area, is a very irregular formation composed of two joined crater rings. Its 
maximum diameter is approximately 100 kilometers (62 miles) and the walls, 
nearly 3000 meters (9840 feet) high, are covered with craterlets. 

TYCHO, SOUTH CENTRAL QUADRANT 

Figure 11 shows the lunar disc with full illumination, accentuating the 
bright ray patterns associated with all the more recent craters. 



FIGURE 11. LUNAR EARTHSIDE HEMISPHERE 
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Tycho is the most spectacular feature iu the picture because of its 
large ray pattern. Tycho is 80 kilometers (50 miles) in diameter with walls 
3500 meters (11 500 feet) above the floor. Crater ray patterns probably are 
important geological features; therefore, the extreme size of Tycho's rays 
qualifies it as a significant singular feature. 

ALTI SCARP, CENTRAL SOUTHEAST QUADRANT 

Figure 12 showing the Alti Scarp, a ridge 500 kilometers (3l5 miles) 
long, is one of the Herbig series of photographs made at Lick Observatory. 





FIGURE 12. ALTI SCARP 

There is a question as to the origin of this scarp and its relationship to Mare 
Nectaris which is northeast, beyond this photo. Interesting features of the 




escarpment include the shelf below, indicating that the feature is very old; the 
height of the features, 4000 meters (13 000 feet) at its highest point; and the 
possibility of exposure of old material layers along its face. 

The southern end of the scarp intersects Piccolomini, a crater 86 
kilometers (56 miles) in diameter, south of Mare Nectaris, with a wall up to 
4560 meters ( 15 000 feet) high. The crater floor has hills, craterlets and old 
rings including a divided central hill. On the south rim of Piccolomini, there 
is a formless deposit of unknown origin and outside the western wall are rilles 
parallel to the Alti Scarp. 


LUNAR EQUATORIAL APOLLO LANDING BELT 


There has been recent interest in short traverses and terrain conditions 
at the landing sites located on or near the equ; tor. Dr. John McCauley of 
USGS in Flagstaff has analyzed this area of the moon by photometric methods 
and color-coded five different grades of surface roughness. There are 14 
singular features and some part of seven concentration areas that fall within 
this belt. An application of this preliminary study has been made to the 
geological sites located within this equatorial region. 

At present 10 landing sites are being studied within the equatorial belt. 
The sites have been selected on the basis of scientific data potential and re- 
strictions imposed by blast off and orbital rejoining criteria. 

Based on the observations and results of this first study, the Hayes' 
Department of Lunar and Planetary Sciences is currently continuing its studies 
of the lunar surface for NASA to accomplish a delineation of major geologic 
problem areas and features and, on this basis, determine major and auxiliary 
geophysical methods best applied to a given landing site. 
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A DISCUSSION OF PROPOSED LANDING SITES 
FOR AAP LUNAR MISSIONS 

By 

John R. Rogers 
| Brown Engineering Company 

ABSTRACT ~ 

N 67-24271 

Five sites on the lunar surface have been recommended for AAP 
missions: Hyginus Crater and Rill, Floor of the Crater Alphonsus, 
Hadley's Rill Region, Censorinus C, and Moltke B, The discussion in 
this paper includes cc ments on the scientific reasons for selecting 
each site and information on the nature of geological problems which 
might be investigated. A possible traverse around Site 2 (Floor of 
Alphonsus) is outlined. Recommendations for the planning of surface 
missions are made, based on new detailed information available from 
Ranger photographs. 
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INTRODUCTION 


Scientific missions at various sites on the lunar surface will be 
directed toward achieving the major geoscience objectives of the post- 
Apollo lunar exploration programs. These objectives include detailed 
studies of: 

• lunar stratigraphy 

• surface and subsurface structure of the maria, highlands, 
and various types of craters 

• lunar internal processes, such as volcanism, isostasy and 
tectonium 

• lunar external processes, such as erosion and deposition 

• thermal anomalies and gaseous emissions 

• distribution of mass in the 'subsurface 

• chemical differentiation of the moon 

Five sites which have been previously selected^’ ^ are analyzed in 
this report (see Figure 1). The major consideration for selection was that 
the sites provide a maximum in geological information, with both strati- 
graphic and structural problems considered. The discussion, in the fol- 
lowing sections, includes comments on the Eciei.Jfic reasons for selecting 
each sit** and ir.*ormafion on the nature of the geological problems that 
might be investigated. 
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1. Hyginus Crater and Rill 

2. Floor of Alphonsus 

3. Hadley's Rill 

4. Censorinus C 

5. Moltke B 


FIGURE 1. LOCATION MAP OF FIVE SITES 



PROPOSED POST -APOLLO SITES 


Site 1 - Hyginus Crater and Rill 

The coordinates of this site are 7°36' N, 6° 18' E (see Figure 2). 
The crater Hyginus is approximately 9 km wide and 880 m deep. It is 
situated near the midpoint of the Hyginus Rill, a roughly NW-SE trending 
rill approximately 175 km long. Situated along the Hyginus Rill are 
numerous chain craters. 

A site on the floor of the crater Hyginus was selected in the 
Sonnett report as a lunar landing site. The reasons for selecting this 
site were to investigate the material of the crater floor and to determine 
whether this is a maar-type crater. The belief that Hyginus is a maar 
crater is also reported in Shoemaker^ who compares Hyginus with a 
diatreme. A diatreme is defined^ as "a general term for a volcanic vent 
or pipe drilled through enclosing rocks by the explosive energy of gas- 
charged magmas". A maar crater, as used here, is a crater caused by 
a volcanic explosion; it is essentially the topographic expression of the 
same phenomenon which causes a diatreme. The most significant reason 
for visiting Hyginus is that volcanic explosion craters of this type are 
known to bring up material from deep within the interior of the planet, 
possible as much as 30-40 km below the surface. Rocks from this depth 
would be very significant in determining whether the moon has undergone 
differentiation into a crust and mantle. If the moon has undergone dif- 
ferentiation, this fact would have far-reaching implications concerning 
the cooling history of the moon and its density and elemental distributions. 
Additional reasons stated in the Sonnett report for studying Hyginus Rill 
are to detect and study any stratified rocks which may crop out along the 
Hyginus Rill and to study discontinuities in the rill and various bordering 
small craters. 




This report differs from the Sonnett report in recommending that 
the landing site be moved from the bottom of the crater Hyginus to its rim. 
The basis for this recommendation is the belief that, with almost any 
mobility system employed, operating on the rim would offer greater 
accessibility to places lying along the rill, e. g. , to the chain craters or 
discontinuities which lie along the rill. The petrographic information con- 
cerning the interior of the moon should be just as abundant in rocks lying 
on the rim as in those on the floor, since the crater origin has been attri- 
buted to a violent explosion, a phenomenon which would have ejected 
material from the bottom outward. 

One commonly used constraint is an 8 km radius of operations 
for the missions. Therefore it would be possible to visit two small 
craters located 4 km south of the rim landing site. Also, there are 
abundant deposits of the Procellarum group located around the landing 
site. Since this site is near the edge of Mare Vaporum, it is likely 
that the Procellarum deposits will not be very thick at this site and 
that a vertical investigation along the crater or rill might reveal bedded 
ejecta blanket deposits derived from several of the nearby maria. The 
superposition of these strata could provide significant stratigraphic data 
for determining the relative ages of the maria. The contact between the 
Procellarum and Fra Mauro formations should be present along the rill 
wall. If good exposures of strata occur along rills such as this, a tech- 
nique should be developed for sampling their walls since it does not 
seem feasible that any manned vehicle would be able to perform this task. 

Site 2 - Floor of the Crater Alphonsus 

The coordinates of this site are 12°45'S, 1°38'W (see Figures 3, 

4 and 5). Alphonsus has a diameter of slightly over 100 km and a maxi- 
mum relief of 3200 m. Rising from the center of its almost flat floor is 
a central peak 1080 m high. 
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FIGURE 5. PHOTOGEOLOGIC MAP OF THE FLOOR !, 
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The floor of Alphonsus was also named as a landing site in the 
Sonnett report^. The reasons stated for selecting this region included 
the following: 1) to study the dark halo craters which are probably rnaars, 

and 2) to study the rills on which some of these craters are located. The 
Sonnett report also pointed out the central peak and the Apenninian ma- 
terial on the floor of Alphonsus as subjects of geological interest. 

On a typical traverse (see Figures 4 and 5) within an 8-km radius 
around the proposed landing site the explorers would be able to achieve all 
of the objectives of the Sonnett report except the study of the central peak. 

The particular site selected in this report would provide the opportunity to 
investigate the crater wall and sample it in several places. In addition, 
the landing site is about 3 km SW of the intersection of a large NW-SE 
trending rill and a large NS trending rill. These rills shouid provide an 
opportunity to study the structural geology of the rills themselves, and the 
stratigraphic sequences exposed along them. The study of the stratigraphy 
coupled with active refraction seismic and gravity experiments should dis- 

P ** 

close the tectonics and structural development of the large crater Alphonsus. 

* 

Furthermore, the information derived from the stratigraphic and geophysical 
analysis would be of particular interest in helping to resolve the isostatic 
rebound theory of craters which have central peaks as postulated by 
Masursky^. 

Numerous small halo craters lie along the rills on the floor of 
Alphonsus. It will be especially useful to investigate these craters to 
determine whether they are maar craters or some other kind of fissure 
eruption. The implications of the maar craters have already been dis- 
cussed concerning Site 1 of this report. If they are of some other type of 
volcanic origin, they may still be very important in providing mineral 
resources and shelter for the lunar explorers. 

The floor material of Alphonsus has been reported to be 
Apenninian^. It will be important to confirm this observation. Also, 
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as in the case of Hyginus Rill, it will be important to determine the 
thickness of the floor material perhaps by finding stratigraphic contacts 
in the rills. The stratigraphic sequence will likewise be a subject of 
great importance in determining the geological history of the Alphonsus 
region. 

As at all the sites, an evaluation should be made of the impact 
process and its effect on the surficial layer. This study would include an 
investigation of some of the impact craters and sampling and mapping their 
related deposits. Since gases have been reported coming from Alphonsus, 
some emphasis should be placed on the problem of gas detec ion, particu- 
larly in planning the instrument program for exploring the chain craters 
and rills. 

The work of McCauley^ at the Branch of Astrogeology offers a new 
dimension to mission planning. McCauley has prepared a geological map 
of the floor of Alphonsus at a scale of 1:100 000 from the Ranger IX photo- 
graphs. This represents approximately an order of magnitude increase in 
scale and resolution over present maps, A map, after that of McCauley, 
and a possible traverse are included (Figure 5) to illustrate the amount of 

detail that is now available to the mission planner. 

# 

Note that the traverse shown in Figure 4 on the untouched photo- 
graph is the same as the one shown on the geologic map of Figure 5. 

Table 1 shows a detailed description of the stations along this traverse. 
Clearly, the choice of stations along the traverse are influenced by the 
surface geologic interpretation. From an untouched photograph alone, 
only a relatively few features appear to be significant. 

Site 3 - Hadley's Hill Region 

The coordinates of this site are 24‘45'N, 2*30'E (see Figures 6 
and 7). It is on the floor of Mare Imbrium, very near the base of the 
Apennine Mountains. The escarpment of the Apennine Mountains rises 
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TABLE 1 


TYPICAL TRAVERSE AT SITE 2 G';.. 30R OF ALPHONSUS) 


Station 

Distance 

(km) 

Major Geological Features 

LS 

0 

Imbrian cratered floor material 

1 

1. 5 

Contact between Imbrian cratered floor material and 
convex rim deposit 

2 

1. 5 

Indistinct crater 

3 

25 

Edge of subdued rim crater 

4 

1. 25 

Contact between convex rim deposit and Eratosthenian 
dark floor material 

5 

1. 25 

Contact between Eratosthenian dark floor material and 
concave rim deposit 

6 

0. 50 

Edge of concave rim crater 

7 

1. 25 

Contact between concave rim deposit and Eratosthenian 
dark floor material 

8 

1. 00 

Contact between Eratosthenian dark floor material and 
crater wall material of Alphonsus 

9 

0. 50 

Edge of crater Alphonsus 

10 

0. 75 

Contact between crater wall material of Alphonsus and 
Eratosthenian dark floor material 

11 

3.25 

Contact between Eratosthenian dark floor material and 
convex rim deposit 

12 

0. 25 

Edge of convex rim crater 

13 

1. 75 

*> 

Contact between convex rim deposit and Eratosthenian 
dark floor material 

14 

0. 75 

Edge of funnel crater 

) 5 

0. 50 

Edge of rill 

16 

2. 75 

Edge of funnel crater 

17 

0.25 

Contact between funnel crater and convex rim deposit 

18 

0.25 

Edge of Satellite crater 

19 

1. 75 

Summit of rim of convex crater 

20 

1. 50 

Contact between convex rim deposit and Eratosthenian 
cratered floor material 

Total 

26. 50 
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Circle Radius 8km 


(Site Number 3) 





abruptly 1800 m above the floor near this site. The exposure of rocks in 
the Apennine Mountains escarpment should pro"ide valuable iniormation 


about the rocks of the Imbrian System and perhaps the pre-Imbrian rocks. 
By investigating this escarpment, it will be possible tc test the hypotheses 
upon which the stratigx*aphy of the Imbrian System is based. Furtherm ore, 
this site should allow us to determine the nature of the basin and range 
margin and could provide valuable data on the tectonic development of 
these two major structural features. 


From the deposits lying on the floor of Mare Imbrium the minera- 
logical and chemical composition, mode of deposition, and geological 
significance can he determined. The single most striking feature of this 
site is Hadley's Rill. This is a sinuous rill which is probably not related 
to faulting or jointing as most rills appear to be. It has been suggested 8 
that this sinuous rill may have been caused fey nuee ardente (a gas charged, 
volcanic cloud) which issued from small volcano-like crater at its / 


southwest extremity. Both the rill and the crate r let ^are within 8 km of 


landing site . Located about 22 km NE of the landing s ite j ss, anothe i* 
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volcano-like feature which is much larger than the one related totherxll. 
This low conical feature, hasac rater 4 or 5 km;acr6ss and is possibly a 

- * ’ '* ° _ . V ( j ° ^ & // 

u caldera.* It is possible that, this feature cjip'hlcT be : reached by a Suhar K 
Flying Vehicle 
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arc three ages of material in the .Columbo region. These ages of material 
correspond to a region depositional blanket of material from Mare Net (ar*s 
(the oldest unit), the deposits of -broken wall polygonaTcratif^. of the 
Gutenberg type (the middle unit), and a regional depositional blanket of 
r mare rial from Mare Tranquillitatis (the youngest unit), - "L-^- 
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/; , 


: : Apparently Censbrinus C is a Gutenberg (broken, wall);type crater 

which has been masked by the ydvinger ejecta blanket from Mare TranquiJ- 
litatis. The, walls of young,' fresh, impact craters in the area of tl>e landing 

0 O w i - 

site might provide a ''window" to study the stratigraphy of the three units 

^ n L ,** * ' * - *- r ' > 

in this region which have been recognized. In addition to the ? regional 
stratigraphy the re, is photographic evide,nc,e: of local, Java ^flows iq, this . . ^ * 
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, acrpss all three of the stratigraphic units ahd are tfi^refore“ younger thani : 
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the units they cut across. ‘Another significant stratigraphic problem at 
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this site will be toicorrelate^these th r efe: re gi.onal s t r 3tig r aphic^uni t s^wi| h 
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the other stratigraphic units on the moon. Pres ejatky^, the units at, this site 
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'* j Th o e :coord^na^S o site are 1°Z0'S, 25^0'E ? (^ee Figure 9 ).^ 

Site 5. was re&omme^d^ by the .fiendix, Systerris Divisibn^ primarily 
* becaus# of rf i|l p rpximitjP tj>~ the b bund ary of a mare and highlands region: V 
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At this site, both- pre-s 



J * material of the mare?? should be withiri’ range. , In addition, it should' be pos- 
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”Sible t tb ? bhrist^uct '&/ subsurface ^apJof .the shallow-' basin ttbqr ' ih^tbia* • • 
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The LAC map of this region indicates several interesting features 
situated within an 8-km radius of the landing site. These features include 
several craters of unknown origin, ' one dome of about 3 km diameter, and 
a set of criss-crossing rills, that intersect almost in the°cent&r of Moltke/B. 

Some of jfche geological problems at this site include the nature 6f f s 
the contact betvi/etJn the mare and highlands and the shallow subsurfapfe of 


the mare, possible stratigraphic contacts, situated on the walls of the rills ~ : * 

' 11 ]' ° ’ ° #* Jj " r k> ■ 1 

or the walls of Moltke^j3> and possible volcanic^ activity or volcari hp deposits 

alid related minerjLli zationjas s 05 iate d with the dome- shaped hill. 
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The proplems at this site may be somewhat simpler than at some 
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. ' ■ of the other sites a;nd perhaps less challenging. Some of?tbe geological 

^ w - ~ „ ■’ w o *> u ^ ' ~ ° 0 lk "' ' ^ >■ 

* \\ “questions posed here may be answered on the first two, or thfe^x^pollo ^ . 
i mis sions . On the other band, the pro b 1 e m s of lunar geology may turn out 
"4 *’L td be more complex thamwe now anticipate? and' a tfite s^h L as this will be 


welcomed.' 



CONCLUSIONS 

o 

It is believed that the five sites art well suited for solvings problems 
of a regional geological natuxc. However, it is realised that there are 
many other sites which are potentially as informative ,aSs these and the se- 
lected Sites s hould not' be considered as having the ultimate in scientific 

" ' , * [;/ 
value. '■ - ‘ r ' t.. 


The results of this study indicate th vt data om, the ..lunar surface at •' 1 

/ - ’ - • " o \t . . .->• 

this time is oTJufficient detail and quality to define a major portion of the 

' ■ 

' lunar exploration problems . and O to discuss methods of solving these prob- •* ■- - 
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' ( lems. However, in view of recent advances in interpreting the Range? VII, 

-T' v. " , • ' ;• - ,r '“'"*** - •> , • 

'VIII, and IX photographs , at the Branch of Astrogeology, it should now be 

, <Js? -■ ■ • „ , r/ __ " " ft o o , 

possible to construct mission profiles "which are much mor,e detailed and ’’ 

• _ -* „ f • . , i •, 5 ■■ 

j significalnt thajk those 1 previously prepared. This investigator strongly 


recommends that a study of pnission profiles be made based oh the large- ^0 
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scale (1:100 000 to 1:1, O' 000} geological maps now available f.yom the , „ 

^Branch^of Astrogeology (Figure 5). % .. ^ " 

° p ^ t O -h t , <*, S 1 •• <•' 

41 “The mission profiles should include a Ts e t o f s ub surface models « 

• . . • \ \ * / ° • * " *\V ' 00 . . •• 

•in additio^to considerations of shrhfce traversers and their^related / . 


sampling, surface probes,- and mobility „problems\ These mpdelS would 
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" ’ predict the subsurface) geometry which is anticipated in the regionshn ' ■> 

<, 0 if | & tt “ " ' „ * ” „o ° *' « ° 

which large-Scale maps are available. -The- Vp'redictfa" subsurface, c-pn- •• s3 ° 
r figuration at, the _s^.me scales w^hich will 1^s°used . on 'actual lunar missions . 
• - could then be us c ed5|to -indicate, the responses anticipaffced from geophysical ) „ 
surveys, grilling; "and subsurface probes. o " ‘ > * « „ >. « - 
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DEVELOPMENT OF A LUNAR PHOTOMETRIC 
FUNCTION FROM EXPERIMENTAL DATA 

By 

£. F. Snyder 
Chrysler Corporation 



The lunar surface possesses some very unusual light reflection character- 
istics. Perhaps the most- striking feature is the preference for moonlight to 
he reflected toward the sun, r always predominately, in the direction from 
which it came. This effect is seen i.s the pronounced brightening of the moon 
at full coon. The need exists for an accurate photometric model' of the lunar 
surface ’$ directional reflection characteristic under all angles of solar 
illumination. Using an empirical approach and some mathematical theory, a 
functional expression is developed to describe the existing lunar directional 
reflectance characteristics, k two-dimensional equatorial photometric function 
is developed Initially, and is expanded into a three-dimensional photometric' 
function. For ten widely distributed lunar craters, and for the entire range 
of phase angles, very good agreement is demonstrated between the experimental 
results end the predicted results using the developed lunar photometric -func- 
tion. The developed three-dimensional lunar, photometric' function has been 
integrated over the face of the moon for many "phase angles. The 'integrated 
brightness at phase angle using the developed photdnietric 'function is in excel- 
lent agreement with experimental results of Rougier arid other investigators, c 




LUNAR PHOTOMETRY - GENERAL 


Lunar photometry is concerned with the investigation of the variation in 
lunar surface brightness seen with changing solar phase angles and locations 
on the lunar surface. Photometry, in general, is concerned with making luminous 
intensity measurements, luminous flux measurements, and brightness measurements 
of an object. The directional reflectance characteristics of the lunar surface 
specify the degree of surface brightness that will be seen from various direc- 
ti.ons, when the lunar surface is illuminated by the sun from any direction. 


THE MOON'S PECULIAR PHOTOMETRIC CHARACTER 


The moon has long been known to possess- some rather unusual photometric or 
reflectance characteristics. The most striking is its strong scatterback 
characteristic, which the entire lunar surface possesses to a degree far in 
excess of any natural terrestrial surface. Most very rough surfaces, other 
than the moon's, possess the characteristics that the incident, radiation is 
scattered or reflected with, nearly equal probability in all directions, in a 
manner similar to that of, the idealized Lambert or perfectly diffuse surface 0 . 

It is a characteristic of A perfectly diffuse sphere illuminated by what is 
essentially a point source, such as the sun, that its illuminated portions will 
appear equally bright regardless of the viewing angle or illumination angle. 

If the moon were a perfectly diffuse reflector, its relative intensity curve 
would be identical to the portion of the- lunar disc area illuminated at. various 
lunar phase angles. This , ' however , ‘ is very far from; being the case. |i Figure 1 
shows the visible lunar brightness integrated over the entire lunar surface „ 
relative to full tnoon brightness a S' measured by Rougier^, compared with the 
^illuminated portion of the lunar disc. Note that the illuminated portion of 
the moon is only about 1/2 as bright at a 30° phasft angle, (<X ~ ,30°),' as this 1 * 
corresponding portion atofull moon, (O' = 0). At first and last,, quarter, 

(a = 90°), the illuminated portion of the moon , is only abou£. 15% as bright as- 
the corresponding portion at full moon. ,, At 20b from -new moon, (a = 160°), the 
Average lunar brightness of the illuminated area is, lower than Che full moon 
„ brightness of the corresponding area by k factor of twenty*' 0 , “ „ > 

O , . - - - ° u - j . - 

_ ' -cj ° •- • rj ' " y ’> J ^ ^ 

^^Sven though the moon exhibits extreme variation in 0 brightness with lupair^ : 
phase . angle , the full <moon surfafea^is uniformly bright (Across it SV face;, except 
for relatively minor lo^al albedo (reflectance) variation. Limb darkening is 

However , both befgre^andafter c fvjll moon, limb 
l t he^e r mina t act closest the^sun and limb darkening is 

>'\£ . Jg, 
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completely absent at full^mooni 
brightening is observed bh-^thj^ 

observed" on the terminator larthest ‘from the sun^ Xhother * observed' effect^ 
that there is relatively 1 litflf variation in brightness,, along the lines' of. 
illumination® , longitude. " * ° ° 1 ° « ** 4 ° 
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Even at Lit brightest^ at full moon 

possessing a full moon albedoof about < 

brightness variation for various sbrface 
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the lunar reflectance is very low, 

sunlight, with relatively small 
Local iTbedos differ by 
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less than a factor of two from the average reflectance; so features of the moon exhibit 
very little brightness contrast, except as provided by shadows. At near shadow con- 
ditions, variation in surface orientation produces very marked contrast. Another 
characteristic is the yellow-grayness and/or lack of color contrast which character- 
izes the moon. The moon possesses a higher reflectance (.17) at the red end of the 
spectrum than at the blue end ( . 07) . This produces a natural yellow tinge which is 
augmented as seen from the earth, due to the deficiency in blue light transmission, 
since blue light cannot pass through the earth's atmosphere without suffering scatter- 
ing. A goal of lunar photometry has been to find a surface that simultaneously matches 
the spectral, directional, and polarimetric reflectance characteristics of the moon. 


SURFACE CLASSIFICATION 

Surfaces may be classified by their directional reflectance characteristics 
to visible radiation. Experience with terrestrial surface reflections reveals 
two common basic types of surfaces, with majority of terrestrial surfaces lying 
somewhere in between. These types are (1) a smooth mirror surface, ^typified by 
a quiet lake, plate glass, waxed table top, automobile finish, or a sheet of 
smooth aluminum foil; (2) a diffuse surface typified by a piece of writing 
paper, a piece of cotton cloth, a wood surface, and desert sands." 

For a smooth surface, Fresnel's reflection laws may be used to predict 
accurately the spectral directional reflectance and polarization. Microscop- 
ically,'’ a rough , surf ace presents a wide assortment of surface normal directions. 
The vast majority of surfaces encountered in nature are semirdiffuse, semi- 
mirror surfaces. That is, they tend to reflect light in a variety of directions, 
with the greatest reflective intensity occuring in' the specular direction. 

" f " o . c ' ■ " . - - 

' The moon, however, does not possess the semi-diffuse-, semi-specular reflect-, 
ance characteristic that is typical of the great majority of terrestrial ^sur- 
faces. Instead it possesses a characteristic that is quite rare among; terres- 
trial surfaces. The moon has a very pronounced tendency to reflect,, Tight back 
in the direction from which it came irrespective, of the direction of inc idChce . 
..This property,, is characteristic of a very porous surface,' composed ofc^vities, 

. recesses, and voids. ” * . 


• : a k general characteristic of convex shapes;, particularly' f ibets^jls that, o 
although diffuse, they still have their peak reflectance in -ttye^Torwerd" or, - 
speculai direction, especially for low grazing, .angles, “^Concave shaped surfaces 
have the opposite characteristic, which is also characteristic of the moon, 
exhibiting much greater scatterback than forward specular reflectance. 
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For a rough surface, the casual employment of Fresnel reflection lews gives 
•j& foot agreement with experiments 1 data because of the ifikAienf multiple ref Tec- jf" 
>” t ions r , involved . Although Fresnel laws remain valid at A microscopic interfacd^^ *" 
andP are , capable of simultaneously ^escribing thedirectionaiand polorimetric o j 
° reflectance characteristics, these/laws are rarely employed to'* obtain adirec-,) 
tional photometric function because 5 of the 'complexity o%<multiple reflection, ° J> 
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LOW GRAZING ANGLES - ABSENCE OF LIMB DARKENING 


At low grazing angles, every smooth plane surface displays a high specular 
reflectance. Many surfaces which diffuse radiation at normal incidence reflect 
specularly or appear smooth at low grazing angles within a degree of the surface 
(89° < i r < 90°). If, at these low grazing angles, (1) the projected features 
of the surface become smaller than one wavelength of the incident radiation, 
and (2) the features that protrude the most are essentially parallel to the 
gross surface, then high specular reflectance is observed, and the surface 
behaves as a smooth surface. At very low grazing anglejs, a rough convex sur- 
face will become smooth; a rough concave surface, however, is prevented from 
doing so by failure to satisfy condition (2). k , 

At full moon, limb darkening is absent; consequently, for low grazing angles 
on the moon, high specular reflectance just does not occur. This strongly 
suggests that the lunar surface, in addition to being highly porous, is pre- 
dominately concave in microdetail. Hapke^ has made the astute observation 
from experimental photometric tests that, for a material to backscatter light 
as strongly as the moon does, it must be extremely porous and composed of inter- 
connected cavities that open to the surface. 


SIGNIFICANT LUNAR PHOTOMETRIC CONTRIBUTIONS 


Significant photometric contributions toward understanding the nature of 
the lunar surface have been made in recent years by Van Diggelen^ and Hapke^ 4 . 
Excellent summaries of the significant lunar photometric knowledge has appeared 
in recent years by M. Minnaert^, Z> Kopal^, and C. A. Pearse?. Russian astron- 
omers have been quite active in the field of observational and interpretive 
lunar photometry, the most notable contributor^ being N. P. Barabascheff®, - 
V. A. Fedorets9, V. G, FesenkovlO, A,, V. Marltp'vll,, N. S ( . Orloval2, N. N. 
Sytinskayal3, apd V, V. Sharonov^. Russian /and American investigators in 
recent years have been largely at odds c concerhing the interpretation and 
emphasis on lunar photometric data, with tjjifej' Russian favoring a surface com- 
posed' of a light porous froth of ^volcanic /slag and the American favoring a > 
meteoric dust; surface, JSJoth groups agree M i'A high-poros ity sur face . At yet, 
no one has; succeeded i^/constructingv a\^i|sfactory physical rondel capable of » 
predicting from theory the observed,, lunary,ff:urface directional reflectance and 
polarization characteristics accurately^ MjJ , S' V - -\ " " 
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A lunar photometric function^ by 
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with “good accuracy the directional r 
surface. This _ photometric function 
„will accurately describe,, the reflect 
located on the lunar surface or in a !y ' 
was developed primarily by an empirii 
amourtt of mathematical thedry. We und| 
in functional forim the three-dimens 
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^.^-Ince characteristics of the lunar 
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orbit. b This photometric function 
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moon as it exists. The technique that was er.iployed might be called "Multi- 
dimensional analytical function development," A one -dimensional function, 
g(a) was developed for equatorial lunar reflectance using certain special 
cases 0=0, p = a, etc). This function g(a) was combined with a mi^pd 
function y(a,p). The mixed function was developed to satisfy certain sym- 
metry relations. The result was a two-dimensional equatorial photometric 
function that, fit well the equatorial observational data of M, Minnaert^ 
(1961), His data had been summarized and averaged previously by N. S, Orloval2 
(1956). The original measurements were made by Fedoretz^ (1952) and by 
Sytinskaya and Sharonov^ (1952). The developed two-dimensional expression 
was then expanded into a three-dimensional photometric function. 


TO IN COORDINATE SYSTEMS / 

/ t - - s*' ' 

ys ' ; ,i '/ 

To describe lunar reflectance, it has proved convenient to employ two dif- 
ferent tj^pes of coordinate systems. An earth-based observational coordinate 
system is defined in Figure 2a, These coordinates consist of a solar phase 
angle, a , a pair of lunar surface coordinates, a lunar surface longitude, p, 
and a lunar surface latitude, 7, where 7 = 0 is the illumination equator and 
P is measured from the sub-earth point. For operations on the lunar surface, 
we will employ a set of lunar surface coordinates. Such a cet of coordinates 
is defined in Figure 2b; these consist of (1) zenith angle, of the incident 
flux relative to the gross surface normal, (2) zenith angle, <J ) , of the exit 
ray (observer direction) and (3) azimuth angle, 0, of the exi* - ray measured in 
the plane of the gross lunar surface. This angle is measured relativeto the 
plane of incidence. The included angle between the incident and exit ray is, 
a, precisely the solar phase angle employe^ in the earth-based ,^oprdinate sys- 
tem. In the lunar-based coordinate system, an important reversible relation 
exists between the azimuth angle, 0 ahd the phase angle, c£ - ’ J<Jj. tf 

" 0 ■’ ' • " ’ <=- 0 U / /f : „ 

co3 a = cos t • cos ♦ *■ sin ty. •-sin ♦ •jtbs 0 , « ' ■ / //(l.e) 
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This relationship is obtained from the dot product of 'the incident and exit 

unit vectors. The significance and usefulness of j:his express ion/fliis in the” 
fact that once t and 0 are chosen, 9 determines ^, and vice ver-sklp 

a J “ °o" ’ « ’■ ' * " ' ' Oa \ " S 'V . 

; o cos 0 = SSL °L Z SSS £ * cos ia, • u . - • - ' ’ \ ’ (€|) 

• • " . cos (t - ♦) - "cos ]Jr • CQfJ •„ y'* - e ' 

Because the lunar photometric funefi^i depends sc strongly upon the pha^ef^ 0 * ' ■ 
angle, a, we chose to. utilizes jis_ the set of 3 independent lunar -ba s ed variables 
t, 1 *, .and a father than t," CCd\,0. , A transition between the two viewpoints * ' 

may be' effected’ by> either "equation’ (1, a) of (1 ,b ) *> ^ 

0 “ ** ^ * > O “ ,J -iVvV-. 
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« .‘‘A., second and third relation may be, derived similarly by taking dot products 
or 5 tbe appropriate unit vectors, «• - ’ *';**.« * ,-° ’ ^ 
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$09+ - cds p/. cos yj? \ - ■ „ * * Z ‘ % '' ' (2) * 

cos - cos (a - p) cos..- 7 »> - (3) 
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Taken as a triplet of equations, (1), (2), and (3) provide the means for 
making a general transition between an earth-based observational coordinate 
system of the moon (a, 3, y) and a lunar-surface-oriented system (t, 9), 

In making a transition from a terrestrial (a, 3> 7) viewpoint to a lunar 
(t, <i>, 0!) viewpoint, we will require only equations (2) and (3), 
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PHOTOMETRIC SIMPLIFICATIONS UPON THE ILLUMINATION EQUATOR 


For points being observed on the lunar equator, y = 0, equations (2) and 
(3) simplify, so that: 


cos <t> = cos 3 
cos ^ = cos (a 


P) 


4 > = ± 3 

t = ± (a - 3 ) 


(2.a) 
(3. a) 


Instead of using three independent variables (f, <t>, a), two Independent 
variables (a, 3 ) will suffice to describe the observed photometric function <=. 
for the lunar equator. Once a two-dimensional equatorial photometric function 
has been obtained, this may be .converted with relative ease into a three- 
dimensional function" in (^, 4>, a) by utiliz*- ■» symmetry relation in ^ and <t> . 


Figure 3, is a plot of the observed two-dimensional photometric function 
of the moon as presented by M. Minnaert^, withsome minor alterations which 


we made after re-examining the data of N. S. ^Orlova^. 
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Figure 3 ^serves as the ; graphical representation of the actual lunar 


equatorial ref lectance characteristics. The two-dimensional equatorial photo- 
metric function that we° developed as our best ? analytical fui>ction fit of these 
^ data is r ^ 
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F 2d ( a > ® 


bos (a - 3) 
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0 , A L " . p ay \ 

’ cos y| * \ + s:m 2 J n 


1 - sin ~ 

2 , 




. a 
cos 2 ^ 


(4) 


DEVELOPMENT OF LUNAR PHOTOMETRIC ^FUNCTION 




Several wel L, known ob.se.rvat iona 1 results may be formalized into a set of 
o general lunar //Statements. .Consistent with the following statements , we develop 
ii/irst a two-dimensional Equatorial Photometric Funct^ori, $2d (P? P) and then 
subsequently expand^ this Into athree-dimensional Photometric Function, which 
is presented in three different fox 
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F 3d <*• ♦» afo«n<*| P 3d <*» ♦»" 0 ). *^ ' F 3d A* 7): 

i - ° t> v J ^ ^ ^ , o t? 

I* At fullf moo'd, the lunar surface u is fuhifqrmil^, bright icross its face,’ 
except for [relatively mi^c,r albedo variations. s This fact will be 
expressed functionally by the* condition itwt" at zero »pha£e angle « 

//*V = A\ l*' 4 n.A«A n 4 MM n 1 - * " M a 0 — ' ■- J 


(a = 0), 



tl| the two-d^m^ns iona 1 "equatorial phot dme tr ic function. 



- <w 


J 

L 




4*. * 

Oo 


■If. 











F(a, P) and the more general three-dimensional photometric function, 
F(a, p, ;) are both unity, irrespective of the value of p and 7, 


F 2d (0, p) = 1, 

F 3d <°» P’ 7) = 1, 


n it 

-2 < P < 2 

it 

0 <7 < 9 


(5) 
(5. a) 


IX. The integrated flux from the moon peaks sharply at full moon, with 

a steep rise prior and up to full moon, and a steep descent immediate- 
ly after full moon. The trace of F(a, p) for a fixed value of p dis- 
plays a sharp apex with a slope discontinuity at a — 0, 


dFfot, p) , dF(0-, P) 

da " da 


( 6 ) 


III, The photometric function is always positive or zero and, except ;• 

’ binder rare conditions, it is confined to ah' upper limit of one, ;ij 

\ •' - - - i:l 

0 < F 2d (a, M <1 v- . (7) 

v'jl " ° " 

0 < F 3d (a, (3, 7) £ 1 „ (7. a) 






F 2d U P)*- 0, 


- 2 - P - a ~2 


J 3d 


(a, p, ?) = 0, - - % p4a-r 



( 8 ) 

* b 

(8.a), 


V. All points that have been observed the lunar surface possess a 
„ “ similar photometric function. The ibtfotypes of lunar features, oi 
„Maria (plains) and Tefrae (highlands), which present the, greatest^ 


o 


visual contrast, possess virtually „ identical reflectance Character 
istics.i Generally * the,, highlands possess agreaterdata scatter 
than the flatter region on the moon. • Large-feature ^hadowing ( prob- 
lems exist in rough terrain at low grazing angles, but otherwise, 
onephotometricfuhction serves to describe the directional ref lect-o 
ance char act efisties'of the entire lunar surface. As a consequence 0 
of , the above statement, certain symmetry conditions will” exist be- .. 
tween different observed points on the lunar surface. A phase angle", 
a and longitude, p symnetry conditidn requires that ^ 


I\?. At a non- zero lunar phase angle >(a # 0), the moon is not full, so 
that part of the lunar surface facing the earth is not directly 
illuminated by the sun. Treating the moon as, a sphere, there exists 
.. a boundry between the illuminated ;;and ’the non- illuminated portions \ 

o| the moon. „ For a >0, this illuminated boundary is defined by the c „ 
simple equation, p = a - where p is the lunar surface illumination 
longitude defined in Figure 2a. At other than fuli^tnoon^ there 
exists a darkened terminator. The photometric function assurites a ... 

. o „ value of zero, at the darkened terminator and for all points on the 
non-illuminated portion of the moon. 


o 
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F 2 d (-<*. P) = F 2d (a, -P) (9) 

F3d <" a » P» 7) = F 3d (a, -p, y) (9„a) 

where a and p may be either positive or negative numbers, 

A latitude symmetry condition requires that 

F 3d (a, P, 7) - F 3d (a, P, -7) (10) 

For non-equatorial points, combining relations 9. a and 10 produce 
fodr points on the lunar surface (at two different phase angles) 
for which the value of the photometric function under these differ- 
ent conditions will be identical, 

' >> ,, 

VI. A reciprocity lav/ involving symmetry in t and <t> states that 1 

F 3 d (^, 0, a) /cos ^ = h (ty, , a) must be symmetric in >|f and <t>. 

h (t, *, a) = h (<f, \|r, a) (11) 

ii 1 t. 

SELECTION OF AN APPROPRIATE SET OF INGREDIENTS '7 

The initial effort of this study whs directed toward developing a satisfactory 
algebraic-trigonometric expression for the lunar equatorial photometric 
function, F 2d (<2, P). It was noted that Statement. IV would satisfied by 
use of, the factor Ccs (Of - p) in the numerator of the expression for F 2d (/ 

(a, P). This factor, cos (a - P) has the characteristic that it is zero aiji\ 
the terminator or whenever p = a / j since "cos j = 0. "~J \ " ' \\V 

The choice for a first ingredient, cos , appears to have been a judicious one \j 

'because it also prepared the way for the easy satisfaction of Statement VI. ■■ „ ■> - d | 

, - v/ 11 - l . x ° C ^ I 

Statements I and IV, which define when the function is one and zero, re- < 
spectively^ are elemental to the functional development. Both tif these state- 
, mfents are satisfied simultaneously by an expression of the form: % * ” 


v t „ [SSSJSULjSLf . 

F (a, p) - cos ( aa . g* (a) ~ 


where g* (a) is a phase angl 0 e dependent function having the chafhc teHSntie-^k. 
g* (0) = 1, where n > 0, and "a M must be a number between 0 ang IT, sudh tp&w 
cos (aa - p) does not assume the value of zero in the range for jwhich - ” // 

F (a, p) is to be define/. 0 . - ■ ■ 


Equation 12 will satisfy statement VI, for n = 1 and a V 1/2 ^because on 
the lunar equator we hive: ' . \-1 •] ” . 


° 

cos (a - p) - cos.i.;.:' 

/(X \ »•»*_ + 

cos r j - p J w cos J — r* 


symmetric in t $nd ♦ 


u -if ' 

& ° . U 



—0 -w J/M 


Statements 1, IV and VI are all satisfied by two-dimensional functions of 
the form: 


F 2d P> “ y (a. P) ' s (a) ( 13 ) 

where 

y (a, P) - 

cos 2 * P 

\ . !> 

1 l /' 

The following expression provides an alternate; means df satisfying State- 
ment I and IV. ! ' 


(14) 


F (a, (3) = [cos (a -p)] ba • g* (a) 

where g* (a) is a different function of a. 

2 ' ■ 1 , 

Expression; of the form shown in equation (14) were examined but were re- 
jected primarily because of their inability to functionally satisfy (reciprocity) 
statement VI. 

, ’ . ft ' 

Equation 13 actually satisfies only the first part of Statement ,IV: that 

pertaining to the terminator boundary. This is easily remedied for /positive 
values of Of by letting: 


*2d («» P> H 




COS 


<R 


■ n « 


. 2 <'l « f P I. < * 
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u >< 
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(15) 
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r For 'negative values "of- phase angle, Of, Equation 9 ofstatement>V will be, 
used to define the function. The y (a, p) portion ofequdtion 15 , already 
^satisfies statement V. Hence, F 2 d (Of, (3) will satisfy ft as” long as g (Qf) is^ r 
symmetric about ■ 


, ' , 

1 ° o " ^ fj 

g (a) = g (-a) V 

■ • " . . \ . 


(«> 



PHASE °ANGlJ DEPENDENT 'FUNCTION, g (Of) ^ & „ 0 

' .. 

The Of, p shaping function, y . (Of^'P) - ' = /a , we ^ ave selected, 

possesses the following characteristics ^ cos - p^ 




» ft * b o , , ^ v .. 

y (a, p). * 1 , , - J < P ^ ^ (Statement I) 


P = Of - & y (a, p) "=oO 


is ' •> 


0> 


; ! ■ 


(Statement IV) 


0 p y (<*. °> - ££S ~§ > o° < a < 9o° i 
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cos 


£ - a 


y (a, f> - 
y (a. jp) = 1 


0 ° < a < 1800 

0° < a < 120° 


y (a, a) = a , o° < a < 90° 

C OS ” 
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Since F, (a, P) - y (a, P) 0 g (a), we may use the above information to de- 
fine g (a) from g (a) = ^ ^ ^ for specific p (a) curves. This results in 


cos — 


gl (a) = 


cos a 


f (O'. , 0), for p = 0 and 0° < (X < 90° 


g (a) = 


g 2 q * f „ P 3 !* 2 j» for P = J ^nd 0° < a < 180° 

h cos 2 ^ <. ' 


83 ( a ) - f (< a. — J, for p = — and 0° < a <f/20°' 


" n 


g4 (a) = cos j • f (a” a), lor p = a and 0° < a < 90° 


Using the experimental photometric curve, figuire 3, we" can select the required 
values for F (a, p) to define g (a). Note that g (a) is being defined redun- 
dantly on all but the/ range, <120° < a< 180°. This will be of assistance in 
eliminating experimental scatter, „ ^ .. 

U ^ ^ Li 

"In Figure 4 and in Table I, we presented a plot and a tabulation, respec- 
tively, of the resulting values, of gl (a) I f . g 4 (a); These different 


. . o j, ^ — * v » » oq. 

g i ( a ) are generally in good; agreement , which serves to indicate that 0 the 
functional choice of y (a,, p) = cos (a - pj>/cos ~ P) wad, a very good choice., 
Also, presented for, comparison purposes are tabular values of the analytical 
expression g (a), which has been developed to fit these experimental gi (a) °' 

• • *,84' data. It is "noted in Figure 4 that g (q) is a mona tonic ally de- 
creasing function such that g (jO) = 1 and g (it) = o. The analytical expression 
that we devel^pg^ to best fit the experimental data for g (a) has the form: 

v VQ'. pSf ■ . •, 

(«) - f . tf\2 &: , * / ' • ■ (18) 

(1 +J sin 7 ) * cos 7 »° ° 0 • „ • 

• A . '% : J ° 2 - 1„ * „ * - • " o - 


COS 7 


'ex ’ > ^ P 


>This -function wa^concbcted Is "follows : 


I. The trignome trie fvn||io^, sin 7 was selected as a suitable monatonic 
Ijinction from Which to algebraically build g (a). It is positive and 
defined for a on the interval (0, «) . \ It is monatonid and assumes the 
values 0 and 1 at the end points of the interval. ’ 0 / «, 

V^”'- 5 -> ,p " ■ » - • ? ' y ; \ .«• “ . ” % 

| o o \), „ © '< A “ • , ,, 

2 ^Different algetic combinations of sin y and 1 were tried, until a 
> er y good wa&h With the experimental data was attained with the « •„ 

^ combination: - /' p* ; , * . » «» ^ * V ", 

* .. „ 0 -• -O y o ' , v . . q , , . * 


-s' f ir 5 -^ 
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0 o’ ^ © 
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the slightly simpler results (equation 4), whenever the situucion may be fully 
described by cc > 0. 


Comparison of Figures 3 and 5 reveal a greater degree of limb brightening 
in the analytical curve, Figure 5, than is indicated in Figure 3„ " However, 
it should be realized that good extreme limb experimental data are difficult 
to obtain,, with photocell detectors because of sensor angular resolution problems. Very 
little data exist in the literature for |3 > 65°, Experimentally, using a 30X 
telescope and with the eye as a sensor, the author has determined that a con- ' 
siderable amount of extreme limb brightening does occur for 85° < p'))< 90° for 
the phase angle range 10° .< a < 50°. This. v is a region for which some authors* 
for example M. Minnaert^, have apparently made the mistake of indicating” ■ 


/ darkening on the extreme^imb near full moon. - 


<2 


>v 


6 




More recently, Pearse?, Herriman, , ! ;Washburn, and Willingharnf^ show, as we 
do, a Ljimb brightening for all value§ of , a , 0, A replot of Figure 5 versus 

the projected longitude, sin (3, rather than the actual longitude, (3, fhveals, 
in Figure 6, the limb brightening predicted by f (a, 3) as it would be obV|rv^ 
across the equatorial region of the moori. What hasp been observed experimentally 
at various ph^e angles, using the eye as a detector, appears to be consistent 
with the degree of extreme limb brightening shown in Figure 6. 

'.V ’■ J , 

INTEGRATED - TWO-DIMENSIONAL LUNAR INTENSITY r 




<S , 

o , 


3 By making the fairly good assumption that the lunar brightness is the same 
everywhere along a lunar longitude line as “it is- on the equator at that longi- 
tude, 0 the equatorial photometric function} I?2i .,{&> (3) may be integrated over 
the face of the moon to produce agood approximation of the actual integrated “ /; 
lunar intensity. This integration -of #2d (Ct, p) .will be symbolized by l2d(CO "°o 

' " r 2 " " " 




l2d «(°0 = A/ F 2d P) <*A, • 2 d I s in 


( 22 ) 
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Table II shows this in tabular ’ form, Compared wit h |Xper imenthl results^ of 
Bullrich 1948. Even for this ana ly tical appr oximabidh, the results are .rather 
good. There is “quite good “agreement for 0 < a/< 6 0°7 moderately good agreement 
-in the range 60° < OC < 1’20° and poor agreement/ in °the range 120° < a <, 180°. 
However, the fact -of thede approx ima t'e o ana ly t icdi^£e suits d if f erin^““f ro0ifi.the 8 
experimental results by more ^than a factor of^ 2.5 in the rafige where .iW; (Q!) 
is Very low would pot be at all apparent in examining a linear plot of ° i 2 d X«) 
such as is shown in Figure 7; a here the agreement look,s rattier good. ‘However 
tfris situation isinotlced in Table II. ° The analytical integrated two- “ * ° a? ° 
dimensional lun^tr Intensity Curve shown in Figure 7 my also be considered a 
plot of the ar/as under the individual curves in Figure 6; 
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of a 3d phoikstric fukctios, i 3d (tf, o, a) 

The no st difficult aspect of developing a three dimensional photometric 
function was the development of a two-dimensional expression which accurately 
described the equatorial photometric function f({v, 3). 



(23) 


The conversion of this 20 photometric function in a and 3 into a 3D 
photometric function in * r , and a, with the requisite symmetry conditions 
involving # and proved not to be very difficult. There exists a finite 
and fairly snail number of ways in which the equatorial identities: 


3 ~ ± $ and a * 0 — ± $ 

can be substituted into equation (6) # 

The combinations investigated were: 

a -» a, ± (v - o) j, , L ' 




i? t <i 'jf - 0 


it-A « a . * . V + 0 , 

2 P i * * ^ * 5 i ± 5 * * 
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The function. £(0, >Jr, a), whan divided by cos T y, must be symmetric in 0 and ty; 
so it appeared plausible to assume that 

H „ o 1 

COS (O’ - 3) ^COS f • 

Two combinations for cos (ex - 3) display the .required symmetry*, in 1/ and 0: 


cos^ - 3^ -* cos 
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The most ^.apparent substitution, for a is a itself. For Substitutions into 
( 1 - sin j) and (1 + sin no, .other substitution shows the required symmetry 

properties. However foir substitution into* cos S, cos (* ~ ^ ) is also a 
possibility. ", '.I 0 * *- * 
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Symmetry considerations reduced to four the possible ways in which the 
equatorial photometric function f(a, p) might be extended into a three-dimen- 
sional photometric function. These possible functional forms for f(^, 0, a) 


fl<#. «) = 


\!f - 


/l - si 


fl + sinSY / cos “ 


o, a) 


cos v 


+ o /, ^ . a\2 

cos — f 1 + sm j 


r qT 

!l - sin ~ 

* . Q 


£3^, *, a) * 


cos 

+ 0 

cos 


1 + sin 


fl - sin S 
1 2 

a 

cos r 


£4(^5 *» a) 


cos 


1 + si 


r- : a 

j 1 - sin 2 
I cos — 


Selection of one from tt^is set of four possible 'extended three-dimensional 
photometric f unc t ions^was accomplished by computer programming the r resulting " 
intensities in a fiviep, 7 grid and drawing the lunar- isophotes for two phase 
angles, 34° andJ?o°, for which lunar observational isophotes are recorded in 
the literatures. The final selection of a three-dimensional lunar photometric 
function wa^f 3 (^, 0, a), 'which henceforth will bie written as F3d0» r , **0!) - 


f l --« sin ? 


a) = % 

n „ 0 17 COS - T - * 

2 - 




It is important to realize that this 1 3 dimensional lunar, photometric function does „ 
resemble, yet is distinctly different from, the photometric function which would re- 
sult from taking the values of tho, equatorial 2 dimensional lunar photoirietyic function 
and assuming that the brightness remains' constaxit along the illushihation loh&tude 
lines/ Jamar photometry is -^potential tool for °%e1inotely mapping the lunar surface; t 
Under certain viewing and illumination conditions , it becomes very, important to use 
the 3D and not a 2D photometric function, because these functions may differ by aa\ 0 
much a§ a factor of 3, for slopes typical of the dunar surface. ^ ^ ” 0 
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3D PHOTOMETRIC FUNCTION , F 3d 9) 

Computet programming of the three “dimensions I lunar directional reflectance 
characteristics, F(*', 0, 9) from a lunar point of view assumed the following 
form: /" " W 

cos :• 1 r sin 2 

F3d(^, 5) = “ 


COS 


H^'(* 


sin |f v / cos f 


where 


sir* 


a 




COS G 


cos 




i + cos a 
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and 


COS Gt - COS * COS 5 * cos 9 V (1 - COS" ^){1 - COS'* 3) 


with 


£ *• £ _ jl + cos 2 * c o s * r - v(l “ cos^ 3)(1 - cos^ v) 

2 V 2 , 


(29) 


cos 

Note that this formulation, although more complicated algebraically, has 
completely eliminated the need for use of the series-, trigonometric routines. 
This expPession is in a form such that cos cos 0, and cos 9 will serve as 
the 1 three independent input variables. The results from equation (29) were 
plotted in three-dimensional form and photographed. Figure 8 reveals the 
three-dimensional directional reflectance characteristics of F 3d (%? 8 > for 

three angles of incidence £ = 15°, # = 45° and = 75°. ^ 

eo 3D PHOTOMETRIC FUNCTION, F 3d (a, 0, >) 

o 

Computer programming the three-dimensional lunar directional reflectance 
characteristics, F' (a, p, 7) from an earth-based point of view, assume the >*> 
following ^orm: 


0 ^ 


F3d (or, p, y) ?= 
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This formulation is also complicated algebraically , however it permits the 
computation to be performed without using the trigonometric and inverse trigono- 
metric routines. This expression is now in a form such that sin S, cos £ and 
cos > will serve as the three independent input variables . 2 


3D INTEGRATED LUNAR INTENSITY, 1(a) 


An integration of i' 3c* (a, £, 7) two-dimens ionally over the face of the 
moon for various phase angles, a, produces the three-dimensional integrated 
relative intensity T (a) curve for the moon. 
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This expression is most readily evaluated, numerically by converting it 
into a double sunraat ion over evenly-spaced values of cos 6 • <16 =- d fsin 61 
and cos* 7 dy ^ ^ ^ ' - • « - v 

d ^m i $ s . 
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rrT^SS r 1 ^^ 1 5 esait >l3d («>, 'from this Evaluation, are shown in &ble 
RoLi^lS i8 ^% 9, Sf* 5!S X? ! 0 ” P f! ed 4 Wit l! he **P**tei**l results of 
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: : 18 ~ , -To \r e ^parea wicn tne experimental results of 
v * 19j ^* ■ Th ? ^lytic^l results I 3d 4 a) are in excellent agreement 
'' ^ With the t experimental Jesuits W Rougier*?* 1953; whereas the approximate 

express ion I 2 g (a) ran low, tyWcally by a factor of about . three in the phdse 
range, 90° < Qf< 18p°* ciomoarison nfT anai«.f^»i rj. : 
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comparison ;o£ the analytical results for I 2 d (a) 
anfl I^d (a) leads to the exclusion that t he higher 8 results of l 3 ^ Va) a* 

co^ared to those of l 2 d ni^pt be due "to polar brightening, since on the" - 
^tor%(o: p) = Fj (a,"p;V ) #0 . „ ^ „• 7 * e on cne ; 
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EARTH BASED OBSERVATIONS. 


The developed lunar photometric function, F3d <?, a), has been expressed 
in the earth-based observational form, (a, 0, 7) and the lunar ^ surf ace 
form, F3d (V, 0, 9). The transition between these two primary viewpoints is 
readily achieved by using the simpler expression F3d (~ r , °, (x) as a common 
element in euch formulation. Comparison between experimental or observational 
results and the predicted results of the developed three-dimensional lunar 
photometric function is most easily accomplished using the earth-based observa- 
tional form, F3J (G, 0, 7). 
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Figures 10 through 13 have been prepared to show the predicted relative 
brightness of various regions of the moon at various phase angles. These data 
are presented in Figures 10 through 13 as lunar isophotes, lines of constant 
brightness plotted across the face of the moon, at the phase angles CX — 34°, 
65°, 96°, 127°, respectively. Also shown in these figures is a set of ten 
lunar craters which are more or less evenly spaced over the face of the moon 
^hat. is ' seen from the earth. Because of latitude symmetry about 7 = 0, < 
isophotes are drawn only for 7 > 0. Consistent with this situation, lunar 
features that are actually in the 4th and 3rd quadrants are shown in the 1st 
and 2nd quadrants respectively, along with those that actually are in the 1st 
and quadrants. \\ :x 


// 






' Figures 14 through 23 show the predicted brightness , Fjj (<2, 0, 7) versus 
phase angle, a for tea d if f e rentlunar craters: Vendelinus, Cleomedes, 

Posidonius, Albategnius, Archimedes, Clavis, Copernicus, Gassendi ,>ocnikhard 
and Grimaldi. Shown on each of these curves are comparative experimental re-- 
suits of ja least one of the following .observers , Blagg and Muller*” Bennett 2 ®, 
Markov and Sharonov 2 * and Fedoretz^ . These experimental values were taken 
from Van Diggeln.’s paper?. The latitude, coordinates for two of the ten lunar 
features were, noted to be recorded incorrectly in reference 3. oWe corrected 
the longitude values of Posidonius (Figure l6) and Schikhard (Figure 22) to 
correspond with values read from the lunar atlas**. Good agreement between 
experimental 0 ' observation and the analytical evaluation of Fj^ (p!,. p, 7) co- 
ncurred for all ten lunar features over 0 the full radge of phase angles. After 
obtaining very’ good agreement between earth-based experimental results and our 
proposed three-dimensional lunar photometric function, it was- decidedto shift 
viewpoints tdbchg. lurtarsurface and see what the relative brightness of the 
moon's! surface looked likieMto an individual standing on the lunar surface. 

0 Figure -24 through 30 were prej^ired to show how the brightness of t^e lunar X 
surface varies in different d^jrectipns. Isophotes, which “are lines of constant 
brightness , have? been plottedXersus the pair . of directional angles (♦, 0) de- 
fined in Figure 2b. Isophotes of F3d 9) areshownln Figures 24* through 


30 for angle of incidence f = 0 C 
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The horizontal base of each of these hemispheres corresponds to obaerva- 
t tonal data taken along the (illumination equator, and represents the directions 
in the plane containing the /surface normal -and the incident flux vector. For. ' 
'’someone standing on the lunar sur face , the outermos t hemisphere at ♦ = 90° 
corresponds to the horizon. The single point denoted as • is that point at 
which sc^tterback is occurring in the direction opposite to that of the ' inci- 
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STRANGE LUNAR REFLECTANCE EFFECTS 


Some rather strange visual effects will be seen on the lunar surface. When the sun 
is at a zenith angle, 4> = 45° , with respect to the lunar surface, the brightest direction on 
the lunar surface will be in the direction away from the sun, looking down 45° below the 
horizon. Standing and looking at this bright spot, you will see the shadow of your head 
masking the brightest zone. The direction at which the bright direct scatterback is 
occurring is denoted by <$ on Figures 24-30. The brightness will fall off from this , 
brightest spot in a pattern resembling eccentric concentric rings. For zenith angles 
greater than 45° there will also exist a darkest l’egion in the general vicinity of the 
specular direction. Superimposed on this effect will be shado\ys due to gross surface 
irregularities, and brightness variation due to surface orientation. 


When the sun is directly overhead, Figure 24, the contrast between the brightest 
nadir point and the darker horizon will only be a factor of 3. The contrast between the 
brightest, and darkest direction upon the lunar surface increases with increasing zenith 
angle. At a zenith angle of 60° , Figure 28, the brightness ratio between the darkest and 
brightest region will be about seven. 

For large zenith angles, such as shown in Figure 30 where 4> ~ 85° and the sun is, 
but slightly above the lunar horizon, the brightness contrast, even on the smoothest 
areas of the lunar surface, wall be exti’eme, differing by a factor of 50. The brightest 
spot reflectance on the lunar surface will have shifted from 5° below the horizon in 
the direction opposite the sun, up to the horizon" in this same azimuth direction. This 
effect should be most prominent on the lunar maria, where, the effect will not be ob- 
scured by the variation in slope of the gross surface features. The information contained 
in the lunar directional isotopes, Figures 24-30, have, by no means been exhaustively 
described; however, it is hoped that a sufficient description has been provided to develop 
an appreciation for the information contained and to acquaint the reader with this 3-D 
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HEMISPHERICAL REFLECTANCE OF -THE LUNAR SURFACE 
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Hemispherical reflectance is that portion* of the incidence reflected in all directions 
from a surface. The hemispherical reflectance of the lunar surface to solar radiation 
is dependent upon the angle of incidence,! Which id the solar zenith angle, ^ for a hori- 
zontal surface. ThehemispheriCal reflectance pan be evaluated by performing an inte- 
gration of the directional reflectance A fm over all hemispherical angles. 

Employing the appropriate weigitiirg function, cos <j>, to account for the projected area 
of the surface in various directions, the expression for the hemispherical refllctance is: 
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Because of variation in full moon albedo, Af m , at various locations on the lunar 
surface, it is convenient to introduce the concept of the directional hemispherical re- 
flectance function, F o ( ip ) , defined as Fo (!» =K Q («)/A m . 

Numerical evaluation of this double integral { 4>) greatly simplified by selecting 
the angles 6 at which to evaluate Fg<j ( ip, (p, 9) so that they are evenly, spaced in 9, and 
by selecting the angles (p such that the values of sin 2 0 are evenly placed. 

The double integral under these conditions reduces to an accumulative average of 
the values of the lunar photometric function Fg^ ( ip, <p, 9 ) , evaluated at a prechosen 
<p, 0 grid of angles. 
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Figure 31 shows the plot of cj^ip) resulting from the computerized numerical evaluation 
of this double integral. This directional hemispherical reflectance R 0 (ip) needed to 
> derive an accurate lunar surface temperature from a radiation heat balance statement. 




It should be noted that-tlie hemispherical reflectance R 0 for all incident angles, || 
remains less than one seventh the full moon albedos As a consequence, the lunar surface 
absorptance, Q! S = 1 - R &(ip ) , to solar radiation is typically only a few percent less" 
than unify. , 
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The directional reflectance effects from the°lunar surface have an influence dnbpth 
the deduced lunar surface temperatures apd upCri an equilibrium surface temperature of ° 
objects located on a luhar surface environmentX To lucbievb surface temperature pre- 
dictions with accuracy greater than 5^o, ± 25* F at 500°R, in a lunar environment, these 
directional reflectance effects and object shadowing effects need to be considered. A • 
detailed evaluation of these effects for spherical and cylindrical Shapes located on the 
lunar surface will be reserved fdjr another paper, which deals specifically with Itihar “ “ 
cryogenic storag_ 
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CONCLUSION AND RECOMMENDATIONS 

A three-dimensional lunar photometric function has been developed in this 
paper that does an excellent job of describing the visible reflectance charac- 
teristics of the moon's surface* For ten craters floors distributed over the 
face of the moon, very good agreement has been demonstrated between the pre- 
dicted results and experimental observations at all phase angles. The predicted 
brightness of the entire moon, 133 (a) versus phase angle is also in excellent 
agreement with experimental results* The isophote pattern lies essentially 
along the lunar longitude and agrees well with observations. The developed 
photometric function is recommended for use in performing accurate engineering 
evaluations for the temperature of objects in a lunar and near lunar orbital 
environment e 
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FIGURE D LUNARilJRFACE COORDIN ATE SYSflli 
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FIGURE 5 LUNAR EQUATORIAL PHOTOMETRIC FUNCTION FOR y 
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FIGURE 14 PREDICTED RELATIVE BRIGHTNESS OF VENDELINUS COMPARED WITH EXPERIMENTAL RESULTS 















FIGURE 15 PREDICTED RELATIVE BRIGHTNESS OF CLEQMEDES COMPARED WITH EXPERIMENTAL RESULTS 
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FIGURE IS PREDICTED RELATIVE BRIGHTNESS OF CLAVIS COMPARED WITH EXPERIMENTAL RESULTS 














PREDICTED RELATIVE: BRIGHTNESS OF COPERNICUS COMPARED WITH EXPERIMENTAL RESULTS 
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predicted Relative brightness of Gassendi compared with experimental results 
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FIGURE 29 LUNAR SURFACE DIRECTIONAL REFLECTANCE ISO?HOTES,\fr= 75° 






















TABLE I 


Phase Angie Dependent Function, g(a). Tabular value of experimentally 
based curves gj(a), g2(c<:), 83(0;) and g4 (a) compared with the developed 
analytical expression for g(a) 


a 

gj/a) 

g 2 («) 

g3<«) 

g4<a) 

g(a> 

0° 

1.000 

1.000 

1.000 

1.000 

1.000 

10° 

.809 

.813 

.809 

.812 

.803 

20° 

.666 

.665 

.665 

.665 

.662 

30° 

.569 

.574 

.572 

.570 

.556 

40° 

.460 

.468 

.46/ 

.455" 

.463 

50° 

.381 

.386 

.388 

.38 

.392 

60° 

.303 

.323 

.324 

.312 

.336 

r 70 s 

.287 

.281 

.275 

.254 

.292 

80° 

.264 

.248 

.233 

. 207 1 

.250 

90° 

- 

.221 

.196 

.166 

.219 

100 01 

• 

.194 

.170 

- 

.192 

no 5 

- 

.178 

.150 

- 

.167 

12.0° 

• 

.168 

.133 

- 

.149 . 

130° 

- 

.154 


- 

aio _ 

140° 

- 

.143 


- 

J 12 _ 

150° 


.116 

• 

- 

.... .095 

160° 

rm 

.092 

- 

- 

.074 

170° 

- 

.057 

- 

- 

.051 

180 s 

- 

.000 

- ' 

«» 

.000 
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TABLE II 


An approximate integrated lunar brightness, I 2 d (&) is obtained for various 
t phase angles, a, by assuming that the brightness of the moon assumes an equa~ 

J torial value everywhere along an illumination longitude line and then inte° 

] grating the equatorial brightness 0) over the projected area of the 

( moon. These analytical results I2d(o0 are compared with earlier results of 

; Rougier 1933 and more recent results of Bullrich 1943 


2D INTEGRATED LUNAR BRIGHTNESS , I 2 d («) 
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TABLE III 


The integrated lunar brightness, ^(a) a t various phase angles, obtained 
by integrate g 3 , 7) for 3 , 7 points on the lunar surface over the 

projected J ar surface area, and then comparing this analytical result with 
the modern experimental results of Rougier, 1933 


3D INTEGRATED LUNAR BRIGHTNESS, I 3tl (a) 


PHASE 

ANGLE 

W\ 

EXPERIMENTAL 
RESULTS 
ROUGIER 
a < 0 

ss 

ANALYTICAL 

RESULTS 

I 3d^ 

EXPERIMENTAL 
RESULTS 
ROUGIER 
a > 0 

0° 

1,0000 

1.0000 

1.0000 

5° 

.938 

.8921 

.920 

t[q£5 

.787 

.7887 

.757 , 

15° _ 

.692 

.6952 

.667 

20° 

.603 

.6109 

.586 

30° 

.466 

.4679 

.453 

40° 

.356 

' .3555 

.350 

__ _5_0_° _ 

.275 

.2682 _ 

.273 

60° 

.211 

.2010 

.211 

i__7Q? 

.161 

.1495 

.1556 

80 6 

.120 

.1102 

.a 107 

9P~ 

.0824 

.0800 

.0780 
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.0568 _ 

.0581 
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ANALYSIS AND CORRELATION OF THERMODYNAMIC 
CYCLES FOR SPACE APPLICATIONS 
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ABSTRACT 


N 67-24273 


a 


Quantitative prediction and correlation of performance 
factors, thermal efficiencies, and optimum cycle design re- 
quirements for any thermodynamic cycle or any individual pro- 
cess in the cycle was possible by the introduction of the concept 
of maximum cycle irreversibility and the utilization bf the 
principles of the Second Law of Thermodynamics. Optimization 
of the thermal efficiency to cycle specific weight ratio was ana- 
lyzed and example problems were included to illustrate the pracr 
tical applications of the concepts presented. The study concentra- # 
ted on the analysis and design requirements of thermodynamic 
^cycles intended for usage as poW'er systems in space where power 
output to specific weight ratio must be optimized. ,v 
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NOMENCLATURE 


Specific heat of gas, Btu/lbm°R 
Enthalpy, Btu/lbm 

Actual irreversibility of cycle, Btu/lbm 

Maximum irreversibility of cycle, Btu/lbm 

Specific weight of cycle equipment, 

m /Q - lbmhr/Btu 
t a 

Total weight of cycle equipment, lbm 
Total energy transferred to cycle, Btu/hr 
Total energy rejected from cycle, Btu/hr 
Entropy; Btu/lbm°R 
Heat source temperature, °R v 
Heat sink temperature, °R " 

p 

Surrounding temperature, °R 

f c 

Net power output of cycle, Btu/hr - , 


Greek Letters 


" c « 

(As. .) - Entropy change of a f , process in the cycle and the 

19(3 ^ s . surroundings, Btu/lbm°R ° o 

T) 8 L> - Thermal efficiency of cycle. , " „ g ° 


'c , .. 


Thermal effieiency<b/fra Carnot cycle 
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INTRODUCTION 


V r ith the advent of large space boosters and vehicles, large orbiting 
space platforms or manned installations .on other bodies such as the moon 
or Mars are feasible but require large output, light weight, and efficient 
internal power systems. To satisfy these requirements, thermodynamic 
cycles such as the Rankine cycle may be employed. 

One of the most important performance parameters Considered in 
the design of cycles for space applications is the ratio of the power. output 
of the cycle, W ne t, to the total equipment weight, mj.. The power output is 
related to the thermal efficiency of the cycle by 


where Q a is the total energy transferred to th;e eye 1 i. The ratio W ne t/m^. 
reduces, to r)/m s where m s = ra^/Q^, the specific equipment weight. The 
ultimate selection of the type and design of the cycle depends on the fol- 
lowing criteria: ■ 


1, . Analysis and under standing of the performance of the complete ' 

cycle and indi vid iial| c^cl e processes., and. the relationship to the- 
, ' thermal efficiency/of the cycle, ’ Q ^ ' '' 

" '■ "M 11 . 

2. Design r eq'uir emem6' for the complete^ cycle and; individual 

V- processes s'hch tha| the thermal ef£iciency^|0 speci-fi|e weight . ^ 
" ratio may be optimized, - » ,« * - ‘ £ 


$ °3, v Correlation of performance and efficiency, characteristics cor 
T a r refepopding to various types^pf cycles, ° ° a ® & 


O '*. o , - S’ O 

« Thermodynamic cycle analyses, have almost "entirely been based^oh « „ °" ,c b 

“ „ * ” coincepts of the First Law of Thermodynamic s^and, as a* result,, the above,'' 0 ., " 

“ ' ° criteria,, cahnot satisfactorily 5 >eyaluated,^ This’ study pr'^hnts a unique » ^ ° » 

i .^approach that<®ela^s, y o|t the basis of the^Second Law -of Thermod^kahiics, „ 

\ \ “ eikiciency and irreversibility q ua&titat£vely for any type of cycle (or process), % 

^ from which significant under standihg, of the above criteria may be deter^iindd, ^ 

\ « ’ 0 ®'V^The eonqept'Cf rnaximum "cycle’ irreversibility f hrta-u^ is introduced and / . <5 " 0 „ 

;• shown tq occhr when the ideal ^tended energy conversion does not'take place o*X 
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ANALYSIS 


Consider a closed thermodynamic system operating in a cycle be- 
tween a constant high temperature heat source at and a constant low 
temperature heat sink at Tj. The temperature-eiatropy diagram, T - s, 
is shown on Figure 1, where 


n 


- number of processes in the cycle. 


dQ = differential quantity of energy transferred to each 

process in system (cycle) from heat source at T, , 

h 

dQ = differential quantity of energy rejected from each 
r 

process in system to heat sink at T^, 

W^et = ne * power output of cycle, and 


surrounding temperature, (not necessarily equal to T^). 


The Inequality of Clausius shows that 


§ 


dQ 

T 


* 0 , 




* 0 , 



*0 


for a reversible or irreversible cycle of n processes operating between 
constant temperature reservoirs (Ref 1), When the total energy transferred 
in the processes, Q a and Q r , is referred to the heat source and heat sink, 
the above equation becomes 


I 


n = 1 



( 1 ) 


where the summation is equal to zero for reversible cycles and greater 
than zero for irreversible cycles. Irreversibilities associated with the 
cycle may be classified in nature as: 
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1 


where 


Then, 


. External to System : energy transfer between the system and 

surrounding across a finite temperature difference and frictional 
shear work crossing system boundaries. 

2 . Internal to System ; molecular friction, sudden expansion, dif- 

fusion, combustion, nuclear reactions, and hysterisis effects. 


A reversible cycle must be internally and externally reversible. 
The Principle of Increase of Entropy (Ref 2) yields 
As. 


isol 


As. 


isol 


As + As 

system surroundings 

As + As + As . , s 0 

system source sink 

the increase of entropy when the system (cycle) 
and surroundings (external heat source and heat 
sink) are considered as an isolated system. 


As 


system 


n 


I 

n = 1 


(As ) = the summation of entropy changes 
for each process in the cycle = 0, 


As 


source 


T / Q a \ 

) ( — — ) = entropy change of external source. 


. h' and 
n = 1 n 


*%ink 


n 


:)■ 


V / Q . 

I ^ J = entropy change of external sink. 


n = 1 n 


The irreversibility of the cycle, I act » may be given as (Ref 2): 


act 


I . = T As. .. 
act o isol 




+ As + As . . i, 

system source sink/ 


> 


( 2 ) 


- T i l + i0 

n = 1 n 
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and 


act 



^ 0. 


(3) 


Thus, the cycle irreversibility, I ac £, is the product of the 
heat sink temperature and the Inequality of Clausius (Equation 1). If 
the interaction of each individual process in the cycle with the surrounding 
is considered separately, then Equation (2) becomes 


act 



n 

V 

L 

n = 1 


As + As . , 

source sink 


n 


(4) 


where (As. ,) = As + As 

isol n n surrounding 

= net increase of entropy of each individual process. 

The maximum irreversibility, I max » of the cycle occurs when 

1. The ideal intended energy conversion does not take place, 
and/or 

2. All internal and external irreversibilities are maximized. 

Consequently, cycle irreversibilities are so large that no energy 
transferred to the system is converted to work and the heat transferred to 
the cycle 


n 


I 



n = 1 


is rejected to the heat sink across the largest temperature difference, 
- T^. From the First Law of Thermodynamics, 


W 


net 


n 


n 


I Q a + I Q r 
n = 1 n = 1 


= 0 


and Equation (3) becomes 


I 

max 



(5) 
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and 


- T l As 


( 6 ) 




max 


source 


The thermal efficiency, 77, of the closed cycle may be defined as 


n 


n 


n 


W 


y q + y q 

L a L r 


1 


7? = 


net 


n 


n = 1 n = 1 _ j n = 1 


Z S 


n = 1 


n 


n 

I 

= 1 


r. 

y q. 


( 7 ) 


n = 1 


The thermal efficiency, Tfc, of a reversible cycle (Carnot cycle, 
for example) operating between the same reservoir temperatures, and 
Tj, is given as »' 


= 1 - 


( 8 ) 


Combining Equations ( 7 ) and (8) and defining ol = V l V Q yields 


n 


O' 


i ° r 

)• < osasl -°> 

I Q a 


n = 1 


Adding 77 and subtracting (1 - T^/T^) gives 


n 


l ** T 

l \ 


n = 1 



and 


1 






Combining Equations (9), (6), (3), and (4) gives 


(9) 


n 



max 


Thus, Equations (10) and (11) relate the irreversibility of each 
individual process in the cycle to the maximum cycle irreversibility, and to 
the ratio of the cycle thermal efficiency to that of a reversible engine 
operating between the same reservoir temperatures, T^andTj. The 
magnitude and significance of the losses in efficiency corresponding to 
each individual process in the cycle may be evaluated or the over-all loss 
in efficiency may be determined directly from combining Equations (4) and 
(10) to give 
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or 


ot 


(12) 



sink 


-) 

I'r 1 n' 


source 


Figure 2 shows a plot of Equation (10). Calculations of a and 
^act^max ^ or thermodynamic cycles must lie on the same straight line. 


DISCUSSION 

It was assumed that the temperature of the heat source and heat 
sink processes, T h and Tj, were constant but, in general, they may vary 
as shown on Figure .< (dotted processes). Under these conditions, Equa- 
tions (10), (11), and (12) refer to constant temperature source and sink 
processes cori esponding to the maximum and minimum temperatures of 
the cycle, and T^. The number of processes in the cycle must include 
the variable temperature source and sink processes in 


n 

n = 1 

Although the irreversibilities between the source processes and 
sink processes do not actually occur, they are inherent in the system and 
represent a loss in efficiency which cannot be prevented unless the reservoir 
processes occur at constant temperature. For example, suppose the 
thermodynamic cycle shown on Figure 1 operates completely reversibly 
between the variable temperature heat source and sink processes. The 
thermal efficiency of the cycle, 17 , cannot equal that of the Carnot cycle; 
i. e. , a < 1, due to the inherent irreversibilities between the heat source 
processes with constant and variable temperatures and between the heat 
sinks processes with constant and variable temperatures. This concept is 
demonstrated in an example given later. 

In the design of any thermodynamic cycle, the total irreversibility 
of the cycle must be reduced to increase the thermal efficiency ratio, a, 
as seen from Equations (10) and (11); i. e. , the irreversibilities of one or 
more of the processes must be reduced. If the irreversibility in a process 
is generated by heat transfer across finite temperature differences, larger 
surface areas can be used to reduce the temperature differences required 
to transfer the same amount of energy. However, internal irreversibilities 
may increase as a result of increased surface friction effects. An optimum 
point usually occurs where the total irreversibility (external + internal) is 
minimized. 

However, reduction in cycle irreversibility to increase the thermal 
efficiency, T), may be associated with an increase in equipment specific 
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weight, m g (ratio of equipment weight to total energy transferred to 
cycle), which is undesirable for space applications. Process redesign 
to reduce irreversibilities may not appreciably affect the weight of the 
equipment such as in the redesign of flow passageways to reduce separa- 
tion losses, turbulence, and frictional effects or in the redesign of 
combustion chambers to promote better combustion efficiency. The 
ultimate design criterion is the ratio of thermal efficiency to specific 
weight, T}/ m g , which was shown earlier to be equal to the ratio of the 
power output of the cycle to total cycle equipment weight, W net /m|.. 

Although it is not the purpose of this study to present design data 
for the method of design to optimize Tj/ m s for the cycle may be 

generalized. Figure 3 shows the qualitative effect of reducing the 
irreversibility [Tj ( ^ s isol)n^ma x ^ a heat exchange process on the 
specific weight of the heat exchanger, (m g ) n . The corresponding increase 
in the thermal efficiency, Tj, evaluated from Equation (10), is included. 
The maximum in 77 / ( m g ) n , if one exists, indicates the optimum heat 
exchanger design. A similar analysis for all processes in the cycle can 
be conducted to yield a maximum ri/m s for the cycle. 

From Figure 3, it can be seen that 




n = 1 


max 


0 


to give a maximum for the entire cycle. Substituting O' = 7 )/ TJb in Equation 
(11) and dividing by m s yields 


where 



Assuming T^ and T^ constant gives = constant, and letting x = lact^max* 
the above equation becomes 

g (x) 
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and 


0 


(rj/m ) 


[. g(x) - (1 - x) g'(x)3 s 
g(x) 2 


which gives a maximum when x satisfies 

g(x) = -(1 -x)g'(x) 

or = _ _1_ 

g(x) 1 - X* 


T 1 I 


Since x = 


n = 1 


the above equation yields 


T I I < is i S0 l>n 



T 1 I l^isol'n 


0 - r^) 

' tviav ' 


and the maximum 7j/m s results when the individual processes are re- 
designed such that the function 


T 1 l 

<-—4 ) 


satisfies Equation (13), if possible. Of course, the function must neces- 
sarily be obtained analytically or experimentally for each process in the 
cycle. 

The same analysis may be conducted for various types of cycles 
and the ultimate selection of the type of cycle for space applications 
necessarily dictates that T )/m # be a maximum. 
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EXAMPLE PROBLEMS 


To illustrate the practical usage of the above concepts and to get a 
better understanding of their significance towards cycle analysis and 
design, the following example problems were considered: 

A. Rankine cycle steam power system installed in an orbiting 
space platform. Figure 4a shows a schematic diagram of the power sys- 
tem and Figure 4b depicts the corresponding temperature -entropy, T - s, 
diagram. Gases, heated within the nuclear reactor (external heat source), 
are cycled through the boiler and transfer energy to the boiler. There- 
after, the gases pass through the superheater, reheater, and economizer 
before returning to the reactor. To improve the thermal efficiency, a 
superheater, reheater, economizer, and two extraction feedwater heaters 
were included in the basic Rankine cycle. The energy rejected in the cycle 
is accomplished by the use of a radiative type condenser where the earth 
is considered as the external heat sink. Refer to Table I for system data 
and values of enthalpy and entropy at each point in the cycle. Entropy changes 
for the constant pressure gas process were calculated from 


s - s. 
e 1 


= — £ c 
“b P 



Neglecting kinetic energies, the temperature changes of the gases were 
determined from 


(T - T.) 
e i 


Q m 

rn b 
c m 


P g 


where Q rn is the energy transferred from the gases to each respective 
process in the cycle. 


The losses in thermal efficiency of the cycle resulting from 
irreversibilities generated in each process of the cycle were evaluated 
from Equation (10) 


ot 


12 

l <Wn 

n = 1 



12 

Z <Aa isol*n 
n = 1 


I 

max 


and are given in Table IL It is noted that a = fj/lfc s 46.25%; i. e. , the 
losses in efficiency amounted to 53.75% for the 12 processes in the cycle. 
The most significant losses result in the following processes: 
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Boiler process--13. 14% loss primarily due to external 
heat transfer across a finite temperature difference. 

2. Condenser process--4, 61% loss resulting from irreversible 

radiation energy exchange with the heat sink (earth). 
Generally, this loss will be larger for applications in areas 
of space void of mass where the effective heat sink tempera- 
ture approaches absolute zero; i. e. , = T q = 0 +o R. 

3. Superheater process -8. 06% loss due to external heat 
transfer irreversibilities. 

4. Gas -nuclear reactor process - -external heat transfer 
irreversibilities yield 8.54% losses, 

5. Turbine c process--8. 23% losses resulting from internal 
frictional irreversibilities produced because of inefficient 
turbine design (turbine blade design, for example). 

6. Low pressure feedwater extraction proc«ss--4. 98% loss 
resulting from external heat transfer irreversibilities. 

Other losses were individually small. Cycle redesign criteria 
may be determined by knowing the types of irreversibilities produced and 
how to eliminate them. Optimization of Tj/m s may be ascertained from 
methods given previously. It is important to note that redesign of 
individual processes may affect the performance of other processes in the 
cycle. Thus, optimization of each process in the cycle must be interrelated 
to each other to insure optimum performance of the complete cycle. 


To determine the maximum value of ot for the cycle, it was 
assumed that the Rankine cycle operates completely reversibly between 
the variable temperature heat source process, Tj, -• T ec , and constant 
temperature heat sink process, Tj. From Table II, 


Ot = 100 

max 


I (minimum) 

ctCT 


max 



= (100 - 8 . 54 )% = 91 . 46 % 



which yields a loss of 8.54% due to external irreversibilities in the gas- 
nuclear reactor process. In addition, suppose that the nuclear reactor 
was eliminated and that gases at temperature were continually supplied 
to the heat source process, -* T ec , from some gas supply (combustion 
process, for example). Calculations showed that the loss in efficiency 
was 8. 54%, the same as with a constant temperature heat source at T^. 
Thus, inherent irreversibilities result from the use of variable tempera- 
ture heat sources (and sinks). The equivalent heat source temperature is 
considered constant at T^ and the gas heating process, R Q , does not 

actually occur but must be effectively considered as the reverse of 
process T^ T ec and included as an additional process in the cycle 
(identical to conditions in example above). 

Although the Rankine cycle shown on Figure 4a is typical of 
modern power systems on earth, some of the components tending to reduce 
rj/m s (feedwater heaters, economizer) resulting from excessive equipment 
weight must be eliminated for use of the cycle in space applications. 

B. Thermodynamic cycles involving combustion processes have 
commonly been approximated by the air -standard cycle (Ref 2). An 
example of this is the Otto cycle shown on Figure 5 where the combustion 
process is replaced by energy transferred to the system from a heat 
source at and the cycle is considered closed (no change in composition). 
The temperature and volume limits correspond to maximum power condi- 
tions of modern internal combustion engines (410 hp). Assuming the cycle 
to be internally reversible, an analysis similar to above yielded 
Of = 7]/ Tic = 0.673. To account for deviations from air-standard approxi- 
mations and irreversibilities such as ’ ^at transferred to the surrouridings, 
Q r , in processes (1 -* 2), (2 -♦ 3), and (3 -* 4), internal frictional irreversi 
bilities, and external friction in bearings and reciprocating components 
(Wfriction = Qy4> en ergv dissipated to the surroundings), the actual W ne ^ 
output of the cycle is given by. 

W net (actual) = e b W net (ideal) 


where e^ — 0. 65 = Mechanical efficiency of a well-constructed engine. 
Thus, the actual value of Of is 


Of = Of e, = 0.438. 
b 


To include these effects in Equation (10), 

4 

T y (As. ) + 

1 L isoln 


«»a a <* - %> 


Of = -2- = 




max 


— ) 
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where processes n = 1 to n = 4 may consider the actual internally and 
externally irreversible processes occurring in the cycle (a < 0. 673). In 
this case, e^ includes only external frictional irreversibilities. 

From Figure 3, it is seen that a = 0.438 is less than 
a = 0.4625 for the Rankine cycle. However, the thermal efficiency to 
specific weight ratio, 77/m g = a T] c /m s , is higher for the Otto cycle. 

Applications of the above concepts may be extended to pro- 
cesses since cycles are merely specialized processes. Detailed analyses 
for processes and other cycles such as gas -particle processes in rocket 
nozzles and combustion driven Brayton cycles have been conducted, 
respectively. 

It was shown that the relation a = f (1 - I act /I max ) can ke 
used for any process (a redefined for processes), including thermo- 
dynamic refrigeration cycles. For the latter, 

1 1 

n 

I <* 

1 .( IfJ 

\T. - T J As 

N h l' source 

where (3= Q a /W net = coefficient of performance and" I max = - (Tj/T^) I max » 

.Concepts of availability of energy were interrelated with the 
analysis given herein. The applications and results will be presented 
in later studies. 

Ackeret (Ref 3) presented a historical background of entropy 
and other examples of the role of entropy in the aerospace sciences. 

CONCLUSIONS 

Significant conclusions of the analysis were: 

1. The concept of maximum cycle irreversibility interrelated 
with the actual irreversibility of the cycle, reflected in 
Equations (10) and (11), successfully determine cycle 
performance and thermal efficiency design criteria 
quantitatively for any thermodynamic cycle. 

2. The performance and efficiency /specific weight of any thermo- 
dynamic cycle may be quantitatively determined, compared, 
and correlated. Graphical representation of Equation (10) is 
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shown on Figure 2 which relates efficiency and irreversibility 
of all cycles. Results (Table IQ for the Rankine cycle shown on 
Figure 4a, and those corresponding to an air-standard Otto 
cycle (Figure 5) are included. 





Heat Sour< 



Figure 1. 
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As . . (max) — 
sink 


T emp er at ur e - Entropy^ Diag ram 


Thermal Efficiency Ratio 



Figure 2. Efficiency-Irreversibility Correlation for 
Thermodynamic Cycles 
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Figure 3. Irreversibility, Efficiency, and Specific Weight 
Characteristics of a Heat Exchanger Utilized in a 
Thermodynamic Cycle 
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Figure 4b. Temperature Entropy Diagram of 
Rankine Cycle Power System 
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T emp e r atur e 



Entropy - Btu/lbm°R 


Figure 5. Temperature Entropy Diagram for 
Air-Standard Otto Cycle 




TABLE I 


RANKINE CYCLE POWER SYSTEM PROPERTIES 


DATA 



e. 

l 


0 ^ 

m g /m b 

m 

x 

m 

y 


Boiler Pressure = 2400 psia 
Condenser Pressure = 2 in. hg.a'us. 

Reheater Pressure = 1200 psia 

Feedwater Extraction Pressure = 1200 psia 

Feedwater Extraction Pressure = 600 psia 

Temperature of Nuclear Reactor Heat Source = 3310°R 

Temperature of Heat Sink (earth) = 520°R 

Exit Superheater Temperature = 1460°R 

Exit Reheater Temperature = 1460°R 

Internal Efficiencies of Turbines a, b, and c = 0. 8 

Specific Heat of Gases = 0.27 Btu/lbm°F 

Ibm of Gas per lbm of Steam Generated in Boiler =2.0 

Mass of Steam Extracted at = 0. 107 lbm/lbm (boiler) 

Mass of Steam Extracted at P^. = 0.263 lbm/lbm (boiler) 


Condensate pump processes assumed reversible. All results are 
referred to 1. 0 Ibm^ of steam generated in boiler. Interconnecting 
line losses were neglected. 



TABLE I 


(CONTINUED) 


State Point 

h-Enthalpy 

Btu/lb 

s -Entropy 
BtuAbm°R 

1 

69.1 

0. 1316 

-> 

lm 

76.27 

0. 1316 

y' (y - y') 

472.0 

0. 666 

y f 

471.6 

0.672 

x 1 ' ;(x ss x 1 ) 

570.0 

0.765 

x f 

571.7 

0.7711 

3 

718.4 

0. 9023 

4 

1101.1 

1.2434 

5 

1461.7 

1.5335 

6 

1371.0 

1.5335 

6' 

1389.2 

1.5485 

7 

1499.2 

1.6293 

8 

909.1 

1,6293 

8' 

1027.2 

1.8556 

9 

1401.0 

1.6293 

9' 

1420.8 

1.646 
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TABLE II 

CALCULATED RESULTS FOR RANKINE CYCLE POWER SYSTEM 




* 
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(100.0 ■ ».75)% = 46.ZS* , - (l 25 * = 38. 8% 
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MATHEMATICAL ADVANCES IN THE THEORY OF ONE DIMENSIONAL FLOW („. 


By 

' Ervin Y. Rodin 
/ Wyle Laboratories , 


SUMMARY 


N67-24274 

This paper contains the derivation of a Burgers 1 equation, as a definite approximation 
of high order, from the one dimensional Riemann pair of equations (with a viscous term) 
which is a description of the motion of waves of finite amplitude. The solution of the 
Cauchy problem for Burgers' equation is presented, together with a discussion of solutions 
by means of infinite products. 



INTRODUCTION 


The ambitious plans of the free world in space travel are predicated to a large extent on 
the ability of scientists and engineers to cope with a demand for knowledge that is rising 
in an exponential manner. The recent great strides in this field are ascribable to two 
important factors: one, the development of large and accurate thrusters and two, the 
appropriate miniaturization of much of the component and support equipment. These, 
in turn, were made possible by basic scientific studies performed earlier in this century. 

There is, however, a very practical limitation both to miniaturization and to largeness 
of rockets. Thus, in order to continue the rapid advance in this field, it is necessary 
to refine many of the existing techniques or theories. For instance, much of the lifting 
capacity of a thruster today is used to put into orbit a considerable mass of ignorance - 
many call it a safety factor - which is necessary; because all calculations used in design 
and operation are approximations of varying exactness. Thus, approximations of higher 
order will reduce this dead mass and refine the processes and equipment used. 

Basically, this is the subject of the present paper. In the hierarchy of approximations 
used for describing various phenomena, such as acoustics of waves of finite amplitude, 
properties of highly conducting gases, etc., exact solutions of the governing equations 
are available only in highly specialized and idealized cases. One uses these exact 
solutions as stepping stones, from which solutions of more real problems can be obtained 
approximately. Thus, the higher in this hierarchy one is able to obtain exact solutions, 
the better the further approximations are . 

One of the more important of these stepping stones is the model equation, proposed in 
1940 by J.M. Burgers, whose general solution was obtained in 1950 and 1951 by E. Hopf 
and J. Cole. This nonlinear partial differential equation, while being a descriptor of 
phenomena only in one space dimension, still includes an important nonlinear feature and 
also takes into consideration effects of dissipation and diffusion . The functional form of the 
equation was shown to be extremely useful in many fields; see the Bibliography. 

The only important boundary/initial value problem solved so far in connection with this 
equation is the one describing the decay of a pattern existing at some initial time. It is 
of great importance, however, to know the behavior of a system with excitation condi- 
tions described at one of the physical boundaries; the solution of this "Cauchy Problem" 
is one of the subjects of this paper. 

The general solution of Burgers' Equation has, usually, the form of the quotient of two 
infinite series of rather complicated terms. Thus, it is not easily adaptable to a qualitative 
analysis, and is very cumbersome quantitatively also. 

In order to cope with this problem, tne second part of this paper obtains same of these 
solutions in the form of a single infinite series. The approach is through infinite pro- 
ducts; a mathematical object which has been used in practice but very little; but one 
which could have an importance in the solution of nonlinear problems that is comparable 
to the utility of infinite series in the theory of linear differential equations. 
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The first section of this paper contains a derivation of Burgers' Equation. It consists 
essentially of the justification for the reduction of a pair of one dimensioned coupled 
equations to the one used and comes from the theory of nonlinear acoustics. This is 
followed by a description of the properties of Burgers' equation and of its solution, 
which is a prelude necessary for the next section, containing the solution of the Cauchy 
Problem . Finally, the last section of the paper contains solutions by means of an 
infinite product approach . 


THE EQUATIONS FROM NONLINEAR ACOUSTICS 


One interpretation of the equations 

2 

(1) v, + vv + r aa = ev 

t x y - 1 x xx 


( 2 ) 


y- ] 

a. + a v + -L-«— av = 0 
t x 2 x 


is that they describe the propagation of plane sound waves of finite amplitude in a 
dissipative fluid (see Lighthill) . Here v is the particle velocity, a a quantity 
proportional to the local speed of sound, y the ratio of specific heats of the medium 
and € is a dissipation constant. The subscripts, as usual, denote partial differentia- 
tions with respect to the letter involved. The system (1) - (2) is a classical Riemannion 
one when e = 0; it is then solvable by means of its characteristics. In its fullness, the 
pair (1) - (2) has not been solved. We shall now present, however, an argument which 
reduces (1) - (2) to the simple equation 


v A + vv = 6v 
t x xx 


where 6 is also a constant . 
The substitution 




a = a. 


d) 


reduces (1) - (2) to the system 

. . .2 y-2 

v + vv + k o' 


(3) 

(4) 

where 


\ + vv^ + k p' p, 

p t + K = 0 


= ev 
x xx 


k 


2 


2 1 
°0 p 0 


r 


is a non-negative quantity. 

Let us introduce a stream function f , defined by 

e> p - ♦„ . v = . 
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Then equation ( 4 ) is satisfied identically by ^ ; there remains to satisfy ( 3 ). 

To this end, we shall assume that 41 satisfies the equation 

(6) d> = c ' <t> 

T t 1 T xx 

for some 0. This, of course, also means that p satisfies the same equation: 


(7) p , - C 1 P xx ■ U - C 1 »xxx = 0 

In fact, we can write for v, from ( 5 ), (6), and ( 7 ): 


( 8 ) 






__ J_t T XX _ ‘X 

v 4 > ~ C] K - ' C 1 P 


Let us use this value of v in ( 3 ); the substitution and subsequent simplifications follow: 


PP~* - P„ P* 


1 2 

P 


xt *x r t 2 P x PP xx p x , 2 y-2 

+ c, X + k p' p = 


'Ip 2 
r P 


-«c, 


P P xxx 2 pp x P xx + 2p x 


p p xt - pp x P ‘ 


t C 1 ( PP x P xx p x | c^ P p x * ^ 


p p - 3 pp p + 

XXX X XX 


2p ‘) 


- 3 ‘) p *x) + ( c i - 2< ) p x - 


p r+1 p =0 

r r x 


( 9 ) 


( e l - p2 P xxx - ( 2 e ! - 3 ‘) P xx + ( e I - 2 ‘) p x * ^ 


y + 1 

P' P. 


The cases 


« , n = 0, 1, 2, 2, A 
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lead to special situations/ corresponding to various physical restraining assumptions; 
we shall not consider them here. 

By assumption (6), we have that p satisfies 


*xt C 1 *xxx ' c l^ ® 


Let us recall that, for increasing values of c^ , the solutions of the heat equation 
approach bounded values in a manner which' is usually of the exponentially damped 
type. Furthermore/ derivatives p of the solution approach 0 for irtcrenriny c. . 

X I 

With this in mind/ let us write (9) in the form 
00) p 

and consider the effect of letting c^ increase. Clearly/ since by the preceding 

argument the numerator of the right hand side of (10) is bounded, it becomes very 
small for c^ large; in fact, it is of the order o (cj ) , where a<2. Therefore, 

choosing 


xxx 


2cj - 3e 
C 1 " € 


c ] ~ 2e 3 !<* 

pp p + . —— p = ' i —• 

rr X r XX Cj - € r X {Cj - i'. C, 


c, " 6>[ 1 ; e; k 2 M^ M J , 

M = max p in some region R 
= max p^ in the same region 

2 3 

pp - 2 pp p + p = 0 

K r xxx "jtxxx K x 

This expression can be integrated; using the letter p for the upcoming functions 
(of t) of integration, we obtain 

pp i w ' 

or 

03) - P,W • 

r 'x 


(id 

where 


we reduce (10) to 
(12) 


i' 1 
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Letting now 


P = 


u 

e 


(13) becomes 


= p, (t) e U u 
x xx H x 

which can again be integrated: 

\\ = -p,(0 e ' u + p 2 W 


and thus 


04) 


= vr 


ij P 2 W “ P,Tt) e‘ U 


One more integration , and a transformation from u back to p gives us 


05) 


P = 


p^W 


. 2 

sin 


yj - P 2 W ‘ 


(VT x + p 3 (»)) 


where we must have Pj (t) (. 0 and pjftK ® 


This function then is the justification for the assumptions made on p a priori, for 

an appropriate choice of the arbitrary functions p.(t) . In particular, let us note 

that p, as given by (15), is a non -negative function; a fact which will further corroborate 

the appropriateness of our choice of p as a solution of the heat equation in later 

developments. 

Our conclusion, however, is the following: if p satisfies 


p = 8p 
xt XXX 


with 6 given by (11), and 


s -5 


r x 
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then equation (4) is satisfied identically. Furthermore, v then satisfies the 
equation 


v. + 2 vv - 8v 
t x xx 


(as vviil be shown in the next section); while the magnitude of the error in using an 
approximation for equation (3) is 



2 ,2 

y-1 t 


E = 

v k exp 

2 y-1 

— /v dx 

+ (6-e) v^ + G(t) 


where G may be freely chosen . Furthermore, we have shown that E is negligibl 
for a 8 chosen large enough, in some region of interest R. 

Our task, then, is to solve the equation 


v. + 2vv = 8v ; 
t x xx 

and it will turn out that the most general solution is 



where 0 is any (positive) solution of the heat equation. 
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BURGERS EQUATION; INITIAL VALUE PROBLEM 


The crucial step in the developments of the previous section was the assumption that p 
satisfies a third order heat equation. However, in order to make that assumption, it was 
necessary to know that its consequences will be desirable ones. In this section we are 
going to show how the heat equation enters the picture independently of any other 
considerations. 

We shall start with Burgers' equation in a slightly more general form: 

(16) v. + Bvv = 8v 

t r x xx 


where (3 and 8 are constants. 

In order to give a measure of physical judgement to this equation, let us multiply it 
by v and integrate on the interval (xj, Xj), to obtain 

x 2 

(17) \ fit ^ dx + ■§■ { y3 ( * 2 • “ y3 = 

*1 


= 8 jv(x 2 ,t) £ £ («,.») ' / (£) * j 

X 1 

Going in order from term to term, we can interpret this equation as one stating that the 
total rate of change of kinetic energy in the system, plus the net flux of this energy out 
across the boundaries, exactly balances the rate ct which work is done at the boundaries, 
less the total dissipation present. 

We now proceed to find the solution of (16); following, substantially. Cole. First, we 
return to the stream function of (5), letting v * f and integrating (16) . If we take the 

constant of integration as 0, we arrive at 

< ,e > ♦. + 4 ♦* - ■ 

Let us now assume that f = f(h); i.e., that g depends on the single function h (of the 
variables x and t) . Then (18) can be written in the form 


n 


V 1 » ■ " -» 




/s - 


1 
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8h + h‘ 

xxf X 


If v/e now assume that h satisfies the heat equation. 


h = 8h 

t xx 


then f has to be a solution of 


f" ; 


which means that 


= f = - -a- In (h + Cj) + c 2 


with Cj and C 2 arbitrary. Thus, 


v = <t> = - 

T x 


28 h x 

7 7 


is the solution of (16). The constant C 2 disappeared and Cj was incorporated in the 
solution of the heat equation. 

Let us remark at this point that no rigorous mathematical proofs are included in this 
paper; however, we shall make a precise statement at each point. Thus, in this case 
we have 

Theorem 1: The following two statements are equivalent: 

i.) v is a solution of (16) in an open rectangle R of the x,t plane 

A 

it.) v, v , andv are continuous in R 
x xx 


i.) For t>0, there exists a function h of the form 


h - C(t) exp (-*/* fr,t) ax ] , 


such thot 

ii.) h>0andh,h , h , h are continuous in R 
x' xx' xxx 

ill.) h « 6h inR. 
t xx 
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It is Interesting to note how ine constant p , representing finite amplitude effects, is 
balanced by the viscosity constant 6: 


( 21 ) 



What we have achieved so far is to reduce the general solution of a nonlinear equation 
to that of a well known linear one. A question that immediately arises is how the initial 
values are related. From (21) we have 


(22) C(t) exp [- jjg/v(x,t) dx J = h(x,t) 

so that, if 


v(x,0) = f(x) 


we obtain the relation 


(23) h(x,0) = C Q exp [- ^ /f(x) dx] 

as the initial condition to be used in the solution of the heat equation. However, we 
shall not pursue this any further, for this question has been studied quite extensively 
both by Cole and by Hopf . Instead, we shall investigate the question of how to relate the 
condition given on a line x = constant for v to h. 



SOLUTION OF CAUCHY PROBLEM 

The problem we shall solve in this section is the following: find a twice continuously 
differentiable solution v of the equation 

v + pvv = 6v , 6, 8,> 0 . 

t r x xx r 


satisfying the Cauchy conditions 


v(x Q ,t) = a(t) 


v x (x 0 ,t) = b(t) 


Clearly, no generality will be lost if we take p = 1 and Xq 
W e know that v is given by 


h (x,0 

= - 26 ‘iM 


where h satisfies 


h = 6h 
t xx 


Let us assume a solution h of the heat equation that is of the form 


(n) , . 2n (n) M 2n+l 

(29 h<*,t) = 23 ir + ^ a ~^‘ ‘fern • 

iTo 8 n “• n = 0 8" ' 2n + " • 

where w and f are arbitrary functions of t and where 


o«0». ££ , Q (0) w 


for any function Q . It can be easily seen that 

8(0,0 = «(th 

•v (0,0 = f(t) , 


154 


Thus 


'-J 0 '*) ■ ?“'<*> 



In order to satisfy the conditions (24), we therefore must have 


and 


Thus 


so that 
(28) 


-26 


1»(t) _ 

wltj 


«(t) 


-26 


«'(t) 

6*»(t) 



o'ft) _ 


u (t) = exp 


^ (° 2 (0 - 26 b(t)) 

[ ^ / (° 2 (0 - 26 b(t)) dt] 


(29) *(0 = - ^ u(t) 

Therefore we see that the conditions (24) "translate" to the heat equation in the form of 
(28) and (29). The precise statement follows. 

Theorem 2. Suppose that each of a(t) , b(t) is analytic for t^O. Then 



+ V (oft) K(t)1 <n) xfl 

TTis g n 2n I 

V [a(t) K(t)) ( n) X 2n + 1 


9 


2n .' 


1 

"26 


(2n+ 1) ! 


= exp [ f (a 2 (t) - 26 b(t)) dt 


is a solution of 


v, .+ vv = 6 v 
t x xx 


such that, for t>0, -oo<x<oo , 

5.) v, v and v are continuous 

X XX 

ii. ) v(0/t) = a(t), v x (0/t) = b(t) 

iii. ) v is unique. 

As an illustration of the solution (30) , we shall give a stead/ state solution of (30), obtainable 
by letting 

v(0,t) = a Q , v x (0,t) = bg , a Q , b Q constants. 

Under such conditions we find that 

c. sinh c« x + cosh c. x 

(31) v = a n c. -T-T — r 

0 1 smn Cj x + c | cos ' 1 °2 x 

where the significant parameters c j and C 2 are given by 


C 1 " ’ 


46 % ~ 28 b 0 


0 

c 2 " ‘ 26 C 1 

which characterizes the highly nonlinear nature of the solution, 
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UTILIZATION OF INFINITE PRODUCTS 

The solution given by (30) is quite typical for Burgers' equation: one infinite series 
divided by another. There are only a very few cases where such an expression can be 
analyzed without further approximations. Therefore/ in this section we shall discuss 
how to improve such approximations and how/ in many cases, to obtain a single series 
representation for the solution v . 

Our approach is going to utilize a simple property of infinite products. Thus, let 


= TT (i+f n M)) 


with the condition 


n = 1 


f" 2 (x,t) <oo 


to insure that the product converges. Let us note that the logarithmic derivative of 
(32) is given by 


h. »,*> _ f, l f n M> l, 

h(x,t) ” 1 + f (*/ 


n = 1 


This can be seen either by talcing the logarithm of both sides of (32) and differentiating, 
or by a factor by factor differentiation of the right hand side of (32), followed by division . 

The functional form of the left h and side of (33) is exactly that of the solution of Burgers' 
equation. The question we have to discuss is therefore the following: how can one 
represent the solution of the heat equation as an infinite product? In answering this 
question we have to remember that only positive solutions of the heat equation need be 
considered (see Theorem 1). 

Let us assume that in the process of solving Burgers' equation, we reached the point where 
we already have the solution of the associated heat equation; let this be hfc,t) . To 
find v, we now have to form the quotient 


v * -28 


h x (x,t) 

h(x,t) 


*57 


However, since h(x,t)^> 0, in many cases it is possible ro write h in the form 


(34) 


h(x,t) = exp [ ]C v n (x,t) ] 


where the series (which can be finite or infinite) is a real one. In such a case, our 
infinite product is 


(35) 


fr ft 

„ = « ' n ' n = 0 


h = d 

and the solution v is given by 


(36) 


00 

= -® £ b V n (x ' ,) 

n = 0 


The situation is not so simple in the general case. Nevertheless, the approach through 
infinite products is possible in many other cases. If the solution h of the associated 
heat equation can be factored at all, then the method applies. Also, in many cases, it is 
much easier - and one obtains much better results - if instead of trying to approximate 
the quotient of two series, one approximates the solution of the heat equation . This can 
be done, for instance, by using the approximation theorem of Weierstrass, according to 
which a function, analytic on an interval, can be approximated to any degree of accuracy 
- on that interval - by a polynomial . This polynomial then can be factored and the 
method applied. 

In addition, one can use Hadamard's factorization theorem for the direct splitting of h 
into its factors; although, to be sure, it does not apply in every case. 

Finally, let us note that one can tackle this problem in yetanother way. To see this, 
we begin the solution of our problem all over again . 

We wish to find a simple series solution of 


We shall let v - 

(37) 


v + vv = 6v 
t x xx 


-26 and integrate: 


♦, - *< - *t, 


XX 
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Let now 


= In £ v (x,t) . 


Then, in abbreviated form, 


L. 6 Fx 

e 6 It 


E t - *E, 


< 39 > E JJ- »„$<»♦) “ 5 E^2 v n ( * ,f) 

3x 

We can now proceed to find the solution of (39) in terms of easily factorizable eigen- 
functions v^; use this in (38) to arrive at 


9 = In £v 0</t) = £ln f (x,t) 


and obtain the solution v as 


[f 0 </t) 1 . 


■ -* ♦* ■ 
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ABSTRACT 

N67-24275 

Radio astronomical observations made from the surface of the 
Earth are restricted at the high frequency end of the RF spectrum by the 
troposphere and at the low frequency end by the ionosphere. Artificial 
Earth satellites provide the means for eliminating these restrictions, but 
low frequency studies require large antennas, so a large stable platform 
is desirable. The surface of the Moon meets this requirement. 

In order to study the sources (Sun, Jupiter, radio stars, cosmic 
background, and perhaps Earth and Saturn) over the frequency range from 
30 kHz to 10 MHz, the radiometer should be able to measure intensity, 
apparent position, apparent size, temporal variation, and polarization. 

The order of experiments will most likely be determined by the antenna 
system complexity and size. Early experiments on the lunar surface 

probably will employ simple broadband "pencil-beam" antennas. Man ! 

I 

will be necessary for antenna installation and initial calibration. j 

t 

i 
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INTRODUCTION 


Over the past three decades a vast amount of astronomical 
information has been gained from observations in the radio frequency 
portion of the electromagnetic spectrum. Earth-bound investigations 
are restricted at the high frequency end of the RF spectrum by the 
troposphere and at the low frequency end by the ionosphere. Artificial 
Earth satellites provide the means for eliminating these restrictions, 
but low frequency studies require large antennas, so a large stable 
platform is desirable. 

In general, the scientific objectives of low frequency (30 kHz - 
10 MHz) radio astronomy in space will be to make observations of all 
possible sources of emission and to derive from these data new knowledge 
concerning source mechanisms and factors influencing the propagation of 
the radio waves. The known sources of low frequency nonthermal radio 
emission are the Sun, Jupiter, the discrete galactic and extragalactic 
sources, and the cosmic background. Perhaps Saturn and Earth will 
also fall in this category. The types of measurements to be considered 
are 1) calibrated intensity, 2) apparent position, 3) apparent size, 

4) polarization, and 5) temporal variation. These studies will require 
large antennas, and man will participate in their erection. The order of 
experiments will be determined by the antenna system complexity and 
size. 

Such investigations should result in a better understanding of the 
magnetospheric systems of Jupiter, Earth, and perhaps Saturn; of the 
distribution of plasma clouds in interplanetary space; and of the structure 
of our galaxy. Very likely many of the results will be unsuspected. 

The author wishes to acknowledge the suggestions and encourage- 
ment of Dr s. A. G. Smith, G. R. Lebo, and T. D. Carr of the Department 
of Physics and Astronomy, University of Florida. Some of the material 


162 






presented herein is the result of studies performed under NASA Contract 
NASR-176 with the University of Florida concerning the development of an 
orbital receiver for the reception of low frequency radio energy from the 
planet Jupiter. 
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ATMOSPHERIC LIMITATIONS 


Tne astronomer is interested in obtaining new knowledge about 
extraterrestrial bodies. The information he obtains at ground level arrives 
in the form of particles or waves and almost none escapes interaction with 
the earth's atmospheric envelope. The difficulties involved in sorting out 
the atmospheric influences are in many cases insurmountable. So, the 
astronomer confines himself to making observations in narrow regions of 
the electromagnetic spectrum where the atmospheric interactions are 
minimal. These regions have been aptly referred to as "windows", and 
are shown in Figure 1. 


I 
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j 

i 
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Those parts of the electromagnetic spectrum which can be advan- 
tageously examined from beyond the sea of gases which encircle our planet 
are, in order of increasing wavelength, cosmic ray, gamma ray, x-ray, 
ultraviolet, visible, infrared, and radio. In the radio region of the spec- 
trum, earth-based measurements are confined to the window between about 
5 MHz and 60, 000 MHz. At the low frequency end of this window, the iono- 
sphere becomes opaque to radio wave transmission. The critical frequency 
for normal incidence, or the plasma frequency of the most dense layer of 
electrons,* depends on the time of day, season, sunspot cycle and other 

variables, but usually ranges from 3 to 12 MHz. Figures 2 and 3 illustrate 
respectively the variation of plasma frequency with altitude in the ionosphere 
and in the magnetosphere. The high frequency boundary of the radio astron- 
omy window is caused by tropospheric absorption. This absorption is as- 
sociated with the rotational energy of the oxygen and water vapor molecules 
and occurs strongly in certain frequency bands. Between these bands, radio 
noise is a serious problem. Even in the spectral region where the earth's 
atmosphere is transparent to radio waves, there are other effects which are 
a nuisance to radio astronomers. The troposphere and the ionosphere cause 
refraction of radio waves resulting in an apparent displacement of the source 
Drifting ionospheric irregularities or clouds of ionization produce scintil- 
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Figure 1. Atmospheric Windows in the Electromagnetic Spectrum (Aller) 
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Frequency as a Function of Altitude from Electron Density Profiles 
Ionosphere at Mid Latitudes (Reference 1) 











Altitude 



Figure 3 . Plasma Frequency as a Function of Altitude from the 
Equatorial Profile of Magnetospherlc Electron Density 
(Reference 2) 167 






lations or rapid variations in the intensity of received signals. The patchy 
ionosphere acts like an optical grating, producing a diffraction pattern on 
the ground. The smaller the angular size of the source, the more pronounced 
are the focusing and interference phenomena. Electrical storms produce 
radio static and, at low frequencies, these natural disturbances as well as 
man-made interference are reflected off the ionosphere and compete with 
the desired signals. For these reasons the lowest frequencies are acces- 
sible only between midnight and dawrr when the electron density in the iono- 
sphere is at a minimum. By judicious choice of the observing station 
location - high altitude and low humidity - atmospheric absorption effects 
can be reduced considerably. 

LOW FREQUENCY SPECTRUM 

The present discussion will concern only one part of the radio 
region of the spectrum - those electromagnetic waves having frequencies 
between 30 kHz and 10 MHz, or wavelengths between 30 and 10, 000 meters. 
This part of the radio spectrum will be referred to as the low frequency 
region. Very low frequencies (f < 30 kHz) and high frequencies (f > 10 MHz) 
are outside the scope of the present endeavor. 

The atmospheric limitations at low frequencies can be avoided by 
making measurements from high earth orbits or from the moon. With no 
ionosphere to contend with, observations can be performed on down into 
the kHz range. At frequencies below the critical frequency in the F 2 region 
of the ionosphere, observations can be made neglecting natural electrical 
disturbances, associated with thunderstorms, and man-made interference, 
for these will be confined in the ionosphere -earth cavity. 

SOURCES 

The sources of low frequency radio emission that we know about 
are the following: the sun, Jupiter, "radio stars", and the cosmic back- 
ground. There are reasons to believe that Saturn might also belong in 
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this group. The radiation belts which encircle the earth probably provide 
magnetic Bremsstrahlung emissions which also could be detected from an 
orbital or lunar radio astronomy package. 

Sun 

The radio spectrum of the sun is one of enormous complexity with 
at least five different types of transient outbursts in addition to the steady 
thermal emission. Study of these outbursts has increased our knowledge 
of the solar atmosphere and, in particular, has led to a better understand- 
ing of centers of activity. Figure 4 shows the dynamic spectra of solar 
radio bursts - electromagnetic phenomena associated with chromospheric 
flares. 

Type I bursts which accompany noise storms are not shown on 
Figure 4. These narrow-band (1% of center frequency), short-duration 
(few tenths of a second) bursts occur in the frequency range 50 - 300 MHz 
and are found superimposed on a broad continuum. All of the noise storm 
radiation, continuum and bursts, is strongly circularly polarized. Its 
most notable characteristic is its long duration; a typical storm lasts 

t 

about one day. 

Type II bursts are the most intense. They are characterized by 

* j 

-I one or mors narrow bands about 50 MHz wide which drift slowly from high 

to low frequencies. When two Type II bursts occur together, as in Figure 4, 
the features are duplicated in two frequency ranges suggesting a funda- 
mental and second harmonic. Type II bursts are rarer, but longer lasting 
than Type III bursts. Some time structure has been detected. 

v! 

M Type III bursts last only a few seconds, are very intense, and often 

$ 

occur in groups. The individual burst appears as a sharp line tilted slightly 
towards the right. The emission at each moment occurs in a narrow fre- 
quency band some MHz wide, which quickly drifts toward long wavelengths. 
Both Type II and III bursts are associated with fast moving distrubances in 

'1 

?! the corona. 
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Type IV radiation has a broad spectrum and great steadiness. This 
radiation, called continuum because of its wide frequency spread, follows 
Type II bursts and may persist for several hours. It is often circularly 
polarized. 

Type V bursts follow Type III. This is also continuum radiation, 
but its spectrum is not as broad as Type IV. It appears only at wavelengths 
greater than a few meters, and no fine structure can be seen. 

Since the solar radio spectrum is continuously changing, attempts 
to extend measurements into the low frequency region which borders the 
top of Figure 4 should be designed to fully reveal the characteristics of the 
bursts in the frequency time domain. The experiment should be capable 
of tuning over a wide range of frequencies in a fraction of a second and 
recording the intensity of the radiation at each frequency. 

Occultation experiments - monitoring the emission from discrete 
sources as they pass behind the solar corona - can provide information on 
the heliocentric distribution of plasma in interplanetary space. Low fre- 
quencies would be especially advantageous here since it is believed that 
the scattering effects of plasma clouds scale as \ 2 . 

Jupiter 

Excepting the sun, Jupiter is the only confirmed nonthermal source 
of radio frequency emission in the solar system. The sporadic bursts at 
dekameter wavelengths are strongly correlated with Jovian longitude and 
with the position of the satellite lo indicating that the sources of radiation 
are linked to the disk of the planet. Whether the mechanism responsible 
is a surface phenomenon or occurs in the exosphere and is coupled to the 
planet through its magnetic field has not yet been determined. Another 
nonthermal component of emission which occurs in the microwave part of 
the spectrum is believed to originate in Jovian Van Allen belts where 
trapped electrons radiate by the synchrotron process. Figure 5 shows 
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the spectral distribution of the average planet-wide flux density from 
Jupiter? The spectral index varies between -5 and -8 over the observable 
dekamef ’ic range, suggesting that the bulk of the radiation may actually 
lie in the low-frequency tail, which is unobservable from the ground 
because of ionospheric cutoff. Furthermore, the apparent structure of 
the Jovian radio sources changes radically at frequencies below 15 MHz, 
and the pulse structure of the emission is unquestionably altered by the 
terrestrial ionosphere. 

Observations of the giant planet at frequencies below 5 MHz could 
establish the spectral distribution of the emission, i. e. , extend the curve 
of Figure 5 to the left. Such information is vital to theories concerning 
the origin of the radiation - total energy involved, particle densities 
required, and efficiency of the conversion mechanism. An examination 
of the pulse structure, unaltered by the terrestrial ionosphere could yield 
information about Jupiter's ionosphere. Dynamic spectra (frequency 
versus time plots) and polarization studies at these low frequencies might 
unravel the mysteries of the planet's magnetic field. Long-term monitoring 
would undoubtedly yield evidence of propagation effects in interplanetary 
space, correlation with sunspot activity, and increase our knowledge of 
another magnetospheric system, perhaps much like our own but on a 
larger scale. 

Saturn 

There have been reports of possible emissions from Saturn at 
4 5 

low frequencies.’ These do not, at present, constitute hard evidence. 

Since modern theories tend to associate a planet's magnetic field with its 
rate of rotation, Saturn might be expected to emit because of its similarity 
to Jupiter in this respect. The fact that no confirmed nonthermal radiation 
from the ringed planet has been detected could be the fault of listening at 
the wrong frequencies. Here again, monitoring at low frequencies, un- 
hindered by an ionosphere and atmospheric and man-made interference, 
could solve the riddle. 
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Figure 5. Spectral Distribution of the Average Planet-Wide 
Flux Density from Jupiter 


Radio Stars 


There are numerous discrete sources on the celestial sphere, and 
these have been dubbed "radio stars" because of their small extent. Actually 
they are not stars, but nebulae, and their emission is continuous in time. 
They can be classed as galactic or extragalactic, thermal or nonthermal. 
Only the discrete sources of our own galaxy exhibit typical thermal emission 
spectra. These sources (galactic, thermal) are thought to be clouds of 
ionized gas - the coronas surrounding bright, hot stars in the spiral arms, 
and are referred to as H-II sources. A typical example is the source 
Cygnus X. All of the galactic, nonthermal sources are believed to be rem- 
nants of supernovae - catastrophic explosions of stars. The Crab Nebula is 
an example and emission is by the synchrotron process. All of the extra- 
galactic discrete sources, which constitute the vast majority of "radio 
stars", are characterized by nonthermal spectra, i. e. , the radiation in- 
creases in intensity with wavelength. Again, the synchrotron mechanism 
offers a reasonable explanation of this emission. These extragalactic 
sources which have an isotropic distribution in space are classified as 
"normal" and "peculiar". The normal sources are bright, nearly spiral 
galaxies with radio outputs of approximately 10 28 kw, i. e. , the same total 
power as the Milky Way. The classic example is the Great Nebula in 
Andromeda. The peculiar sources are much more intense. For example, 
the source Cygnus A, believed to be two galaxies in collision, has a radio 
output of 10 34 kw. 

All of the discrete radio sources emit a fairly smooth continuum 
of radio frequency energy. Measurements of the intensity of these sources 
at low frequencies will provide additional spectral evidence which has a 
bearing' on the stability of the emission mechanism and on the energy spec- 
trum of the particles responsible. Surveys at frequencies below 5 MHz 
showing the positions in space of the discrete sources would be valuable 
additions to the information now in hand. This mapping, along with indi- 
vidual source studies to ascertain angular dimensions and intensity 
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distributions, will require high resolution and narrow beamwidth. Such 
investigations are vital to attempts at optical identification of the radio 
. . sources. 

i , 

| Galactic Background 

i i 

By background radiation, we mean the emission that has a continu- 
ous spatial distribution as opposed to the discrete sources. Most of this 
radiation comes from our own galaxy. and two types can be distinguished, 
the galactic continuum and the hydrogen-line emission. Since the emphasis 
here is on the scientific worth of extending the radio astronomy spectrum to 
^ lower frequencies, the 1420 MHz emission from neutral hydrogen will not 

be discussed. 

Galactic continuum rveys or maps of radiation intensity show the 
following: average intensity varies over the celestial sphere but appears 
to be constant for a given direction; intensity is greater in the galactic 
plane and is a maximum toward the center of the galaxy; surveys at dif- 
ferent frequencies indicate that both thermal and nonthermal processes 
are operative. The bright band of emission in the galactic plane is thermal 

f 

\ emission from clouds of ionized hydrogen (H-U -regions). The extended 

component of the background emission is synchrotron radiation by cosmic 
iij ray electrons in the galactic corona which is spherical about the galactic 

f- 

- center. This emission is most intense at lower frequencies as indicated 

1 

in Figure 6. There is synchrotron emission from the galactic disk also, 
but it is absorbed by H-H regions. At long wavelengths (dekameter and 
| longer) where the nonthermal emission predominates, galactic observa- 

tions permit a study of the structure of these regions in absorption. 

-i 

The points in Figure 6 show the measurements of the galactic 
f background noise. One group of experimenters believes in a rather 

? drastic decrease in intensity between 3 and 1 MHz. Another group sees 

a leveling off in this region. 
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Figure 6. Observed Cosmic Noise Intensities Below 10 MHz 


The intensity values reported by Benediktov, et al. at 0. 7 MHz are 
an order of magnitude smaller than the values reported by Huguenin, et al. 
Perhaps what is being sampled is earth noise. This might be expected to 
possess large, order of magnitude, variations whereas the cosmic back- 
ground would not. Measurements made outside of the magnetosphere could 
perhaps resolve this enigma and give a true reading of the intensity of cosmic 
radio emission below 2 MHz. 

The shape of the spectrum of the galactic continuum at low fre- 
quencies has consequences on other lines of endeavor. If the intensity of 
the emission levels off or increases, it will be difficult to detect other 
sources above this interference. On the other hand, if the emission falls 
off, the detection of weak signals in this part of the spectrum will be pos- 
sib’e. Whatever the shape, the spectrum definition will provide input to 
questions of energetics - populations of cosmic ray electrons, and acceler- 
ating mechanisms. The spatial distribution of galactic radiation at low 
frequencies will undoubtedly lead to a clearer delineation of the spiral 
structure of the Milky Way. 

Earth and Interplanetary Medium 

The Van Allen radiation zones probably provide magnetic 
Bremsstrahlung emissions which could be detected from the moon. Oc- 
cupation experiments with the aim of detecting propagation effects caused 
by the interposition of the earth's ionosphere or magnetosphere between 
the source and the receiver can yield valuable data on the terrestrial 
environment. 

It is presumed that interplanetary space contains filaments or 
blobs of plasma that have been spewed out of the sun. From the moon it 
should be possible to ascertain the distribution of these solar plasma clouds 
by monitoring celestial sources at low frequencies over a time span of one 
year and scrutinizing the data for propagation effects such as the scattering 
of radio waves by irregularities in electron density. The moon if well 


177 



suited for such an investigation for, besides providing a large platform 
for low frequency observations, it is outside of the terrestrial magneto- 
sphere a majority of the time and thus enables a distinction between propa- 
gation effects in interplanetary space and those operative at the magneto- 
spheric boundary. 

TYPES OF MEASUREMENTS 

* 

Certain types of measurements are possible utilizing the tech- 
niques of radio astronomy. These are intensity, apparei.l; source position, j 

j 

apparent source size, temporal variation and polarization. Intensity I 

measurements at fixed frequencies will extend source spectra into the ; 

f 

region inaccessible to earth-based radio telescopes. Swept frequency 
measurements will yield dynamic spectra,' i. e. , how the frequency of i 

emission drifts with time. Measurements of source positions or "mapping" 
will be a definite asset to studies of galactic and extragalactic structure. 

Source size is a clue to the mechanism of emission. Synoptic observations 
will provide data for time variation investigations, correlation with other 
space parameters, and cause and effect analyses. From this standpoint, 
unattended operation of the experiments would be desirable to provide a 
long time span. Polarization data is necessary to infer the character of 
magnetic fields existing at the source and certain propagation effects. 

I 

INSTRUMENTATION 

i 

A radio astronomy package designed for the space environment 
should include antenna, receiver, calibration source, and programming 
circuitry. Figure 7 shows thp relationship of these essential components. 

Power supply and telemetry are considered supporting equipment. The 
calibration source and programming circuitry are standard items. Various j 

types of receivers have been developed for low frequency experiments in 
space. The greatest effort will be required in the field of antenna design. 

In the detection of relatively weak and small sources, very large antennas 
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Figure 7. Essential Components of a Radio Telescope 
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will be required. There are two reasons for this. First, the effective 
area must be sufficiently large to capture enough energy to be measurable, 

Second, because of diffraction, the angular resolution of a radio telescope 
improves with each increase in the ratio of the aperture to the wavelength. 

When the antenna is swept across a point radio source, radiation is re- 
ceived over an angle of approximately X/D radians, where D is the effective 
width of the antenna. Two small sources will appear as a single source if 
their angular separation is much less than X/D, and they cannot be resolved. 

i 

At low frequencies, cosmic noise exceeds the receiver noise and j 

f 

55 

becomes the factor which limits the sensitivity of a radio telescope. [ 

I 

Figure 6 shows the increase in cosmic noise intensity with decreasing j 

frequency. High resolving power (paid for in terms of large antenna size) 
is necessary to insure that the integrated cosmic noise does not drown j 

out the desired signal from some discrete source. t 

,* i 

TYPICAL EXPERIMENTS 

Advances in science usually await advances in instrumentation. I 

It is only natural that the capability to leave the earth and its atmospheric ; 

envelope be exploited to extend the frontiers of science. Low frequency ! 

i 

radio observations of extraterrestrial sources can be performed from 

above the ionosphere, that is, from earth-orbiting platforms at altitudes 

greatei than 2000 km and from lunar bases. Earth- orbiting platforms are 

less expensive, thus early efforts will be in this direction. However, for 

even moderate resolution at these long wavelengths, the antenna arrays 

become quite large. There are limitations to the size and stability of I 

antenna systems deployed in space. When these limitations are too severe 

for the experiment being considered, then the lunar base becomes attractive. 

Earth -Orbiting Observatories 

• ' 

Low frequency radio experiments are planned for the Orbiting 
Geophysical Observatory (OGO) and the Radio Astronomy Explorer (RAE) 
series. Both are unmanned, so the experiments are entirely automatic. 

* 1 

180 | 




The OGO spacecraft is designed to accommodate about twenty different 
experiments on each flight. Two classes of orbits are planned - one 
circular, polar w . apogee of 1000 km, the other elliptical, 31 ° inclination, 
with apogee of 133 , 000 km. For the reception of low frequency radio signals 
from extraterrestrial sources, the latter class of orbits is more desirable 
since both the outer ionospheric electron density and the magnetospheric 
electron density decrease with altitude. Experimenters at the University of 
Michigan intend to measure the dynamic spectrum of solar noise bursts 
using a swept frequency radiometer and a dipole antenna. The radio 
astronomy group at the University of Florida hopes to determine the 
spectrum and pulse structure of Jupiter's low frequency emission on one 
of the later OGO flights. The receiving system for this experiment is a 
four-channel superheterodyne total power radiometer and the antenna is a 
240-foot dipole. Figure 8 shows the antenna patterns at the four operating 
frequencies: 4* 2, 1, and 0.5 MHz. This simplest type of antenna provides 
very little directivity. 

The Radio Astronomy Explorer satellites will go one step further 
in antenna directivity. Using two 750-foot VEE antennas, observations of 
terrestrial and celestial noise will be made from a 6000-km circular orbit. 
At this altitude, the satellite will be within the magnetosphere, but obser- f 
vations should still extend down to 0. 5 MHz, the local plasma frequency. 
Figure 9 shows the satellite configuration. 

The high-altitude Apollo flights (synchronous orbit) provide the 
capability for the next advance in low frequency radio astronomy. The 
addition of man to the experiment makes it possible to employ large direc- 
tive antenna arrays and guide’ them. By extravehicular maneuvers, the 
astronaut can assemble these more complex antennas and point them in 
preferred directions using the spacecraft control system. Figure 10 shows 
a large log-periodic array sensitive over a frequency decade. The availa- 
bility of man also permits maintenance and alteration of system parameters 
based on the results obtained, thus making the experiment more reliable 
and more flexible. 





Figure 8. Antenna Patterns for a Free Space Dipole of Length L 






Unfurlable Log Periodic Antenna for the Frequency Range 0.7 - 7.0 MHz 







Lunar Observatory 


Knowing the advantages that a high altitude manned orbiting plat- 
form provides, why attempt radio astronomical measurements from the 
moon? The lunar surface affords the advantages of a large stable platform. 
Very large antennas can be deployed from a roving vehicle and orientation 
of the antenna beam can be accomplished electronically. Once the lunar 
base is a reality, maintenance will be less unfamiliar under lunar gravity 
than under the zero-g conditions of an earth-orbiting spacecraft. For 
those experiments demanding high resolution, the lunar horizon can be 
used as an occulting edge, allowing the precise positioning of the sources 
of low frequency radio emission. On earth, the atmosphere rules out 
this possibility, and the guidance of an artificial occulting disc in proximity 
to an earth- orbiting spacecraft would be extremely difficult. 

In order to design a flexible system with research capabilities in 
many of the areas mentioned earlier-, the characteristics of the essential 
elements of a radio telescope must be considered. For the antenna, these 
characteristics are bandwidth, beamwidth, gain, and steerability. Since 
we are interested in a wide range of frequencies (30 kHz - 10 MHz) the 
antenna should be broadband. This can be accomplished by constructing 
the array of broadband elements. The beamwidth should be reasonably 
narrow, allowing the resolution of closely spaced sources. This is 
determined by the length of the array. Narrow beam implies high gain 
or the ability to detect weak sources. The antenna beam should be 
steerable in elevation so that a large portion of the celestial sphere can 
be monitored. For a broadband array, this can be accomplished using 
delay line techniques. To measure intensity, it is best to use a large 
but simple antenna; however, such an antenna does not give the position 
and diameter with high accuracy. Interferometers are best from this 
standpoint, but they cannot make reliable measurements of intensity due 
to the overlap of the interference fringes of neighboring sources. Single 
beam antennas provide the best plots of distributed radiation, i. e. , a 
general survey of the sources and the background radiation. 


There are many different types of radio receivers, and these will 
not be discussed here. Briefly, fixed frequency channels are better for 
absolute intensity and pulse structure measurements, while swept fre- 
quency receivers provide information in the frequency-time domain and 
are better for locating sources of interference. 

The radio telescope shown in Figure 1 1 is a compromise which 
will achieve many of the objectives mentioned earlier. It is designed to 
operate over the frequency decade 0. 2 - 2. 0 MHz. The resolution of this 
array is 10° at 0. 2 MHz. At higher frequencies the "pencil beam" be- 
comes narrower, but there are sidclobes. At frequencies below 0. 2 MHz 
the beamwidth is greater than 10°. A multichannel receiver at the center 
of the array can perform fixed and swept frequency observations. The 
east-west baseline is 10. 5 km long and contains 22 dipoles each a half 
wavelength at 0. 2 MHz (X/2 = 750 m). The north arm is 6 km long and 
contains 7 dipoles. The spacing between the dipoles can be less than 
750 m since the lunar surface is assumed to be a lossy medium. Such a 
medium cuts down on sidelobes and makes the antenna more broadbanded. 
The array has been deployed on the moon by a roving vehicle. Surface 
irregularities up to a tenth of a wavelength can be tolerated. To change 
the elevation of the beam, switches are activated which insert delay lines 
between the separate elements and the receiver. The astronaut will be 
needed primarily for the initial installation. Operation of the experiment 
will be automatic and should collect data for a period of at least one year. 
Proximity of the antenna to the Emplaced Scxentific Station will facilitate 
use of the latter's data storage, power, and telemetry systems. 

With stfch a radio telescope, we can expect to learn more about 
the structure of our galaxy, the distribution of plasma clouds in inter- 
planetary space, the magnetospheric systems of Jupiter, earth, and 
perhaps Saturn, and source mechanisms in general. Later experiments 
will involve larger, more complex antenna .systems for finer directional 




resolution and for polarization measurements. Besides the experiments 
described, there are other worthwhile areas of investigation which have 
been omitted for the sake of brevity. These include millimeter observations, 
beacon transmissions between the earth and the moon, lunar ionosphere 
sounding, point-to-point communication on the moon, and electromagnetic 
probing of the lunar surface and subsurface. 
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Our basic concept is to employ the cross correlation of two narrow 
light beams for remote sensing. The problem with most optical methods 
such as schlieren, shadowgraph systems or spectrometers is that the 
received signal is always integrated along the line of sight. However, 
we want information for a particular spot along the line of sight. To 
this purpose we adjust a second beam to intersect the first at the point 
of interest. The integration along the two beams is then eliminated by 
a cross correlation of the detected fluctuations. This paper indicates 
how this new concept of combining standard spectroscopy with a statistical 
cross correlation analysis might be used for dynamic and kinematic studies 
such as the formation and motion of water vapor or dust clouds, the 
dynamics and production of the ozone layers, and the generation and 
motion of radiation belts. 


The concept has been demonstrated successfully by measuring the 
turbulent fluctuations in subsonic jets. The main problems in AAP 
applications are formulated in terms of questions, such that the 
background information becomes apparent which is needed to establish 
the feasibility of particular experiments. 



INTRODUCTION 


The cr jss correlation of optical signals has been used for remote 
sensing in jet shear layers [l]. Concentration fluctuations of a water 
vapor fog could be resolved with good accuracy in space and time [2]. 

In principle, the method should also work on an atmospheric scale. 
Experimental investigations about the generation and motion of atmospheric 
constituents such as clouds, nuclear debris, ozone, radiation belts, etc., 
offer an exciting challenge. We will now give a simplified version of 
our new remote sensing concept and outline some of the problems, which 
one would have to overcome in atmospheric applications. 


THE CROSSED BEAM CONCEPT 

Remote sensing of distant light sources like planets is based on a 
fundamental theorem of geometrical optics, which states that the spectral 
radiance, I, of a source and its image are equal, if one neglects light 
extinction and emission along the line of sight. The spectral radiance 
of a distant source is therefore accessible through the spectral radiance 
of its image on the exit slit of a monochromator, as shown in fig. 1. 
Source area, AA, and wave length interval, AT, are set by the optical 
filter. The solid angle, Ai>, into which the exit slit radiates, is set 
by the photodetector. Therefore, the received radiative power is directly 
proportional to the spectral radiance, I, of the distant source: 
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However, if the line of sight goes through the atmosphere, the detected 
spectral radiance Id(t, T) is changed due to the light extinction or 
emission in the detector's field of view. In the two limiting cases of 
dominant extinction or dominant emission, the change in spectral radiance 
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may be expressed by integrating the extinction or emission coefficient 


1 cSX . . 

K(x, t, 7 ) = - s— = K(x*, 70 + k(x>, t, 7 ). (3) 
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In astronomy and optical communication, these effects are a nuisance. 

From an atmospheric science point of view, they can be used to study the 
thermodynamic properties of the atmosphere, since the extinction coeffi- 
cient K is uniquely determined by the thermodynamic properties. 

The problem with local measurement of atmospheric properties is how 
to retrieve the local information, that is buried in the integrated 
signal. A solution is only possible by repeating the remote sensing 
as a function of the remaining independent variables, such that some 
mathematical transformation can be used to extract the local extinction 
coefficient . 

Several authors [3-8] used the remaining space variables y and z 
and had to assume axial symmetry. This assumption can be made in the 
study of single rocket exhaust and flames but not in the atmosphere. 

Another group [9 - 15] used the dependence on optical wave length 7 . 

They subdivide the line of sight into as many zones as there are isolated 
molecular bands, and need some advance knowledge about the spatial distri- 
bution of temperatures and concentrations. The method worked on a package 
of controlled burners [16]. In the atmosphere, a larger number of isolated 
bands which are free from scattering, and the advance knowledge are hard 
to find. However, the "zonal approximation" should be considered as a 
backup of our cross correlation method, since it does give mean value 
information, which we cannot obtain. 

We use the time aspect to retrieve the local information. Consider 
the experimental arrangement on fig. 2. Two lines of sight intersect in 
an atmosphere-fixed point. These lines define a plane which we shall 
call "turbulence wave front" and a normal in the intersection point 
(x, y, z) which is called the wave normal. We then measure the spectral 
radiance l(t, 7 ) (radiative power/solid angle, area, and wave length 
interval) along both fields of view, split off the time average signal, I, 
cross correlate the fluctuation!!, i(t), and calculate the quantity, G, as 
described in fig. 2. Let us assume that the two detecting telescopes 
have been aligned to the y and z axis of a Cartesian, atmosphere-fixed 
coordinate system. Each line of sight is denoted by an y or z subscript 
and each point on the line of sight by the distances q or f relative to 
the intersection point. The quantity G is then equal to 


(4) 
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The interpretation of the measured quantity G follows by relating 
it to the local fluctuations, k, of the extinction coefficient. We have 
shown that G is equal to an area integral over the space time correlation 
function [ 17] , 


-^-CO 


g(s, t, a>-> = 


k(x, y, z+t, t, Ax) k(x+f., y+r), z, t+x, A 2 ) d^df, (5) 


provided that the following assumptions are met: 

1. The number of statistically independent atmospheric disturbances 
is small enough to avoid a cancellation of the integrated 
fluctuation. 

2. The integrated fluctuation of the extinction coefficient is small 
enough to permit linearization. 

3. The light source fluctuations and the atmospheric fluctuations 
are statistically independent. 

4. The fluctuations of the two light sources are uncorrelated. 

5. Either extinction or emission dominates the change of the 
detected radiative power. 

The first assumption is not very restrictive since the integrated mean 
value would still be large. Furthermore, the local fluctuation, k, does 
not have to be small; only the integral over statistically independent 
fluctuations is assumed to be linearized. In the very few cases where 
these weak restrictions are not met, one can possibly choose another wave 
length, A, where the second assumption is met. The other assumptions can 
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mostly be met by a suitable experimental arrangement. Therefore, the 
result of equation 5 appears to be sufficiently general for atmospheric 
applications. 

Two important conclusions can be derived from equation 5: 

1. For a pair of statistically independent light sources, the source 
characteristics do not have to be known. The method would work for any 
source, such as ground beacon;-:, the sun, or reflected sunlight. 

2. The random nature of turbulent fluctuations assures that integrand 
drops to zero over distances t>.hich are comparable to a typical eddy size. 
This explains, qualitatively, why the cross correlation partially eliminates 
the usual optical integration along a line of sight. 

According to equation 5, the crossed beam correlation works whenever 
the wanted signal is common to both lines of sight. The way in which the 
wanted signal is then pulled out of the integrated signal may be described 
best by splitting each of the integrated signals in two parts. The first 
describes the contribution of the "correlated volume" where the integrand 
of equation 5 does not vanish. The second part describes the rest of the 
integrated signal. Multiplying the fluctuations of the two integrated 
signals then leads to four products. Three of these change between 
positive and negative values in a random fashion. If we now average 
these products over time, their mean value will go to zero, whereas the 
deviation from the mean will increase with the square root of time. 

The fourth product describes the contribution of the correlated 
volume, which is common to both lines of sight. The associated product 
is always positive and its summation over time should increase linearly 
with time. Therefore, it will be the dominant one, if one goes on adding 
long enough. Even small contributions at the beam intersection can be 
pulled out provided (1) the integrated signals show detectable fluctuations 
and (2) the combination of light scarce fluctuations, shot noise and 
instrument noise is not orders of magnitudes larger than the root mean 
square value of the integrated signal. 


MEASURABLE TURBULENCE PARAMETERS 

The optical integration is restricted to the wave front. Along the 
wave normal, no integration takes place and local information inside a 
correlated volume may be obtained by repeating the crossed beam experiment 
for several beam separations along the normal as shown in fig. 3. The 
space separation, f; describes the minimum beam distance, which defines 
the wave normal. Also a time separation is introduced electronically by 
a time delay unit. As will be discussed later, the measured signal, G, 
then approximates a two-point product mean value, which would have been 
measured by point probes on the wave normal. 


195 



G( £ ) ^ f 


i z (t) iy(t+r) (^^=0 


G( 0, 0, *)£ (i z (t) i y (t)) 0 (Iy)^ 0 


k(x, y, z, t, 7Q k<x+P., y, z, t-Hr, 70 
k(x, y, z, t, 7\) k(x, y, z, t, 70 


(6) 


The position of these imaginary point probes is indicated on fig. 3 
by two dots. Tims, all turbulent properties, which are commonly derived 
from two-point product mean values, may also be approximated with the 
crossed beam method. 

The local integral scale of turbulence is defined by averaging the 
area under the covariance curve G(|, t= 0) and its graphical evaluation 
is illustrated on fig. 3. 

Plotting the beam correlation coefficient against the time separation 
instead of the space separation, fig. 4, allows bulk convection speeds 
to be read. For a given space separation, the correlation will reach a 
maximum at the time that the turbulent eddies need to travel from the 
upstream Co the downstream beam. The time lag, which is indicated 
by the maximum is therefore a direct measure of convection speeds. 

Power spectra and cross power spectra are defined by a Fourier 
transform of the measured autocorrelation, g ** 0, and cross correlation 
| j* 0, curves. They may be calculated using the digital computer programs, 
developed for the analysis of mechanical vibrations [18]. 

A study of. the envelope to the time correlations, r(r). provides a 
measure of the temporal rate of decay of turbulent eddies as seen by an 
observer moving with the convection speed. Therefore, the decay of this 
correlation curve to 1/e of its original value is used to define the 
average eddy lifetime. The Fourier transform of this envelope would 
Indicate the dominant frequencies which are felt by the moving observer. 

Sisnmarlzing the graphical evaluations sketched in figures 3 and 4, 
ve find that the crossed beam method remotely senses Integral scales of 
turbulence, convection speeds, eddy lifetimes and spectra for arbitrarily 
chosen atmosphere-fixed points at Which a beam intersection can be achieved. 


196 



The above turbulent properties are only approximations, since the 
optical integration over the wave front does catch fluctuations which 
are not felt by the point probe on the wave normal. However, our experi- 
ments in a subsonic jet have shown that weighting of contributions close 
to the wave normal is sufficiently strong to yield a good spatial 
resolution. One may thus expect a good spatial resolution on a wide 
range of other turbulent motions which are of practical interest. 

The absorption (or emission) spectrum of extinction coefficients 
is commonly used to find species concentrations and temperatures of both 
stagnant gases and small scattering particles. Adjusting the optical 
wave length of the crossed beams to chosen atomic lines and/or molecular 
bands, one might therefore extend the spectroscopic methods from stagnant 
media to inhomogeneous turbulent flows and study multicomponent or two-phase 
flow phenomena. 


EXPERIMENTAL VERIFICATIONS 

The most revealing and instructive test of the crossed beam concept 
is to compare the optical approximation of point measurements with known 
hot-wire measurements. The initial portion of an axisymmetric free shear 
layer was chosen for the experiments since (1) the turbulence is aniso- 
tropic and inhomogeneous and (2) convection speeds, turbulence scales 
and eddy lifetimes are well documented in the literature. 

All measurements were taken in a subsonic (M = 0.2) air jet exhausting 
through a one-inch diameter nozzle into the atmosphere (fig. 5). Light 
extinction was achieved by spraying a small amount of liquid nitrogen 
into the settling chamber. This produces small water droplets in the 
exhaust flow which attenuate the crossed beams by scattering. Later, 
the nitrogen injection was replaced with water spray which offers a 
better control of the tracer concentration. 

A diagram of the optical and electronic hardware is also shown in 
fig. 5. A plane mirror and a spherical mirror collect the light from a 
powerful mercury arc and transmit i.t as a parallel beam through the jet. 

A similar mirror combination projects the arc image on the entrance slit 
of a grating spectrograph. The aperture stop in front of the spectrograph 
is set at a beam diameter of 2 mm. The photomultiplier then scans the 
first order spectrum of the grating. In this way, the wave length interval 
A?^ and the wavelength ? are adjusted by the width of the monochromator 
slit and the rotation of the grating, respect Ivcly. By exchanging light 
sources, gratings, and photodetectors, the system is able to cover the 
spectrum from the vacuum ultraviolet (). * 1200A) to the infrared (, k * 25 ?r). 
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The above light source and detector arrangement is used to generate 
two beams, one in a horizontal and one in a vertical direction (fig. 5), 

A lathe bed provides a mechanical support which allows us to align the 
beams parallel to the nozzle exit plane. A beam separation iu the stream- 
wise direction gives the longitudinal turbulance scales and convection 
speeds. Transport of the whole system allows us to repeat the measure- 
ments for all axial and radial traverses of the jet. 

i ' 
i ' 

A typical sample of optically measured two-point product mean values is ; 

shown in fig. 6. The beam intersection was adjusted to the center of 
the free shear 1 yer three exit diameters downstream of the nozzle exit. 

The general appearance of the temporal correlation functions is indistin- ; 

guishable from published hot-wire results. The turbulent convection 
process is clearly indicated by the displaced maxima for finite beam 

separation. The monotonic decay of the envelope shows the decay of the j 

moving eddies. Measuring the time delay for which this curve falls to 
1/e of its initial value yields a value of 880 microseconds. This is 

in excellent agreement with the hot-wire data of Ref. 19 from which we j 

estimate that the value of the time scale for our experimental conditions 
should be 920 microseconds. 

Fig. 7 shows the same experiment under different conditions. The 
initial intersection point of the beams was six diameters from the jet 
exit, (i.e., at x/D = 6.0) instead of at x/n - 3.0. Instead of using 
liquid nitrogen to generate a water droplet fog tracer, a small nozzle 
was mounted in the jet settling tank which sprayed wa**er droplets into 
the flow to form the tracer fog. The curves do appear to fall onto or 
close to an envelope as required. Unfortunately, the range of beam 
separations used is not sufficient to meri the extrapolation of the 
curve to 1/e of its initial value. There t.ore, the moving axis auto- 
correlation curve has been obtained from the hot-wire data of Ref. 19 and 
is shown as the dotted line superimposed on the measured cross correlation 
curved 

The spatial resolution of the crossed beam method eems yery encourag- 
ing, considering that the envelope from the hot-wire data exactly matches 
the individual curves from the crossed beam method. For the moving frame, 
the optically measured autocorrelation agreed with the hot-wire curves, regard- 
less of whether nitrogen or water spray was used as a tracer* 

A quantitative comparison of individual cross correlation curves 
G(x)f, has been based on the indicated convection velocities. These 

velocities have been evaluated graphically following the procedure of 1 

fig. 4. To increase the accuracy, we have always used more than one 1 
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curve. The chosen set of beam separations was plotted as a function 
of the time lag values i^, which were indicated by the common tangent 
of envelope and individual cross-correlation curve. The result is a 
straight line which is shown in the upper part of fig. 7. The desired 
convection speed follows from the slope of this line. 

The results of all convection velocity determinations made using the 
crossed beam correlation system to date are summarized on fig. 8, 
where they are also compared with "hot-wire data from Reference 19. 

It is shown that in the outer portion of the shear layer (t) > 0) the 
convection velocity determinations agree within experimental accuracy. 
For T| < 0, however, the crossed beam results tend to follow the mean 
velocity profile (U/U 0 ) rather than attain the constant maximum value 
observed for hot-wire data. This was not an unexpected result. All 
convection velocity results presented in Reference 19 were taken at 
axial stations where the potential core still exists. At these posi- 
tions, the large difference between the mean velocity and convection 
velocity has not been completely explained However, Davies [ 20] has 
discussed the subject recently, and it does appear that the closer 
correspondence between mean and convection velocity observed here is 
to be expected when a measuring technique other than the hot-wire 
anemometer is being used, and/or the potential core no longer exists. 

The optically measured results in fig. 8 were only 0. 2 inch 
apart. This is much smaller than the geometrical shear layer thickness 
( 8") which is indicated by the hot-wire results. Thus, convection 
speeds have been resolved inside a correlated volume by moving the 
intersection point. 

REMOTE SENSING FROM LOW AND HIGH EARTH ORBITS 

The concept of remote sensing through cross correlation of optical 
signals is not restricted by scale and the method should, at least in 
principle, work also over atmospheric distances. Mounting the two 
photodetectors on satellites instead of on the ground or on airplanes 
has the advantages of a larger field of view, quicker scanning, and 
avoiding vibration problems. In addition to these benefits, which 
were discussed during the recent remote sensing symposium, the satellite 
detector can use reflected sunlight, whereas the ground detector is 
restricted to scattered sunlight. The orbital detector can thus use 
absorption spectroscopy for a detailed study of atmospheric and iono- 
spheric constituents, whereas the ground detector is restricted to 
refractive index fluctuations, the interpretation of which is very 
difficult. 




Sometimes the fluctuations oL r atmospheric constituents may be 
detectable with the radiative power, which is provided by reflected 
and scattered sunlight or by emission inside the detector*s field of 
view. In such a case, the optical hardware which would be needed 
for the crossed beam method is quite similar to that which is 
developed for the photographic surveys of the earth and the moon. 

Bo-th applications require a diffraction-limited telescope coupled with 
an optical lilter which selects the ultraviolet, visible, and infrared 
regions of the spectrum. The main difference consists in replacing the 
camera with a sensitive phototube and to fly two satellites instead of 
one. 


Two orbital detectors could adjust their lines of sight to intersect 
on an atmosphere-f ixed spot, by continuously turning the telescope 
during the fly-by. Remote sensing of disturbances around the inter- 
section might be feasible for both the lower and the upper atmosphere. 
Typical lines of sight are sketched on figures 9 and 10. 

We suggest to study the feasibility of remote sensing in the lower 
atmosphere with respect to the following AAP experiments: 

1, The generation and motion of nuclear debris using 
emitted radiation during the night cycle. 

2. The motion or generation of clouds, dust, and iodine 
using reflected sunlight (fig. 9). 

Remote sensing in the upper atmosphere is intriguing, if the use of 
ground beacons can be avoided. The optical path between the upper 
atmosphere and the orbital telescope is then short. Also, oxygen has 
extremely high absorption peaks in the ultraviolet, such that the 
motion of very diluted oxygen might still be detectable. We suggest 
to study the feasibility of remote sensing the upper atmosphere and 
ionosphere with respect to the following AAP experiments: 

1. The motion or generation of ozone using direct or 
scattered sunlight, fig. 5. 

2. The motion or generation of radiation belts using the 
emitted electromagnetic radiation. Both vertical and 
horizontal motions might be studied using the same 
arrangement sketched for the upper atmosphere (fig. 10). 
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The main problems with the crossed beam method are as follows: 

1. To find an atmospheric constituent that produces 
detectable fluctuations of the integrated signal 
along accessible lines of sight. 

2. To reduce light source fluctuations, shot noise 
and instrument noise to root mean square levels, 
which have the same order of magnitude as the 
integrated signal. 

3. To have the lines of sight intersect on an atmosphere- 
fixed spot for an observation time, which is sufficiently 
long to pull the local signal out of the integrated 
signal and the combined light source noise, shot noise 
and instrument noise. 


These problems are formidable and they may prohibit some of the appli- 
cations that were outlined in the previous section. However, since 
we can pick the detected radiation from any region of the electro- 
magnetic spectrum and since the orientation of the lines of sight may 
be chosen arbitrarily, we believe that there is at least one interest- 
ing AaP experiment where the above problems may be overcome through 
special test arrangements. Such an experiment presents an exciting 
challenge, since the dynamic behavior of the atmosphere, or radiation 
belts could be added to the "static" mapping on oceanic, geographic, 
and atmospheric surveys. 


Several facts about atmospheric constituents and optical AAP 
hardware must be collected before we could decide the feasibility of 
particular experiments. These facts have been identified in the pre- 
vious discussions and may be summarized by the following questions: 

1. Which constituents of the lower atmosphere, such as 
clouds, dust and nuclear debris, produce detectable 
fluctuations within a response time of 0.1 second 
along a line of sight which traverses the entire 
atmosphere? 

2. Do the usual constituents of the upper atmosphere 
and ionosphere such as ozone produce detectable fluc- 
tuations along a line of sight which traverses the 
outer portions of the upper atrosphere. 
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3. Do reflected or scattered sunlight and/or radiation 
belts put enough radiative power into the orbital 
detector's field of view to reduce the shot noise 
below the expected signal levels, considering the 
telescopes and optical filters which are developed 
for other AAP missions? 

4. Does the continuous motion of the detector's field of 
view lead to fluctuations of the reflected sunlight that 
are comparabLe to the signal fluctuations? 

5. For which observation time is it possible to have the 
narrow fields of view intersect at a fixed point on the 
lower and the higher atmosphere and what are the beam 
diameters and pointing accuracies? 

6 3 Is it possible to extend these applications to the complete 
day and night cycle by using ground beacons besides reflected 
and scattered sunlight? 
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7 SOME CONCEPTS FOR AN EMPLACED SCIENTIFIC 
STATION ON THE LUNAR SURFACE 

By 

/ 

; Russell D. Shelton / 

/ Marshall Space Flight Center j 

INTRODUCTION 

N67-24277 

The concept for an emplaced scientific station (ESS) has been 
studied at the Marshall Space Flight Center for the past year and at 
the Westinghouse Corporation under contract NAS8-20245. Tht fig- 
ures used in the following discussion were taken from the Westinghouse 
Final Report, and the concepts discussed resulted from a joint study 
effort. 

The ESS concept was based on the following general scientific 
objectives (emphasizing that observations were to be made for a year): 

1. To acquire data on the lunar interior by means of 
seismic detectors and probe measurements 

2. To monitor the lunar atmosphere 

3. To measure fluxes of nuclear and meteoroid particles 

4. To study the moon as a possible base for improved 
planetary and stellar astronomy 

The physical requirements of the ESS are such that it can be stored 

3 

in a volume of 0. 6 m , weigh 290 kg, or less, and packaged suitably 

for astronaut deployment. For seismic purposes, the ESS array would 

consist of 4 stations: a central station for communications with the 

earth as well as for experimentation and three identical satellite sta- 

o 

tions deployed at 120 -arc intervals on a circle of 8 -km radius. 

The central and satellite stations would be a standardized package 
designed to make possible a wide latitude in the selection of experi- 
ments and to support a variety of scientific instruments. 

Westinghouse group directed by Vernon B. Morris under Contract 
No. NAS8 -20245. 
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This particular philosophy seems desirable in a situation in which data 
are being acquired continuously, instrumentation is being improved, 
and new ideas are evolving. 

The packages would use local signal conditioning and passive 
thermal control to conserve power. The central station would include 
an orientable platform for a number of observations. The general 
appearance of the deployed central and satellite stations of the ESS 
related experiments are shown in Figs, 1 and 2. 

• The design of the ESS was based on the considerable interest in 
the measurement of physical quantities which vary in time and on 
experiments which depend on simultaneous readings from points widely 
separated on the lunar surface. Considerable emphasis was given to the 
discussions held at Woods Hole and Falmouth, Massachusetts in the 
summer of 1965 [l, 2] . The satellite stations transmit data to the 
central station where it is stored for subsequent transmittal to earth. 
Astronauts will deploy the central and satellite stations; therefore, 
their safety and work limitations must be considered. Presumably, 
the ESS will use man's capabilities to obtain effective data for appli- 
cations in future mission planning with a minimum of effort. 

Experiments yielding information on the lunar atmosphere and 
surface rank high in priority with emphasis on the experiments which 
might possibly suffer from contamination by radioactive materials or 
exhaust gases from landing rockets. An attempt will be made to study 
the lunar subsurface by active and passive seismic techniques and the 
surface structure by radiometric and probe techniques. 

There are a number of observations of the earth which can be 
made effectively from a lunar base, particularly of phenomena which 
are seasonal and which cover a large part of the surface of the earth 
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FIGURE 2 - DEPLOYMENT OF A SATELLITE STATION 


at one time. Examples of such phenomena would be ionospheric dis- 
turbances, atmospheric activity, oceanographic properties, and items 
of geodesic interest. It is possible also that a number of observations 
of interplanetary space, the sun, and the planets could be made more 
effectively from the moon. 

It is probably true that no disciplines of the physical sciences 
have been able to do more with small clues than astronomy, cosmology, 
and astrophysics. Consequently, it could be of great interest to con- 
tinue on a long-term basis the improved stellar, interstellar, galactic, 
and intergalactic observations possible from the lunar surface. 

Naturally, by the time an emplaced scientific station on the lunar 
surface is possible, there will be considerably more added to the 
background of instruments available for space measurements. It is 
a widely accepted axiom that maximum use should be made of proven 
instruments which, if man-operated, can give the largest amount of 
significant data at the lowest cost. Instrument groupings which pro- 
vide redundancy, continuity, and inter -related measurements are of 
course desirable. Naturally, the planning for any experiment involves 
some estimate of weight, size, power, and environmental control. 
Operational constraints such as pointing or handling, instrument sensi- 
tivity to RFI and environmental changes, command and control require- 
ments, etc. , must also be considered. 

RADIATION MEASUREMENTS 

There are a number of radiation experiments which can be done on 
the lunar surface. Of almost universal interest is the total radiation 
dosimeter which provides an indication of the astronauts' exposure. 
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The vibrating -reed ionization chamber is probably a good choice for 
such an instrument. The nuclear emulsion, which would be returned 
with the astronauts, still provides one of the moat reliable means of 
identifying many particles and establishing their energy and directional 
distribution. Other instruments of interest on the lunar surface would 
include: 1) a magnetometer, which is highly developed at the present 

time; 2) an electric field meter, which presently is the subject of 
much discussion; and 3) a VLF meter, which could measure some 
of the interactions of the solar wind with the moon. 

LUNAR BODY MOTIONS 

The lunar body motions can be studied by means of a three-axis 
passive seismometer, by geophone, tidal gravimeter, or UV tele- 
scope. The lunar surface can be studied from the viewpoint of radio- 
meters, borehole temperature probes, and perhaps lasers for map- 
ping and ranging. 

The seismic studies will benefit from the largest possible base- 
line array for information on direction, distance, and depth of lunar 
tremors. Although a central station and two satellite stations may appear 
to be sufficient from a theoretical viewpoint, it is felt that, from the 
viewpoint of reliability and ease of data interpretations, three satellite 
stations are needed. The stations should be roughly 8 km apart, as 
shown in Fig. 3. 




TELESCOPE APPLICATIONS 


The observations of the earth and the astronomical studies will 
probably require the pointing of multiple instruments for simultaneous 
imaging, and spectrometers will require simultaneous studies at dif- 
ferent wavelengths. The most important ingredient of the orientable 
platform planned for the ESS would be a Cassegrainian-type 12" re- 
flecting telescope. 

In assessing the value of a telescope, there are a number of con- 
siderations. It is of little use to design a large expensive telescope 
if there is not a parallel quality in optical accuracy, pointing accuracy, 
platform stability, and tracking capabilities. The specific capabilities 
of the ESS should include a lunar seismic and gravity network, a solar 
astronomy observatory with a 12" diameter telescope with a receiving 
capability from 10 microns to 1/10 micron, and a radio observatory. 

A lunar-based telescope could enjoy several advantages, such as 
a very long day for thermal equilibrium, no atmospheric absorption, 
a platform so stable that many good pictures could be taken even if 
the tracking devices failed, and a good means of tracking lunar land- 
ings and astronaut traverses. Image intensifiers for TV transmission 
and other advantages may make it possible to realize the theoretical 
capability of a small telescope on the lunar surface. One use of a 
lunar -based telescope would be to study the lunar motions, some of 
which are illustrated in Fig. 4, which shows the apparent motion of 
the earth in the lunar sky. 
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In the central scientific station the instruments listed in Table I 
should be available. The listing reflects a crude order of priority 
which, of course, can be attacked from several viewpoints such as 
the area of most experimental importance, the interrelationships of 
various measurements, the reliability and availability of instruments, 
etc. For example, during the night or day the electric field meter 
measurement should be associated closely with the solar wind measure- 
ments and the magnetometer measurements. 

GENERAL CONSIDERATIONS 

The radioisotope thermoelectric generator appears at this time 
to be the best choice for providing power to the ESS, although it is 
possible that a regenerative fuel cell plus a solar cell array will be 
a close contender. The advantages of the RTG power source are that 
it provides continuous power, is light in weight, occupies a small 
volume, is reliable, and is a source of heat. The disadvantages are 
that it presents a heat and radioactivity hazard during handling, re- 
quires special arrangements to dissipate (before deployment) the waste 
heat power, which is roughly 20 times the electrical power output, is 
of limited availability without long-range program planning, and is 
quite expensive. 

A weight, volume and power summary of the ESS is shown in 


Tables I and II. 


TABI -E I ESS WEIGHT AND VOLUME 
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1 SIMULATED WEIGHTLESSNESS UNDER NEUTRAL BUOYANCY (, 

By 

Robert D. Dean 

t 

and 

Robert P. Langan 7 1 ' 

I The Boeing Company , 

N67-24278 

The use of neutral buoyancy is a relatively new zero-gravity simulation 
technique that permits measurement at operator performance in full-scale mock-ups 
with six degrees of freedom. This includes the use of a large water-filled tank 
in which subjects and work materials are ballasted to approach weightlessness. As 
a simulation of zero gravity, it is limited by the viscosity of water and resulting 
hydrodynamic damping effects, which precludes its use when rapid, gross body 
movements are required of the subject, or when very large objects must be moved. 

It is also limited by the density of water ( 800 times greater than air) , which pre- 
cludes gross changes in depth when compressible materials are used. A further 
limitation is that during deflection, the body entrains the water for a period of 
time after the initial force of movement disappears. This is equivalent to an 
apparent increase in mass and results in further deflection. 

In a true zero gravity environment, there are no damping effects, no 
density of the surrounding medium , and no apparent increase in mass. The sig- 
nificance of these limitations will be fully understood only when data obtained in 
neutrally buoyant environments can be correlated with data obtained in free space. 
These comparisons will be possible when more complete performance data are 
available from future manned space flights. In the meantime, the simulation 
technique is good so long as the limitations are recogniz 

, i 

In spite of the inherent limitations of neutral buoyancy, experience indicates 
that it can provide an inexpensive and reliable means of evaluating designs for 
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space-vehicle hatches, compartments, airlocks, restraints and transfer devices. 
Comparisons between performance under one-g and neutrally buoyant conditions 
can provide designers with an understanding of the qualitative effects of zero 
gravity; some tasks such as hand tool manipulation become more difficult in 
the neutrally buoyant environment, others such as hatch traverse become easier. 
Neutral buoyancy can also provide an invaluable training aid, in that it permits 
practice on performance in an environment that approximates true zero gravity. 

Boeing work with neutral buoyancy was begun in September, 1963. Five 
major studies were conducted from September, 1963, to August, 1965. Early 
studies were conducted in an underwater facility located in a small lake outside 
Seattle. Studies are now conducted in a 25-foot diameter tank specifically de- 
signed for neutral buoyancy work. 

The first study was completed in November, 1963. This study evaluated 
the effects of neutral buoyancy on the time required to traverse 20- and 28-inch 
circular hatches. Subjects wore wet suits and scuba equipment ballasted for 
neutral buoyancy. Average hatch traverse time was approximately half that in 
the one-g environment. 

The second study, in July, 1964, was designed to evaluate the effects of 
neutral buoyancy on man's ability to manipulate hand tools (wrenches and screw 
drivers). One task, fastener removal, was used throughout. Again, subjects 
wore wet suits and scuba equipment modified to achieve neutral buoyancy (Fig. 1) . 
There was a 50 percent increase in average time required to remove fasteners 
compared to one-g control data. 

The third study, completed in February, 1965, was concerned with extra- 
vehicular assembly. The subject wore a full pressure suit ( 3. 5 psig) and life 
support backpack. The suit and equipment were neutrally buoyant ( Fig. 2) . 
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FIGURE i. SUBJECT PERFORMS FASTENER REMOVAL 


Lighting was arranged to provide a single source of high-intensity illumination 
to simulate sunlight. The task was to assemble a 12-foot length of triangular 
truss structure. The results suggested that assembly operations could be per- 
formed with greater ease while wearing a 3. 5 psig pressure suit under water 
than under shirt-sleeve conditions in a one-g environment. A significant result 
of this study demonstrated the need for self-locking connecting devices to prevent 
them drifting away. A need was also indicated for diffuse reflecting surfaces on 
the gloves and for semi-rigid non-snagging tethers. 



FIGURE 2. EXTRAVEHICULAR ASSEMBLY 

The fourth study, completed in April, 1965, was to assess a variety of 
restraint concepts. These included hand rails in traversing a 30-inch circular 
hatch (Fig. 3) , toe rails (Fig. 4) , lap belts (Fig. 5) and Velcro tape as 
restraints for performing console operations, and Velcro tape as an aid in tra- 
versing the length of a compartment. The subjects wore coveralls designed for 

* 

neutral buoyancy simulating shift-sleeve clothing that may be worn in future space 
vehicles and a standard scuba regulator connected to a specially designed scuba 
chest pack. The results indicated that although the hand rails had no effect on 
the time required to traverse the hatch, they did provide greater control over 
body movements in preventing inappropriate actions such as stepping on console 
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, FIGURE 3. USE OF HAND RAIL 
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id controls. Toe rails proved to be very helpful in positioning the body. Both lap 

4'i 

belt and Velcro tape were unacceptable as restraints when performing console 

\s 

\ operations. The lap belt made it difficult for subjects to change position, and 

the Velcro tape gave insufficient support. Velcro tape increased average com- 
7-j partment traverse time by a factor of five. 

i I 

The fifth study, completed in August, 1965, had two objectives. The 

", r 

; first was to assess a prototype airlock design to determine whether it was large 

enough to permit the astronaut to turn around inside the airlock ( Fig. 6) , to 
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FIGURE 4. USE OF TOE RAIL 
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FIGURE 6. ASSESSING PROTOTYPE AIRLOCK DESIGN 


permit exit and entry through an extravehicular hatch ( Fig. 7) , and to permit 
donning and doffing of a simulated Apollo suit. The second objective was to 
evaluate the effects of weightlessness on ergometric performance in terms of 
push, pull and torque (Fig. 8). The study results indicated that the airlock 
design permitted performance of all three tasks. The ergometric data indicated 
that weightlessness degraded the ability to push and to torque, but facilitated 
the ability to pull. Ability to pull was apparently improved as a result of greater 
freedom to position the legs. Comparisons between data obtained under different 
suit pressures indicated that increasing suit pressure irom 0 to 3. 5 psig re- 
sulted in poorer performance. 
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FIGURE 8. SUBJECT PERFORMS TORQUE TEST 
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The neutral buoyancy experimental approach used at Boeing has been one 
of repeated measures on a small number of subjects in a wide variety of experi- 
mental conditions. Despite the obvious difficulty of generalizing to a larger 
population, this approach has four distinct advantages. First, a large sample 
of subjects , trained to work under water in a full pressure suit, is not available. 
Second, equipment of sufficient reliability and quantity to permit rigid experi- 
mental research is still under development. For example, a minimum of seven 
pressure suits, modified for underwater operations, would be required to con- 
duct a study involving five subjects ranging in anthropometry from 5th to 95th 
percentile; Boeing presently has three such suits (Fig. 9) . A third advantage 
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FIGURE 9. VIEW OF PRESSURE SUIT 
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is that critical problem areas can be identified without incurring the high cost 
of obtaining normative data on large subject samples. The fourth and most 
significant advantage of the present experimental approach is the use of neutral 
buoyancy as a design tool. When used as a design tool, rather than as a re- 
search device, neutral buoyancy becomes a part of an iterative process in which 
a design concept (for example, a restraint system) is initially formalized and 
then tested (Fig. 10). The initial concept is rejected, or modified as a result 




FIGURE 10. CHECKOUT OF DESIGN CONCEPT FOR A RESTRAINT SYSTEM 

of the test data, and tested again. This process of successive design modifica- 
tion and testing is continued until an acceptable design is developed. When the 
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neutral buoyancy simulator is used in this manner, there is rarely sufficient 
time between successive iterations to obtain data on more than a few subjects. 
Experience has shown that intelligent redesign can typically be undertaken after 
only a few preliminary observations in the neutrally buoyant environment. When 
this occurs, as it usually does, there is little justification for continuing an ex- 
tensive evaluation of a design concept that is already in the process of modifi- 
cation. 


Perhaps, in summation, the advantages of time and cost should be 
mentioned. Compared with Keplerian trajectory flights, in which weightlessness j_ 

is achieved only for some seconds, the time advantage is tremendous. Testing j- 

j' 

in the neutral buoyancy facility can be conducted during a normal eight-hour \' 

day or on a multiple shift basis. Since aircraft and flight crews are obviously 

not necessary, there is a very respectable cost reduction in equipment and 

training, as well as a significant reduction of h.zard. Further, materials and 

designs to be tested are limited in size only by the volume of die facility used. i 


The potential as well as the immediate value of this technique as an aid 
to engineering analysis and design has been proven. It is relatively simple and 
economical to use in comparison with other methods, but the returns can be 
great. 
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. DEVELOPMENT OF A DEVICE FOE PHYSICAL 
CONDITIONING DURING WEIGHTLESSNESS 

By 

( Gary Chase, Caswell Grave, John Goode, Gary Graham ) . 
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ABSTRACT 


From analytical studies it was deduced that an exercise device, which 
would impart pulsatile head-to-foot acceleration to human subjects to maintain 
their gravity tolerance, would be useful in space vehicles. The criteria re- 
quired that the subjects receive plus and minus 2G longitudinal acceleration, 
each acceleration to last three milliseconds. The completed device resembled 
a double trampoline, the subject vibrating on a cart between the impacting 
surfaces. A concept for a light-weight model for eventual application in a space 
station was also developed. 
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INTRODUCTION 


The maintenance of the functional Integrity and physical well-being of 
astronauts during exposure to a weightless environment has generated 
considerable Interest In the Aerospace community. It is generally 
agreed that some physical activity (exercise) will be desirable on short 
missions fluid necessary to prevent physical debilitations on longer, 
perhaps, 30 -day missions. 

Development of apparatus to achieve the goal of maintenance of well- 
being requires consideration of the undesirable physiological changes 
which would contribute to loss of physical capabilities of the astronaut. 
Experience has shown that prolonged lack of exercise leads to reduction 
of physic flil strength and the ability to do nuous physical tasks. 

Among other things, research has revealed a decreased ability to trans- 
port oxygen in the blood stream to the muscles, as well as calcium loss 
from the bones. An equally significant finding is that humans In simula- 
tions of weightlessness tend to lose their ability to stand erect In the 
earth's gravitational field without fainting. This ability Is termed 
orthostatic tolerance. Apparatus development should thus provide methods 
to flex the muscles, stress the skeleton and provide a force or treatment 
to substitute for gravity la the normal one-g environment. 

So work has clearly shown that there is an efficient substitute for voluntary 
exercise In muscular and skeletal maintenance although many techniques 
such as massage, electrical stimulation and passive exercise have been 
employed. \An acceptable device for maintaining well-being would thus 
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Incorporate the ability to engage in voluntary exercise. Many approaches 
have been used, and such devices as bungee exercisers. Isometric devices, dyna- 
mometers, roving machines and bicycle ergometers, all of which can be used In 
space vehicles, have been developed. A shortcoming of this group Is that 
they do not prevent the development of some degree of deconditioning. 


The development described here was designed for testing In an experimental 
simulation - bed rest - of the zero-g environment. It provides conditioning 
for the muscular, skeletal and cardiovascular systems . j This Is done by 
Imparting cyclic, pulsatile accelerations along the long axis of the body. 

A runner-mounted carrier supports the subject in a face-up, horizontal position 
and Is free to move back ai._ forth between two vertical, opposed trampolines 
(Fig. l). The man oscillates by bouncing between the trampolines over a 
distance necessary to develop the desired acceleration-time profile. Such 
treatment provides a substitute for the conditioning normally afforded humans 
in their daily, one-g living In which gravitationally induced hydrostatic forces 
in the circulatory system are thought to play a significant role. It was 
originally hypothesized (Grave et al, 1964^) that oscillating linear accelera- 
tion would maintain normal cardiovascular reflexes if the duration of the 
acceleration was sufficient to cause fluid shifts In the circulatory system.* 


*Tbe hypothesis was based on observations by Whodon et al (1949) 2 ^ that 
oscillation of bed rest subjects on a Sanders bed largely prevented the 
development of orthostatic intolerance during simple bed rest, as reported by 
Deitrick et al (l9i<6)(3). Other pertinent observations by Graveline ( 1962 ) 1 *) 
showed that Intermittent venous occlusion by cuffs on the leg and arm prevented 
a similar development of orthostatic hypotension during water Immersion. It 
was further suggested that fluid shifts which accompanied operation of the 
oscillating had and occlusive cuffs could be produced by pulsatile accelerations 
(vibration) with much shorter time factors thsn were employed with the oscil- 
lating bod (105 seconds psr cycle) or .cuffs (120 seconds per cycle). Studies /ft . 
by 8 tapp (l96l)(*/, Xbmhauser (is 6 l)w), Thompson ( 1962 )'?), and Oaggs (1950)1®) 
showed that the threshold for fluid shifts probably lies in the range from 
70 to 3$0 milliseconds. Accordingly ws dasignort our device to furnish tbs /■ 
subject an Impact +g_ for a duration of at least 300 milliseconds, with a 
peak fores of +2 g. 
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We have conducted laboratory studies In the Biotechnology organization of 
The Boeing Company over the past two years using long tern bed rest as the 
method of simulating the effects of reduced gravity on human systems. During 
these tests, subjects were given dally exposures to passive and/or active 
peak rebound accelerations of 2 g's at each end of Hie bed travel. Only 
the active mode provides a significant amount of Isotonic exercise. This 
exercise provided test subjects with the physical conditioning required to 
complete tests of orthostatic tolerance without major 111 effect. A detailed 
report of the procedures and results of these tests will be presented 
elsewhere. 


DETAILS OF MACHINE CONSTRUCTION 

The carrier Is mounted on precision ball bushings running on hardened steel 
shafts. The shafts are supported between canvas tarpaulins suspended by 
shock cord on frames at either end of the shafts, resulting In a double 
trampoline device (Figure l). The platform Is padded for comfort and fitted 
with adjustable restraints to allow for a variety of positions and sizes 
of test subjects. The bed can be propelled either by a standing operator 
(passive exercise) or by the legs of the test subject (active exercise). 

The quantity of exercise can be adjusted by changes to the mass of the cart 
and regulation the duration of the ride. The accelerometer provides 
a continuous readout on a Brush recorder so that by monitoring the peak 
acceleration, tbs efforts of the operator or subject can be adjusted to 
produce the desired acceleration profile. From time to time It was also 
necessary to replace or adjust the tension of the shock cord supporting 
the tarpaulins since this affected the time profile of the rebound 
accelerations. 
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HARDWARE ADAPTATIONS OF METHOD TO SPACE VEHICLES 
During laboratory tests periodic engineering design studies were being 
conducted. Some ideas for the application of the "trampoline' 1 device 
are shown In Figure 2. 

It is believed that the configuration labelled Concept A is probably the 
best and can easily fit the dimensions of most space laboratories. 

In order to provide a complete physical conditioning system in one location 
it was deemed desirable to locate other exercise, and test equipment on or 
around the trampoline . 

CONCLUSIONS 

Only by means of orbital laboratory studies utilizing the trampoline or some 
similar system of physiological support will the true need for a conditioning 
program become clearly established. Quite possibly the magnitude of the 
deconditioning problem will Jpe proven Inconsequential and further consideration 
will be totally unnecessary. However, nobody can m a te a certain prediction 
with the information at hand, and we must anticipate the worst possibilities 
until the critical data became available. Manned space flight and the 
emphasis here rides on the man component - obviously imposes a different 
and more demanding set of objectives on the total mission, as contrasted to 
unmanned flight. We must plan for a safe round trip of our astronauts and 
recognize that not only are they an Important part of the scientific package, 
but their viability also carries an acknowledged and very large political 
impact. These considerations suggest analysis and test of the trampoline 
device along with other methods for prevention of deconditioning of the 
astronaut.! 
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FIGURE 1. EXPERIMENTAL VERIFICATION 



PROPOSED APPLICATIONS 
FIGURE 2 
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HYPERSONIC NONEQUILIBRIUM FLOW AT RE-ENTRY OF 
SPACE VEHICLES AND ITS THERMODYNAMIC RELATIONS ^' 2 

By 

3 ' 

. Rudolf Hermann ‘ i 1 ■ 


I University of Alabama Research Institute 


INTRODUCTION 


N67-24280 


With the availability of the Saturn/Apollo vehicle, it will be possible to perform 
many different types of space missions with a broad range of applications. In many cases, 
the space vehicle has to return safely to the Earth's surface after completion of its mission. 
Even though re-entry and soft landing of space capsules, such as Mercury and Gemini, 
have been technically mastered, analytical calculations of flow fields at hypersonic 
velocities, as they are encountered during re-entry from outer space, are still in the 
beginning phase and include many uncertainties. Therefore, all present investigations 
make use of simplifying assumptions. 

In this paper, the following two restrictions have been employed: 

First, with respect to the body shape, the calculations presented are carried out 
for a circular cylinder, the prototype of a, blunt body. The next logical step is the 
calculation for a sphere, which is a first approximation for the front portion of the Apollo 
Command Module during re-entry. This task has been attacked by various laboratories 
including our Institute. 

Second, with respect to the flow field, this paper deals with the inviscid flow region 
up to the equator of the body in such a range of temperatures and pressures where onl y 
oxygen dissociation and, of course, molecular vibration occur. Later calculations should 
include nitrogen dissociation and formation of nitric oxide. After the inviscid flow region 
is known, the calculation of the boundary layer and the wake can be started. 


^ The preparation of this paper was supported by the National Aeronautics and Space 
Administration under research grant Ns G -381 . 
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NONEQUILIBRIUM HYPERSONIC FLOW FIELDS AND RE-ENTRY CORRIDORS 

Flow Fields and Flight Regions During Re-Entry 

The configuration of the Apollo vehicle with the schematic flow field is shown in 
Figure 1 (Ref. 1). It is a typical blunt body which, in high speed flight, generates a 
detached bow shock. The front region between the shock wave and the outer edge of 
the boundary layer, called the inviscid shock layer, has a subsonic portion in the vicinity 
of the stagnation point; farther downstream the flow is supersonic. Besides the subsonic- 
supersonic inviscid shock layer , there is the boundary layer, and behind the body, a 
viscous separated region. The latter two regions interact in the viscous mixing region 
which poses almost unsurmountable difficulties for a theoretical analysis. Unfortunately, 
because of the location of antennas, this is also an important region as far as electro- 
magnetic wave propagation is concerned. 

For axisymmetric flow at zero angle of attack, flow fields are symmetric with I 

respect to the body axis; and the inviscid portions can be calculated by presently avail- 
able methods even though they require considerable computational effort. In the case ’ 

i 

of Apollo re-entry flight, the capsule itself is still an axisymmetric blunt body, but 
typically it descends at some angle of attack (to c. maximum of 33°) which is varied | 

for control purposes. The flow field is now unsymmetric, thus, adding a major complication 
! to the problem, and analytical methods are not readily available. 

Before an analysis of the boundary layer or the wake can be made, the inviscid 
; flow field must be known. The initial effort must, therefore, be directed toward the 

determination of the inviscid flow field. Since air is a rather complicated mixture of 
gases, especially when dissociation and ionization must be considered, a truly exact 
! representation is not possible at the present time. In order to properly select a model , 

it is advantageous to study the conditions that are encountered along a typical re-entry 
trajectory, usually plotted in a velocity-altitude diagram. 

In Figure 2, the cross-hatched region indicates the re-entry corridors of manned 
space capsules from circular orbit around Earth with about 7,9 km/sec and for lunar return 
with about 11.3 km/sec. Also included are equilibrium conditions behind a normal shock; 
thus the temperature, the pressure, and the density as they occur in the stagnation point 

I 

i 
i 

i 
i 
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region of a blunt body re-entering the Earth's atmosphere can be read from the graph. 
Thermodynamic data for the graph were taken from References 2 and 3. 

It is interesting to observe that a major portion of the space vehicle trajectory 
is approximately parallel to a line p $ = const. The temperature lines at low velocities 
are practically vertical, that is, the temperature depends only on the square of the 
velocity. For higher velocities, the temperature depends on both velocity and altitude. 
The lower density of high altitudes has the effect of increasing the degree of dissociation 
which in turn results in a reduced temperature since the dissociation process absorbes 
energy. Nevertheless, the stagnation temperature reaches very large values, up to 
] 1,000 °K in the case of lunar return. 

Nonequilibrium Flow Features 

In most hypersonic applications, the flow field in the shock layer will generally 
be in some nonequilibrium state, because there is insufficient time to obtain thermody- 
namic equilibrium. Depending on free stream velocity, density, temperature, atmospheric 
composition, and the absolute size of the body, two limiting cases exist; either the flow 
may be almost frozen, or the flow may reach very nearly equilibrium. 

A gas at rest is, by definition, in thermodynamic equilibrium, if a particular volume 
of the gas has sufficient, -or better, infinite time to bring all its internal modes of energy 
in equilibrium with the translational energy of the molecular motion. For our considera- 
tion those modes are molecular vibration, dissociation, electronic excitation, and 
ionization. Considering flow processes of a gas, it is obvious that equilibrium flow is 
only one limiting case, namely, when the changes of the state of the gas flowing along 
a streamline are so slow that at any point equilibrium is obtained, or stated more exactly, 
equilibrium is very closely approached. 

At hypersonic velocities, the time available is, in general, too short for the gas 
particles which are undergoing rapid density, temperature, and composition changes to 
reach thermodynamic equilibrium. Hence, in general, there will be nonequilibrium 
flow. The degree of molecular vibration, the degree of dissociation (chemical compo- 
sition), and the degree of ionization will still change from point to point along the 
streamline but will not reach thermodynamic equilibrium at any point. 


\\ 
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The second limiting cose occurs when the gas moves so fast that the internal 
energy modes have no time to follow the changing density and temperature with the 
result that the vibrational energy, the energy in dissociation, and the energy in 
ionization stay very nearly constant. This flow is called frozen; Jhe gas might be 
vibrationally frozen, and/or chemically frozen (frozen dissociation or no change in 
degree of dissociation), and/or the gas has frozen ionization. 

The departure from equilibrium of a flowing gas depends upon the ratio of a 
characteristic translational time needed by the macroscopic motion of the particles to 
cover a typical distance, relative to the magnitude of the respective relaxation time 
which is different for each mode of excitation. Since the characteristic translational 
time is dependent on the body size, and since relaxation times are dependent on the 
temperature, density, and to some extent on local composition, it is obvious that 
nonequilibrium flow fields are generally not similar for geometrically similar bodies 
even at completely equal free stream conditions; this is a very important fact. 


THERMODYNAMIC RELATIONS INCLUDING REAL GAS EFFECTS 

Simplified Air Model, Degree of Dissociation 

At the present time it is not possible to consider simultaneously all possible 
individual reactions which can occur between all components in air. Therefore, we 
must concentrate our efforts on the study of those components and reactions that are 
of greatest engineering concern today. 

The simplified air model used in the calculations presented consists of oxygen 
and nitrogen only. The chosen composition in undissociated state is 21 per cent Oj 
and 79 per cent Nj, by volume. Furthermore, we will restrict the range of application 
expressed in the appropriate ranges of temperature and pressure such that only oxygen 
may dissociate. 

The mole fraction of oxygen molecules in our air model, before dissociation, 
is expressed as 

"/v» 

b - 0.21 (1) 


02 


n 02 + n N2 


where nj denote the number of particles of the ith qwcies per unit volume of the gas. 
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In general, at a given pressure and temperature, only a certain fraction of the 
molecules are dissociated into atoms. The degree of oxygen dissociation, a , can be 
defined as the ratio of the number of dissociated oxygen atoms per unit volume to the 
initial number of undissociated oxygen atoms per unit volume which are, of course, in 
molecular state. 


* 


Thermal Equation of State, Internal Energy, and Enthalpy Equation 

Assuming that the individual species obey the perfect gas law, the total pressure 
of the mixture is the sum of its partial pressures, or 


p = ?p. = ?p. R.T 


( 2 ) 


? ’ i y ’ I i 

Combining expressions for p, , Rj , a , and b, we finally obtain 

p = (1 + ba) p R T - ZpRT (3) 

where Rj is the specific gas constant of each species, R is the gas constant of the 
undissociated mixture, and Z * (1 + ba) is called the compressibility factor. 

The internal energy of the gas mixture, E , is the weighted sum of the internal 
energies of each species, eg , due to the various modes of excitation (translation, 
rotation, vibration) and of the dissociation energy, D , of the dissociating component, 

E M'l*l + T 0 O2 < 4 > 

whore f. denotes the mole fractions for the dissociated mixture. 

i 

The enthalpy of the gas mixture is defined in the same way as the enthalpy of 
a gas with only one component 


h * E + ^ • 


(5) 


After determining the mole fractions, and obtaining the internal energy of the mixture 
from statistical thermodynamics, the enthalpy equation finally becomes ^ 


baDJ 


Q2 + f baT + jT + (I - a) b 


0 


02 
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r- + o 




®N2 

- b) 7^T. • 


Here 0j is the characteristic temperature of vibration (2270 °K for 3390 °IC for Nj), 
and Dqj l» *bo characteristic temperature of oxygen dissociation (59366 °K). 
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Rate Equation, Dissociation Rate Constant, and Equilibrium Constant 

In general, a large number of collisions among the particles is required to bring 
molecuior vibration, dissociation, and higher degrees of excitr ‘ion in equilibrium with 
the local translational temperature. This means that a finite amount of time (relaxation 
time) is needed for the gas properties to approach thermodynamic equilibrium. The 
departure from equilibrium of a flowing gas is characterized by the magnitude of this 
relaxation time relative to some translational time (residence time) needed by the particles 
to move over a characteristic distance on a body within a nozzle. Since only oxygen 
dissociation is considered, only one chemical rate equation is needed. Molecular vibrations 
are assumed to be in equilibrium. 

In ?he process of dissociation, an oxygen molecule must collide with another particle 
having enough energy to break up the oxygen molecule. The recombination process is the 
inverse of the dissociation. In order to occur, it is necessary that two atoms collide at the 
same time with a third body, which is able to carry away the energy that the two separate 
atoms must release to form a stable diatomic molecule (three-body collision). 

In a flow field with a local velocity vector q and a varioble degree of dissociation, 
characterized by grod a , the following relation con be derived, called the rate equation: 

q • grad a * F (7) 


where 


CVZK, 


(1 -<<) 
TCp 


2 * 

a_ 

” K 

c.. 


The expression for F is called the source function for oxygen atoms. Among others, it 

contains the number of oxygen atoms per unit mass of the gas, C , furthermore the 

dissociation rote constant, k^ , which depends strongly on the type of the colliding body 

and some other gas properties, and finally the concentration equilibrium constant, K , 

which is a function of temperature cnly. In order to evaluate the rate equation (12), we 

still need expressions for k . and K . 

o c 

The dissociation rate constant, k rf , cannot be predicted by ordinary kinetic gas 
theory. Wray (Ref. 4) has given a comprehensive review of recent work carried out in 
that field, and hts values will be used In the present investigation. In order to account 
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for the influence of the catalytic species, we are using a dissociation rate constant, which 
is averaged with the respect to population. This finally results in the expression 
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It can be seen that the dissociation rate "constant" is actually a function of temperature, T , 
and composition, a . Fig. 3 shows the range of this equation for a between 0 and 1 , 
together with the values of other authors for various colliding bodies, k^ has the dimension 
particle ^ sec ^1 . 

Equation (8) was obtained through el imination of the recombination rate constant, k^ , 
by introducing the concentration equilibrium constant, K . This is a function of temperature 

C . . -3 

only and results from statistical thermodynamics like k^ . has the dimension [particle m ] 
and its final expression for the numerical calculation is: 
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Only a few important relations needed for the following flow calculations have been 
given in this section. More details can be found in Reference 5. 


NONEQUILIBRIUM FLOW AROUND BLUNT BODIES WITHOUT AND WITH 
DISSOCIATION IN THE FREE STREAM 

Direct and Inverse Methods for Flow Calculation 
The importance of blunted shapes in the hypersonic flight of a vehicle re-entering 
the atmosphere from outer space has led to many investigations of hypersonic flows with 
detached shock waves, using both inverse and direct methods. At the inveise method, a 
certain shock shape is assumed or prescribed, and then the governing equations are integrated 
through the flow field, and the associated body shape follows from the calculation. The 
author and his associates at the University of Alabama Research Institute in Huntsville have 
applied a direct method, where the body shape is given and the flow field, including the 
shock shape, is a result of the calculation. This method uses "integral relations", first 
described by Dorodnitsyn (Ref. 6) In 1959 for the solution of two-dimensional boundary value 
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problems, later applied by Belotserkovskii (Ref. 7) to the calculation of supersonic flow 
of a perfect gas past a circular cylinder. 

Up to now a limited number of investigators obtained results for real gas flows in 
connection with the integral method. Recently Hermann and Thoenes (Ref. 8) have given 
new extended and improved results for nonequilibrium hypersonic flow of air past a circular 
cylinder, which are presented in this section. They have also treated, for the first time 
to the author's knowledge, the case of dissociation in the free stream. Such conditions are 
encountered either in the atmosphere at high altitude (above 90 km), or in the nozzle of 
high temperature hypersonic flow facilities designed to simulate re-entry flight conditions, 
where appreciable frozen dissociation may occur. 

Basic Equations for Inviscid Flow Past a Circular Cylinder 

We are restricting our analysis to inviscid flow and, in addition, we are neglecting 
heat conduction and radiation. Figure 4 shows the cylinder with the selected polar 
coordinate system . The basic equations of motion for steady adiabatic flow, after trans- 
formation from the vectorial form to polar coordinates, yield the following first three partial 
differential equations. The fourth equation, expressing the conservation of oxygen atoms, 
had to be added since the gas under consideration is reacting. F is the source function of 
oxygen atoms, Equation (8). 

Conservation of mass: 




k 4>u ) + Jr (pv,) ’ 0 

(ID 

0-momentum: 


8u 9u . 1 9p _ A 

u 9S + v, 5r + uv + F55"° 

(12) 

r-momentum: 

8 

8? 

(puv) + [<P + P» 2 )rl - (p + P- 2 ) " 0 

(13) 

Rate equation: 

• 


u f|* vr *_ Fr ,0. 

(14) 


In order to solve these equations for the flow field variables, they must be supplemented 
by the equation for the conservation of energy, containing an expression for the enthalpy. 
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Eq. (6), and by the thermal equation of state, Eq. (3) . Note that the first three equations 
(conservation of mass and of the two momenta) are independent of the particular gas model, 
while the last three equations (rate, energy, equation of state) do depend strongly on the 
selected gas model . The conservation of mass and the r-momentum equation are given in 
the so-called divergence form, which : s required for the application of the integral method. 

Boundary Conditions 

The condition for flow tangency on the body surface is v^ = 0 . The conditions 

behind the shock are obtained from the conservation of mass, momentum and energy across 

the shock. For the present nonequilibrium flow calculation, it will be assumed that the 

chemical composition of the air does not change across the shock, i.e. , a. = a • 

I s 

Various relations between the velocity components in front of and behind the shock which 
will be needed later are obtained from Figure 4. From the same figure, one can also 
derive the following relation between the local dimensionless shock wave distance from the 
body, e , and the shock wave angle, a , both being functions of 0 : 

gg- = - (I + e) cot h * 9) . (15) 
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Application of the One-Strip Integral Method for the Cylinder 
Due to the boundary condition v^ = 0, it can be seen that the 0-momentum equation 
and the rate equation may be used in their exact forms. Hence, only two equations of the 
set, namely the conservation of mass and the r-momentum equation, have to be approximated 
by assuming a linear variation of certain integrands across the shock layer. Thus, an 
integration of the equations in the direction of the radial coordinate can be performed, 
which in turn results in two ordinary differential equations with the tangential coordinate 
as independent variable. 

The result is the following complete set of governing equations: 


Continuity: 
du. 


du dp 


vipi uu up | . r "I 

% w + % sr * p , w + ", W- " -r 1 L (p b “b ' p > * p.vJ 


( 14 ) 


0-momentum: 


p b“bsr 'ar 0 


( 17 ) 
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The five ordinary differential equations (15, 16, 17, 18, 19) and the two algebraic 
equations (20,21) constitute a system of seven equations for the seven unknowns, which 
are the five parameters along the wall of the cylinder % / Pjj / Pj, / , a^ , and the two 

parameters for the shock shape, e and <x . Special subroutines during the process of the 
solution are necessary to calculate those variables in the above equations that are functions 
of some unknowns (mainly a^ and T^) and of free stream conditions. 

Before the numerical integration of the system of equations for the flow around the 
circular cylinder can be started, the stagnation point parameters, serving as initial values, 
must be determined. For this purpose, the governing equations are specialized for the 
stagnation streamline, where 0 = 0 and u L = u = 0 . 

D S 

Henceforth, the further calculation of the flow parameters consists of two principal 
sections; first, the calculation of the flow along the stagnation streamline including the 
stagnation point itself; second, the calculation of the flow variables along the cylinder 
surface and behind the shock. Both parts still require considerable analytical effort in 
getting the final equations for the computer program. This procedure is described in more 
detail in Reference (8), 
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In the first section, the resulting system of six equations can be solved for any 
particular chosen value for the stagnation point shock detachment distance e 0 . To find 
the solution, a stagnation point pressure p^ must be assumed as a first trial, then the 
system is integrated and iterated on p^ , until the assumed value p^ agrees with that 
resulting from the integration. 

After the stagnation point conditions have been obtained, the integration around 
the body is started with an unknown, but assumed shock detachment distance e 0 . The 
equations for the variables along the cylinder surface are strongly dependent on the 
selected value c 0 , and smooth transitions of all variables from the subsonic to the 
supersonic flow regime are only obtained for a correctly chosen e 0 . Depending on the 
case, between ten and fifteen iterations were needed to determine e 0 to four significant 
figures. 

A fixed-step Runge-Kutta technique of fourth-order accuracy was used for both 
the integration along the stagnation streamline and for the integration along the body 
surface. One thousand steps were chosen for the stagnation streamline from the shock 
to the stagnation point. In the second part, the step size was fixed at 0.002 radians, 
and approximately 800 steps were used for the integration from the stagnation point to 
the equator (0 = ir/2). All calculations were performed on a UNIVAC 1107 high speed 
digital computer, which is located in the University of Alabama Res' „ ch Institute. 

Discussion of the Results for the Cylinder 

For chemical nonequilibrium flow, the degree of dissociation and the temperature 
are particularly interesting parameters which are shown in Figs. 5 and 6 along the stagna- 
tion streamline between shock and body. The free stream velocity is 4300 m/sec at an 
altitude of 30 km, corresponding to Mj *14.2. It can be seen that, depending on the size 
of the body, quite different regimes of nonequilibrium flow are encountered. For relatively 
large bodies, a characteristic flow time is large compared to the local chemical relaxation 
time, which causes tf|e flow to reach the state of thermodynamic equilibrium close behind 
the shock. During this process, the degree of oxygen dissociation increases up to approxi- 
mately 90%, while the temperature decreases correspondingly within a small region behind 
the shock. On. the other hand, for small bodies, a characteristic flow time Is very short if 
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compared to the relaxation time. In this case, the flow remains essentially frozen and 
equilibrates only near the stagnation point, where the velocity approaches zero and the 
local residence time of a flow particle again is large compared to the relaxation time. 

It was found that the flow always reaches thermodynamic equilibrium at the stagnation point. 
Note that with Tj = 225 °K, the temperature scale goes from 4500 °K to 7650 3 K. For 
small bodies (nose radius< 0.1 cm) the nonequilibrium temperature in the main part of the 
shock layer near the stagnation streamline is between 2000 °K and 3000 °K higher than 
the equilibrium temperature which exists for a sufficiently large body (nose radius = 10 cm) , 

In both figures, the present results from the integral method are compared with data 
calculated by Conti (Ref. 10), who used an inverse method, also a basically different 
approach. The air model was the same as ours, but the reaction rate constants are some- 
what different. The comparison is very gratifying. 

Calculated shock shapes up to the equator of the cylinder (body radius 10 cm) are 
shown in Fig. 7 for Mach numbers 3.0 to 14.2, without and with dissociation in the free 
stream. With increasing free stream Mach number, the bow shock moves much closer to the 
body. This trend is already well known from perfect gas calculations and is retained in 
chemical nonequilibrium flow. The figure also indicates the location of the sonic point 
on the body surface, which moves toward the stagnation point with increasing Mach number. 
The shock shape deviates considerably from a concentric circle, even where the velocity 
in the shock layer is still subsonic. 

The shock detachment distance at the stagnation streamline as function of the free 
stream Mach number is presented in Fig. 8, where our nonequilibrium flow results with and 
without free stream dissociation (a^ = 0.5 ; <X| = 0) are compared with perfect gas results 
(Y = 1 .4) from References (11) and (7). Figures 7 and 8 indicate clearly that dissociation 
of the free stream, keeping all other free stream parameters unchanged causes the bow shock 
to move away from the body. One reason for this effect is that, for a dissociated free stream, 
the density behind the shock is lower than for corresponding conditions without free stream 
dissociation. The effect is seed to increase with decreasing free stream Mach number. 

From Fig. 8 it is observed that the presen* calculations for an undissociated free stream 
yield a stagnation shock detachment distance which is much smaller than the values obtained 
from perfect gas. calculations. Responsible for this effect, in the range of Mach numbers 
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approximately between 4 and 7, is the fact that our calculations include the energy of 
molecular vibration/ which is assumed to be in equilibrium. For Mach number 6.6/ the 
present calculation/ which includes dissociation/ does not show any degree of dissociation 
in the flow field (Fig. 10). Beyond Mach number 7, the influence of the dissociation on 
the shock detachment distance becomes more pronounced. At free stream Mach number 3/ 
where the temperature caused by the bow shock does not yet cause appreciable molecular 
vibration in the shock layer/ the present calculation furnishes a value which is very close 
to the known perfect gas results. Finally, in the higher Mach number range/ at M = 14.2 
our results agree closely with those of Conti (Ref. 10)/ which were obtained by an entirely 
different approach. 

The velocity distribution along the surface of the cylinder as evident from Fig. 9 is 
almost linear up to the sonic point. A distinct effect of free stream dissociation can be 
seen. 


The degree of oxygen dissociation along the noncatalytic body surface is presented 
in Fig. 10. For zero free stream dissociation/ it is seen that at M| = 6.6 no dissociation 
occurs along the body. The degree of dissociation increases strongly with increasing Mach 
number and reaches about 0.95 at Mach number 14.2. Especially for higher Mach numbers/ 
it is seen that the recombination process dominates in the subsonic regime as the flow 
expands around the body. The degree of dissociation decreases slowly until/ in the supersonic 
regime, the local residence time of a particle becomes so small compared to the relaxation 
time that the flow freezes. 

Finally, Fig. 11 shows the temperature distribution along the surface of the cylinder 
for selected free stream Mach. numbers with and without dissociation of the free stream. As 
expected, free stream dissociation has a strong effect on the temperature. Note that for 
Mj = 6.6 the presence of free stream dissociation raises the stagnation temperature about 
1000 °K. It should be mentioned tho! the pressure distribution, not shown here, along the 
body surface is similar to the "Modified Newtonian" distribution, and that free stream 
dissociation has practically no effect on the pressure. 


Comparison with Flow Around A Sphere 

From an engineering point of view, the sphere or th4 spherically copped cone has a 

*) 

greater interest than the cylinder. However, only a few data for the sphere are available 
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for nor. -jcju! librium flow. Fig. 12 shows a comparison of the shock detachment distance 
between cylinder and sphere with the results for the sphere taken from Ref. 12. Note 
that the Mach number is nearly identical; however# the altitude is markedly different. 

Both are valid for zero free stream dissociation. It is seen that the shock wave for the 
sphere is much closer to the body and has less than one half of the distance of the shock 
wave for the cylinder. Results with exactly equal free stream conditions for both cases 
are not available at the present time. The fact that the shock wave of the sphere is 
much closer to the |j}ody than for the cylinder should be true also for identical free stream 
conditions; this is known to be the case for supersonic flow of a perfect gas. It has been 
found that the location of the sonic point is very sensitive# thus one cannot compare sonic 
point locations of two bodies with different free stream conditions. 

SUMMARY 

Analytical investigations of hypersonic nonequilibrium flow around blunt bodies are 
reported as contribution to the Interdisciplinary Symposium on Apollo Application Programs. 
Flight regions and equilibrium values of flow parameters behind the normal shock of a space 
vehicle# such as Apollo capsule# re-entering the Earth's atmosphere are presented. Basic 
features of equilibrium# nonequilibrium# and frozen flow are discussed. 

Real gas effects are treated using a simplified air model . With restriction to oxygen 
dissociation only# the thermal equation of state# the thermodynamic expressions of energy and 
enthalpy# and the rate equation for the net production of oxygen atoms are shown. Hypersonic 
inviscid flow about the forward portions of cylinders up to the equator has been numerically 
calculated in the range of oxygen dissociation# using the direct one-strip integral method of 
Dorodnitsyn. The investigations include the influence of the absolute body size and the 
effects of the free stream dissociation on the shock detachment distance and the flow field. 

For a cylinder flying at 30 km altitude# real gas effects can be detected beginning 
above M » 3 due to molecular vibration# and beyond M = 8 # due to molecular vibration and 
dissociation. Dissociation of oxygen only occurs in the range between Mach numbers 8 and 14# 
approximately# while at higher Mach numbers# nitrogen dissociation must be considered. The 
shock detachment distance for a real gas is much smaller than for a perfect gas. Free stream 
dissociation causes the bow shock to move away from the body. 

it 
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LIST OF SYMBOLS 


b 

c ' 


D 

D‘ 


e, E 


f: 


h 

h t 

k 

k 


M 

"l 



P 

q 


r, 0 


R 


R* 


T 


Constant l - ), def, by Eq. (1) 

Number of oxygen atoms per unit 
mass of the gas [ particles kg" 1 ] 

Dissociation energy [ J kmol * ] 

Characteristic temperature of 
dissociation [ °K 1 

Internal energy [ J kg * ] or 
[ J kmol - * 1 

Mole fractions [ - ] 

Source.function for oxygen atoms 
l sec - * 1 , Eq . (8) 

Enthalpy ( J kg”* ] 

Total enthalpy [ J kg”* ) 

Boltzmann's constant l J particle”* °K .*J ** 

Dissociation rate constant 
-1 -1 

( m° particle sec 1 ] 


u, v 

x* y 
z 

a 

Y 

L 


e 


Velocity components in x,ey 
direction, respectively, l m sec ] 

Coordinates ( m ] 

Compressibility factor l - 1 

Degree of oxygen dissociation ( - ] 

Ratio of specific heats l - ] 

Local shockwave detachment 
distance l m ] 

Dimensionless shockwave distance 
from body l - ] 

Characteristic vibrational tempera** 
ture l °K ], and body angular 
coordinate [ radians or degree ] 

-3 

Density [ kgm J 

Oblique shock angle [ radians or degree ] 


Concentration equilibrium 
constant { particles m”^ 1 

Mach number l - J 


Number of particles of ith 
species per unjt volume 
( particles m"^ J 


Avogadro's number 
[ particles kmol - * ] 

_2 

Pressure [ N m ] 

Velocity vector E m sec”* ] 


Radial and angular coordinates 
t m; radians or degree ] 

Gas constant of undissociated gas 
l J kg” T °K“' I 

Universal aas constant 

I J km«r’ °K-’ J 
Temperature [ °K ] 


Subscripts 

1 Free stream 
b Body surface 

s Behind the shock 

n Normal component 

t Tangential component, also 

total conditions 

N2 Molecular nitrogen 

02 Molecular oxygen 

O Atomic oxygen 

o Standard conditions, also 

stagnation streamline conditions 
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FIG. 1. HYPERSONIC FLOY/ REGIONS (SCHEMATIC) FOR APOLLO VEHICLE, REF. 1. 
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FIG. 5. DEGREE OF OXYGEN DISSOCIATION ALONG STAGNATION STREAMLINE FOR 
VARIOUS BODY RADII (FLIGHT SPEED 4300 m/sec, ALTITUDE 30 km, M, = 14.2] 
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FIG. 6, TEMPERATURE ALONG STAGNATION STREAMLINE FOR VARIOUS BODY RADII 
(FLIGHT SPEED 4300 m/s#c, ALTITUDE 30 km, M, * 14.2). 
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FIG. 11. TEMPERATURE DISTRIBUTION ALONG BODY SURFACE 
' (r. * 10 cm, ALTITUDE 30 km). 
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